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1.  liJTRODUCTIOil  m SUMMARY 


This  report  smimarizes  the  analytical  and  experimental 
investigations  of  the  infinite  array  radiation  and  coupling 
properties  of  bifurcated  twin  dielectric  slab  loaded  rec- 
tangular waveguide  dual  frequency  array  elements con- 
ducted under  Contract  P19628-75-C-0197 . Specifically, 
the  report  presents: 


•The  complete  analysis  of  the  element  in  infinite 
array  configurations. 

Theoretically^  determined  element/grid  design 
trade-off  conclusions,  leading  to  a proposed 
configuration  for  operation  over  15%  bands 
centered  at  4GHZ  and  8GHZ. 

.Experimental  verification  of  the  analytical 
results  and  the  examination  of  a unique  strip- 
line fed  notch  antenna  mode  exciter. 

.The  computational  details  and  computer  programs 
developed  during  the  study. 
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The  bifurcated  twin  dielectric  slab  loaded  rec- 
tangular waveguide  dual  frequency  array  element  shown 
in  Figure  1 is  a unique  concept  for  providing  simul- 
taneous aperture  usage  at  two  widely  separate  frequency 
bands.  At  low  frequency  both  upper  and  lower  halves 
of  the  waveguide  are  excited  in-phase  with  equiamplitude 
signals.  For  moderate  slab  loading  (assuming  relatively 
thin  slabs)  this  array  will  behave  similarly  to  a rec- 
tangular waveguide  array  excited  in  the  TE^^q  mode  with 
scan  behavior  associated  with  these  elements  in  the 
basic  lattice  (either  rectangular  or  triangular) . At 
the  high  frequency  the  first  odd  and  even  half-waveguide 
modes  can  be  independently  specified  such  that  four 
phase  centers  are  defined;  within  a single  low  frequency 
cell  the  fields  are  confined  predominately  to  the  slab 
regions . 

For  practical  array  designs,  the  low  and  high 
frequency  lattice  cells  are  identifical,  and  the  prin- 
cipal element/grid  design  trade-off  is  to  minimize  the 
number  of  phase  centers  (or  phase  shifters)  while  main- 
taining main  beam  purity  and  gain  over  a specified  scan 
volume,  particularly  in  the  high  frequency  band.  Con- 
sequently the  lattice  is  selected  such  that  element 
spacing  is  not  optimum  for  either  band,  but  presents 
the  best  compromise  of  element  count  versus  grating 
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lobe  free  scan  volume  in  the  high  frequency  band.  At 
scan  conditions  for  which  high  frequency  grating  lobes 
are  entering  or  present,  the  multimode  aperture  excita- 
tion is  adjusted  (slightly)  to  cancel  the  grating  lobe. 
These  considerations  are  treated  fully  in  Section  2 
which  gives  details  of  the  analysis  of  array  performance, 
and  section  4 which  summarizes  design  tradeoffs  and  ex- 
perimental results. 

The  application  of  the  twin  dielectric  slab 
loaded  element  to  dual  frequency  aperture  sharing  follows 
from  the  unique  propagation  properties  of  the  inhomo- 
generously  loaded  structure.  The  symmetric  loaded  guide, 
shown  in  Figure  2,  is  inherently  wideband.  At  sufficiently 
low  frequency,  a single  guide  mode  (the  mode)  propo- 

gates  and  has  an  electric  field  distribution  somewhat 
broader  than  the  homogeneously  loaded  guide  dis- 

tribution. As  frequency  is  increased,  the  LSE2Q  mode 
enters,  having  electric  field  distribution  similar  to 
the  TE2Q  distribution  of  the  homogeneously  loaded  guide. 
Concurrently,  the  LSE^^q  distribution  begins  to  develop 
a minimum  along  the  guide  mid-plane.  At  sufficiently 
high  frequency,  the  distributions  become  essentially 
identical,  except  for  symmetry  about  the  mid-plane,  and 
the  ratio  of  longitudinal  wavenumbers  approaches  unity. 
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By  appropriate  selection  of  guide,  slab  permittivity, 
and  operating  points  such  that  only  the  two  modes  propa- 
gate, the  guide  will  simultaneously  support  a single 
propagating  low  frequency  mode  with  phase  center  at  the 
midplane,  and  a conglomerate  high  frequency  distribution 
with  two  independent  phase  centers  at  (roughly)  the  slab 
centers.  The  analysis  of  propagation  in  the  inhomogenously 
loaded  guide  is  given  in  Section  3. 

In  general,  he  dispersion  in  the  inhomogeneously 
loaded  guide  is  not  linear  in  frequency.  Consequently,  the 
use  of  a bidirectional  exciter  requires  load  terminations 
at  the  back  of  the  guide  to  ensure  the  proper  aperture 
field  phase  at  both  frequencies,  and  results  in  a 3 dB 
power  loss.  This  difficulty  is  alleviated  by  a unique  uni- 
directional exciter  consisting  of  three  stripline  fed  flared 
(2  3 4) 

notch  antennas ' ' ' ' inserted  into  the  back  of  the  feed- 
guide  in  sue  h a manner  as  to  provide  a minimum  of  25  dB 
frequency  band  isolation.  Preliminary  experimental  investi- 
gation of  this  exciter  design  concept  was  begun  during  this 
study,  and  is  discussed  in  Section  5. 

Four  appendices  are  included.  Appendices  A,  B,  and  C 
give  the  details  of  the  analysis.  Expxicit  expansions  of 
modal  coupling  coefficient  integrals  are  given  in  Appendix  A. 
In  Appendix  B,  the  derivation  of  the  differential  equations 
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relating  feedguide  modal  fields  is  given.  And  in  Appendix 
C,  explicit  expressions  for  feedguide  mode  orthonormaliza- 
tion integrals  are  given.  The  remaining  appendix  gives 
complete  listings  of  all  programs  required  to  reproduce  the 
numerical  results  given  in  this  report. 

A time  dependence  e^^^  is  assumed  throughout. 


2.0  ANALYSIS  OF  INFINITE  PHASED  ARRAYS 
OF  BIFURCATED  TWIN  DIELECTRIC  SLAB  LOADED  RECTANGULAR 
WAVEGUIDE  DUAL  FREQUENCY  ELEMENTS 


In  this  section,  the  formal  solution  for  the  radiation 
properties  of  the  element  in  infinite  array  configuration 
is  presented.  The  unique  property  of  the  bifurcated  twin 
dielectric  slab  loaded  rectangular  waveguide  dual  freq- 
uency array  element  is  that  it  possesses  five  phase  centers: 
one,  associated  with  the  low  frequency  band  operation;  and 
the  remaining  four,  with  a high  frequency  band.  By  appropri 
ate  exciter  design,  the  element  can  simultaneously  operate 
over  both  bands. 

The  basic  element  is  shown  in  Figure  3.  Arrays  are 
formed  by  stacking  these  elements  in  rectangular  or  tri- 
angular grid  configuration.  The  element  consists  of  a 
rectangular  waveguide  bifurcated  in  the  E-plane  by  a 

septum  of  thickness  4.  Outer  dimensions  are  D and  D , ana 

X y 

inner  dimensions,  A and  B’,  where  D and  A are  associated 

X 

with  the  .X  coordinate.  Four  half  height  lossless  dielectric 
slabs  of  thickness  6 and  relative  dielectric  constant 

r 

are  located  at  distance  a + 6/2  {on  center)  from  the  narrow 


I: 


I 


Figure  3.  Dual  Frequency  Element 


9 


walls.  At  low  frequency,  a phase  center  is  maintained  at 
the  element  center  (i.e.,  over  the  septum)  by  exciting  the 
LSEj^q  mode  equally  in  the  two  half  guides.  At  high  freq- 
uency, four  independent  phase  centers,  at  roughly  the  four 
slab  centers,  are  formed  by  appropriately  exciting  the 
LSEj^q  and  1'SE2q  modes  in  each  half  guide. 

When  the  elements  are  arrayed  in  a rectangular  grid 
the  element  lattice  vectors,  s^i  and  £2  are  as  shown  in 
Figure  4a,  provided  the  septum  thickness, -6  , is  not  equal 
to(Dy  -B 'j/ft.  When  4 =^1^ -B'^l^^the  low  frequency  lattice  is 

as  shown  in  Figure  4a,  and  the  high  frequency  lattice  is 
as  shown  in  Figure  4b,  When  arrayed  in  a triangular 
grid,  the  lattice  vectors  are  defined  as  in  Figure  5, 
regardless  of  operating  frequency  or  septum  thickness. 

In  section  2.1,  the  formal  solution  for  active  array 
element  pattern  is  given.  The  method  of  solution  is 
similar  to  that  developed  by  Lewis,  ^ for  the 

analysis  of  a parallel  plate  array  with  protruding  dielec- 
tric. In  the  present  work,  the  formalism  is  extended 
to  two  dimensional  array,  cells,  and  the  unique  dual  freq- 
uency unit  cell  geometry  is  accounted  for.  In  section  2.2 
numerical  results  are  presented  and  particular  attention 
is  given  to  the  disposition  of  high  frequency  band  grating 
lobes.  Discussion  of  grating  lobe  suppression  is  deferred 


(o)Septum  Thickness, s=  Oy-B',  H}gh  ond  Low  Frequencies 


(b)Septum  Thickness  gs=Dy-B',  High  Frequency 


Figure  4.  Lattice  Definitions  for  Rectangular  Grids 
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Figure  5.  Lattice  Definitions  for  Triangular  Grid 
Either  Frequency  Band 
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to  section  4.1.  In  section  2.3,  numerical  results  are 
checked  against  published  data  for  several  limiting  geom- 
etries . 

2 . 1 Active  Array  Element  Pattern 

The  active  array  element  pattern  is  determined  from 
a unit  cell  formulation  of  scattering  at  the  feedguide 
free  space  interface.  The  interface  is  taken  as  coincident 
with  the  2 = 0 plane,  with  the  array  elements  occupying 
the  z < 0 half  space.  The  scattering  matrix,  S,  which 
relates  feedguide  modal  voltages  to  the  modal  voltages  of 
the  space  harmonics,  in  the  manner  indicated  by  the  network 
in  Figure  6,  is  obtained  by  matching  the  transverse-to-z 
fields  in  the  cell  across  the  interface.  The  field  matching 
is  accomplished  via  Galerkin's  method,  from  which  the  scat- 
tering formalism  follows  directly.  Active  array  trans- 
mission coefficient  is  then  obtained  from  the  network. 

In  the  following  discussion,  the  assumed  cell  configuration 
is  that  shown  in  Figure  4a.  The  extension  of  these 
results  to  either  of  the  other  two  cases  is  straight 
forward . 

For  the  configuration  of  Figure  4a,  the  unit  cell 
perimeter  may  be  taken  as  coincident  with  the  element 
outside  perimeter.  Thus,  the  unit  cell  overlays  two 
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Figure  6.  Network  Representation  of  Unit  Cell  Discontinuity 


independent  aperture  regions  and  the  modal  representation 
of  total  transverse- to-z  electric  and  magnetic  fields  at 
z = o”  is  given  as 


E^(s)  = U(y)p^^e>(s)+U(-y)p<jef  (s) 
H^(s)  = U(y)EIj^^h^(s)+U(-y)EI^jhJ(s) 


whe  re 


(3) 


U(£;)  = 


1,5>0 

O.C>0 


The  subscripts  > and  < are  used  to  distinguish  the  two 
aperture  regions.  V . and  I. . are  modal  voltage  and 
current  coefficients,  and  are  related  by 


(4) 


where  Y ^ is  the  modal  admittance  of  the  i^^  aperture 
< i 

mode,  and  the  single  ordering  index  i is  used  to  simplify 
notation.  The  mode  functions  e and  h are  given  in 
section  3.1. 
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At  z = 0 , the  unit  cell  guide  representation  of 


(11) 


k = ksin9^cos(j)_+pti  + qta„ 


(12) 


^ypq  = ksine^sind)^  + pti^  + qtz^ 


(13) 

t^*  Sj  = 2H6 

(14) 

c = 1 SiXSa 1 

is  Kronecker's  delta,  are  the  lattice  vectors, 

as  shown,  for  example,  in  Figures  4 a,b,  and  5,  and 

0^  and  (f>^  are  the  spatial  look  angles.  The  mode  function 

normalization  is  taken  such  that  V I*  is  power,  and 

apqr  apqr 

such  that  the  modal  voltage  and  current  are  related  by 


(15) 


I - Y V 
apqr  apqr  apqr 


where 


is  a modal  admittance,  given  as 


(16) 


apqr 


’"zpq'^o''o' 


1 

2 


k^  = 211/ A is  the  free  space  wavenumber  and  is  the 
free  space  impedance  377  ohms.  The  indices  pqr  and  pq, 
which  appear  explicitly  in  equations  ( 5)  through  (12) 
will  henceforth  be  replaced  by  the  single  index  a. 

Matching  transverse  fields  at  the  aperture  plane  of 
the  unit  cell  gives 
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(17) 


U(y)ZV^.eHs)+U(-y)ZV^  .e>(s) , 

£>11  j D 

0,  elsewhere  in  the  aperture 


and 


(18) 


ZI  h 
_ ao— aa 
0 


(s)=U(y)  (s)+U(-y)  ^I<.h<  (s) , 

i j 

in  the  aperture 


Approximate  solutions  for  the  parameters  of  the  network  in 
Figure  6 are  obtained  when  the  modal  series  are  truncated. 
As  a result  of  the  truncation,  the  continuity  equations 
(i.e.,  equations  (17)  and  (18)  can  no  longer  be  exactly 
satisfied,  and  vector  error  terms  Ai  and  A 2 must  be  inserted 
to  restore  the  equality.  These  error  terms  are  given  as 


(19) 


^ I J 

^ V,,,®=.n(s)-U(y)  ZV.  .eWs)+U(-y)S  .e<(s)  , 

0 ao— au  — 

in  the  apertures 

ZV  e ^(s),  elsewhere 
„ a0— a0  — 


and 


(20) 


M I J 

A2=na0ha0(s)-U(y)  p^£h>(s)  -U(-y)Zl^  ^h^  (s)  , 


in  the  apertures. 


Fisr^i- 


It  is  now  required  that  the  projections  of  Aj  and  A2  onto 

the  appropriate  modal  spaces  be  zero. 

*4" 

The  domain  of  definition  of  E^(s)  is  over  the  entire 
unit  cell,  and  the  domain  of  e”(£)  may  be  artificially 
extended  over  the  metalic  portions  of  the  cell.  Thus,  the 
domain  of  Ai  is  the  unit  cell,  and  the  modal  subset  span- 
ning the  space  are  the  h^^ ( s ) . Requiring  orthogonality  of 
Ai  to  the  ^(s)  and  performing  the  inner  products  over  the 
cell  results,  after  manipulation,  in 

(21)  = l\  * g\ 

> 

is  a matrix  of  coupling  coefficients,  the  elements 
are  given  as 

. = /dAef (s) • (h* (s)xz  ) 

(T , 1 y —1  — — — — o 

1/2  cell  ^ 

The  elements  of  the  column  vectors  are  the  feedguide 

< 

modal  voltages. 

The  domain  of  definition  of  A2  is  over  the  aperture 
only.  Therefore,  the  appropriate  basis  spemning  this  space 
is  formed  by  the  concatination  of  the  trunicated  modal 

— _ > 

'^Complete  expressions  for  E . are  given  in  Appendix  A 

0 f 2. 


where  | 
of  which 

(22*) 
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(26) 


> > 

= E^* 

i,a  0,n 


It  is  convenient,  then,  to  define  the  partitioned  vector 
V such  that 


(27) 


and  the  partitioned  matrix  E as 


(28)  £ = (£^li^) 

Then,  the  voltage  and  current  equations,  (3-21)  and  (3-24), 
respectively,  take  the  form 

(29)  V,  = E V 

— ”” 

(30)  I = |*% 

where  the  asterisk  denotes  conjugation  and  the  t denotes  the 
transpose  operation. 

The  vectors  V,  , I,  and  in  equations  (29)  and 
(30)  are  total  modal  voltages  and  currents.  Using  the 
conventions  established  for  the  network  in  Figure  6, 
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assuming  that  all  external  sources  are  zero  (i.e.  V^  = 0) , 
and  manipulating  equations  (29)  and  (30)  results  in 
an  expression  for  feedguide  reflected  field  voltage  coef- 
ficients, V , in  terms  of  the  active  modal  excitations, 

V , gxven  as 


(31)  v“  = {2[Y  + i*^E]“^Y-l}v‘*' 


where  is  the  identity  matrix. 

Let  the  scattering  matrix  of  the  network  be  defined 
by 


(32) 


-’R"  ’1 

\S  2 1 S 2 if 


where  the  scattering  blocks  have  the  usual  meaning.  Then, 
from  equation  ( 31) 


(33) 


S.i*  2[Y  + E*^  y n"^Y  -1 


and  from  equation  ( 29) 


(34) 


%zT  + Sii) 
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with  V 
—a 


0. 


In  more  standard  array  configurations,  the  feedguide 
and  aperture  are  designed  for  single  feedguide  mode  propa- 
gation in  a single  frequency  band.  In  these  configurations, 
the  active  array  reflection  coefficient  is  simply  the  one 


f 


element  of  Sj i corresponding  to  reflections  in  the  domin- 
ant mode.  Defining  that  matrix  element  as  T(6o  <j>o)/  norm- 

9 

alized  active  array  element  gain  pattern  is 

? 

(35)  gg(6o/4>o)  = (1-  lr(6o,(J)o)|)  COS0O 


f- 


provided  no  grating  lobes  have  entered  real  space.  Follow- 
ing the  entrance  of  the  first  grating  lobe,  it  is  necessary 
to  track  the  propagating  beams  individually,  and  the 
relative  power  in  the  beam  due  to  the  single  excited 
mode  (i=l)  is 


(36)  | 


where  T^(0,(J))  is  the  (a,i)th  element  of  the  partitioned 
block  S21,  and  0 and  ({>  are  the  actual  location  angles  of 


Active  array  element  pattern  is  defined  in  a different 

manner  for  the  array  of  bifurcated  twin  dielectric  slab  loaded 

rectangular  waveguides.  In  this  instance,  the  aperture  is 

always  considered  to  be  multimode.  In  the  low  frequency 

til 

range,  the  mode  (ordered  as  the  i ) is  excited 

equally,  in  amplitude  and  phase,  in  both  regions  of  the  aper- 
ture. Consequently,  for  an  I-mode  feedguide  aperture  field 
approximation  in  both  upper  and  lower  regions,  the  power  in 
the  main  beam,  a=m,  is  given,  for  principle  plane  scan*, 
as 

(37)  Pjj^(6o,4)o)  = l/2|s^{^(6o,<)>o)  + ( 0 o , 4>  o ) 1 " 

where  S^i^(9,<{))  is  the  (m,i)th  element  of  £21,  and  the 
r feedguide  mode  is  the  propagating  mode.  By  the 

definitions  in  equations  (27),  (32),  and  (34), 

(0o/H>o)  is  the  voltage  transmission  coefficient  for  coupling 
from  the  i^  mode  in  the  upper  aperture  region  to  the  m 


Equation  (3-37)  is  strictly  va.1id  only  for  6o5^0  and 
in  the  principle  planes  ((()=0,II,  or  <p=Jl/2,  3IT/2)  . For  all 
other  scan  planes  (and  at  broadside) , the  total  power  in  the 
main  beam  is  the  sum  of  the  powers  in  the  dominant  (p=q=o)E 
and  H modes. 
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beam  in  free  space;  and  is  the  voltage  transmission 

from  the  i mode  of  the  lower  aperture  region  to  the  m^ 
beam* . 

Active  array  reflection  coefficient  is  also  obtained 
via  superposition.  For  the  multimode  aperture  configuration, 
it  is  necessary  to  independently  track  all  propagating  waves 
in  the  feedguide.  Hence,  for  low  frequency  excitation, 
the  total  reflected  power  in  the  mode  of  the  upper 

aperture  region  is  given,  for  any  scan  angle,  as 

(38)  R^(6o/4>o)  = l/2|S^'^(6o,4>o)  + («  o ,<!>  o ) | 

For  the  lower  aperture  region,  the  reflected  power  is 

(39)  R^(eo,<()o)  = l/2|stt^'^(0or(})o)  + ^ (0  o , <J>  o ) | ^ 

X"  t 

where  SiJ  (9o/<J>o)  is  the  voltage  scattering  coefficient 
from  the  r aperture  region  mode  to  the  t^  aperture  region  mode. 
It  is  evident  from  these  two  equations  that  the  reflected 
powers  in  the  two  regions  are  not  necessarily  equal. 

Indeed,  it  is  found  that  for  low  frequency  excitation, 

R^(6o/<l>o)  equals  R*^(0Of<J>o)  only  in  the  H plane  of  scan. 


*It  is  assumed  that  the  mode  ordering  in  the  two  aperture 
regions  is  the  same.  This  assumption  will  carry  through 
the  remainder  of  the  report. 


As  an  example,  reflected  power  is  shown  versus  E plane 
scan  angle  in  Figure  7.  The  element  is  a WR187  guide  with 
.250"  thick  slabs  of  £^=9  dielectric  located  .450",  on  center, 

from  either  narrow  wall  and  with  a .032"  septum. 

The  operating  frequency  is  2.5  GHz.  There  is  considerable 
difference  in  reflected  power  between  upper  and  lower 
regions  throughout  the  scan  range  .17  <sin6<.95.  Consequently < 
low  frequency  excitation  of  the  upper  and  lower  regions  of 
the  element  from  a common  post  phase  shifter  feed  point,  as 
is  desireable  tor  several  low  frequency  feed  concepts,  will 
produce  an  imbalance  at  the  outputs  of  the  power  divider  net- 
work. Since  the  impact  of  this  effect  on  feed  and  exciter 
design  is  beyond  the  scope  of  this  study,  it  will  be  given 
no  further  consideration  in  this  report. 

At  the  high  frequency  band,  the  element  is  excited 
such  that  four  independently  controllable  phase  centers  are 
distributed  in  the  aperture.  To  maintain  this  phase  center 
distribution,  four  propagating  modes,  two  in  each  region, 
are  excited.  The  modes  are  LSEj^q  and  LSE2Q.  For  sufficiently 
high  frequency  and  dielectric  constants,  these  two  modes  have 
nearly  equal  dispersion.  In  addition,  to  a crude  approxi- 
mation, the  modal  field  distributions,  e";  and  e" 


differ 


only  in  symmetry  and  behave  like  | s i n ( 2 l!x/A}  j and  sin(2IIx/A), 
respectively.  By  exciting  the  bSE^p  and  bSE2Q  modes  of  the 
upper  aperture  region  with  voltages 

(40)  Vi  = Vi  = cos  ( . 25k^D^sin0^cos(J)^) 
and 

(41)  V2=Va  = jsin(.25k  D sinG  cos(p  ) 

0X00 

respectively,  and  the  modes  of  the  lower  aperture  region 
with  voltages 

(42)  Va  = Vi  = vf  exp[j  0.5)^  B'^in0^sin(p^] 
and 

(43)  V4=  Va  = Va  exp  [ j 0 . 51^  B feinG^sinc))^] 

the  four  phase  centers  are  established  a^  roughly,  x = +A/4 

in  each  region.  The  beam  is  scanned  to  (0  (}>  ) . 

Of  o 
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As  for  low  frequency,  the  normalized  power  in  the 
high  frequency  propagating  beams  is  determined  via  super- 
position (i.e.,  using  equation  (29)).  The  total  power 
delivered  is 


2 2 

P = 2[Yi  |V^  + YaiVal  ]=  2[Yj+Y2] 


The  power  in  the  beam  is  therefore,  given  as 


P^(9  ,<l>  ) = 
a a a 


where  the  ordering  of  the  elements  of  V has  been  altered 
to  simplify  the  equation.  Taking  the  same  liberty  with 
mode  ordering,  the  power  reflected  in  the  jth  mode 
(j  = 1,2, 3, 4)  is  given  as 


R ■ (9  ,<P  ) = 

J 0 ^0 


2 S^I^(e^,4  )Vj  Y./P 


O "O  1 ' J 


2 . 2 Numerical  Results 

In  this  section,  numerical  results  are  presented 
for  several  element  and  grid  geometries  to  illustrate  the 
principle  performance  characterio tics»  uf  tue  bifurcateu 


*See  footnote  to  equation  (37) 


twin  dielectric  slab  loaded  rectangular  waveguide  dual  freq- 
uency array  element.  Element/grid  design  is  discussed  more 
fully  in  section  4. 

For  purpose  of  discussion,  it  is  convenient  to  present 
performance  data  in  a somewhat  unusual  format.  Rather  than 
present  realized  gain  pattern,  (i.e.,  normalized  directive 
gain)  power  transmission  coefficient  is  given  for  each 
radiating  beam.  The  advantage  gained  by  this  form  of  pres- 
entation is  that  it  allows  a direct  comparison  of  the  power 
in  the  radiated  beams.  If  P^(0^,<j)^)  is  the  power  associated 
with  the  beam  when  the  main  beam  is  scanned  to  (0o/<j>o)/ 
then  the  directive  gain  of  this  beam  is  proportional  to 
P^j(0^^/<j)|j)cos0^.  Consequently,  for  a given  scan  angle,  com- 
parison of  beam  directive  gains  includes  the  comparison  of 
projected  aperture  at  the  various  beam  locations  . 

Data  are  presented  for  rectangular  and  triangular  grid 
configurations.  Two  elements  are  discussed;  a WR187  guide 
with  .250",  e^=9  slab  loading;  and  a WR137  guide  with  .125", 
4.75  slab  loading.  The  former  is  operated  at  2.5  GHz 
(k_A/2=1.246)  and  6.0  GHz  (kA/2=2.990)  . The  later  is  operated 
at  4 GHz  (kA/2=1.463)  anvd  8 GHz  (kA/2=2.926)  . The  dimen- 
sions of  the  elements  are  given  in  Table  1.  In  the 

following,  the  elements  will  be  distinguished  by  the  WR 
number  of  the  rectangular  guides. 
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WR187  WR137 

£^=9  e^=  4.75 


2.000 

1.000 

1.872 

.872 

.032 

.420 

.325 

.361 

.250 


1.500 

.750  or  .960 
1.374 

.622  or  .832 
.032 

.295  or  .400 
.281 
.281 
. 125 


Table  l 

Dimensions  (in.)  of  WR187  and  WR137  Elements 


The  grating  lobe  diagrams  for  the  four  grids  are  shown 
in  Figures  8a,b,c,d.  Each  diagram  shows  the  near  in 
grating  lobe  locations  for  the  two  operating  points.  Low 
frequency  grating  lobes  are  indicated  by  solid  boxes, 
and  high  frequency  lobes,  by  solid  dots.  What  is  immediately 
obvious  from  the  figures  is  that  the  triangular  grid  provides 
a grating  lobe  free  scan  region  for  all  directions  in  the 
plane  at  high  frequency.  One  consequence  of  this  is  improved 
high  frequency  broadside  match,  cis  is  shown  below. 

Figures  9 through  12  show  power  transmission 
coefficientin  the  principle  planes  at  low  frequency  for 
each  of  the  four  grids.  In  the  scan  range  sin6o<.95,  no 
grating  lobes  enter,  as  can  be  seen  from  the  grating  lobe 
diagrams  in  Figure  8.  In  general,  the  performance  of  the 
four  configurations  is  the  same.  There  is  some  improvement 
in  scan  coverage  of  the  WR137  element  over  the  WR187  element, 
but  the  difference  is  not  large  enough  to  show  up  on  the 
scale  of  the  figures.  In  the  E-plane  (sin4i=  1.0),  the 
fall  off  is  nearly  cos^'^^O,  out  to  60°  - the  WR187  element 
shows  slightly  greater  scan  loss.  In  the  H-plane,  the  scan 
loss  exceeds  cos  '^^6  by  approximately  .5  db  at  9=60°.  As 
might  be  expected,  the  performance  of  the  rectangular  and 
triangular  grid  configurations,  as  measured  by  power  trans- 
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Figure  9.  WR.187  Element  Power  Tr;in«;mi  ssion  Coefficient 
Rectangular  Grid,  Low  Frequency  (=2.5GHz) 
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Figure  12,  WR137  Element  Power  Transmission  Coefficient 
Triangular  Grid,  Low  Frequency  (=4  GHz) 


mission  coefficient,  is  not  distinguishable,  one  grid  from 
the  other. 

In  Figures  8a  and  8c,  the  high  frequency  grating 
lobes  of  the  rectangular  grid  configurations  are  shovm  re- 
siding in  real  space  for  the  no  scan  condition.  In  this 
situation,  a broadside  power  loss  occurs,  and  a high  farout 
sidelobe  condition  is  created.  To  a certain  extent,  the 
power  delivered  to  these  lobes  can  be  reduced  by  introducing 
a complex,  multiplicative  correction  for  the  voltage  excita- 
tion coefficients,  Vz  and  V4,  given  by  equations  (41) 
and  (43)  f Ideally,  this  correction  is  independent  of 
scan  and  frequency  (for  small  bandwidths)  . Typical 
grating  lobe  levels,  with  and  without  the  correction,  are 
shown  in  Figure  13  for  a thin  walled  element  in  a rec- 
tangular grid.  Without  correction,  grating  lobe  levels  reach 
-16.7  db.  With  the  multiplier  1.164  - j.291*,  the  maximum 
grating  lobe  level  is  -20.4  db,  and  beyond  sin0o=.3,  the 
level  is  below  -25  db.  Differences  in  the  main  beam  due  to 
the  two  excitations  are  too  small  to  be  represented  in  the 
figure . 


*The  case  shown  here  is  supplied  by  R.  Mailloux,  the  multi- 
plier 1.164  - j.291  was  also  found  suitable  for  the 
WR187  elemertt  in  a rectangular  grid. 
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By  configuring  the  WR187  and  WR137  elements  in  triangular 
lattice,  the  difficulties  encountered  due  to  the  location 
of  the  high  frequency  grating  lobes  are  largely  eliminated. 
Figures  8b  and  8d,  show  the  disposition  of  high  freq- 
uency grating  lobes  for  triangular  grid.  No  grating  lobes 
enter  real  space  along  any  cut  plane  for  sin9o<.391. 

In  the  principle  cut  planes,  only  the  main  beam  is  propa- 
gating out  to  sin6o=.820,  or  nearly  60^  scan. 

For  wide  angle  scan  applications,  the  triangular  grid 
configuration  does  not  fully  eliminate  the  necessity  for 
the  type  of  excitation  correction  described  above.  Consider 
Figures  8c  and  8d.  The  (1,0)  harmonics  of  the  structures 
follow  the  same  track  for  E plane  scan,  and  while  they  are 
in  real  space,  one  is  the  image  of  the  other,  reflected 
about  the  /k^  axis.  Hence,  for  the  same  voltage  correc- 
tion term,  loughly  the  same  amount  of  power  will  be  dumped 
in  the  triangular  and  rectangular  grid  lobes  when  the 
E plane  scan  sines,  sin0^  and  sin6^,  respectively  are 
related  by 

(47)  sine^  = A/Dy  - sin9^ 

This  is  illustrated  in  Figure  14  for  the  WR137  element 
operating  at  8.64  GHz.  Both  E and  H mode  space  harmonics 
are  shown.  It  is  clear  from  the  high  power  levels  shown  in 
the  figure  that  grating  lobe  control  is  required  for  either 
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8.64  GHz  A = 1 .374  in. 

0.278  in.  B = 0.295  in. 

0.278  in.  = 1.500  in. 

0.132  in.  Dy  = 0.750  in. 

4.75  in.  SEP  = 0.032  in. 

CUT  PLANE  SIN^=  1 .000 


MAIN  BEAM  - TRIANGULAR  GRID 

MAIN  BEAM  - RECTANGULAR  GRID 

1 MODE  - TRIANGULAR  GRID 

1 MODE  - RECTANGULAR  GRID 


Figure  14.  Grating  Lobe  Levels  Obtained  from  Rectanguiar  and 
Triangular  Grid  Configurations  of  WR137  Elements 
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grid  type,  even  though  the  cos6^  beam  broadening  factor  is 
less  than  .4  throughout  the  scan  range. 

Because  of  the  high  grating  lobe  levels  obtained 

for  the  rectangular  grids,  there  is  significant  (approximately 

.1  db)  main  beam  loss  at  broadside  relative  to  main  beam 

levels  obtained  for  the  triangular  grid.  For  the  WR137 

element,  the  difference  is  .13  db.  A comparison  of  principle 

plane  main  beam  levels  for  the  two  grids  is  shown  in  Figures 

15  and  16  . The  operating  point  is  8 GHz.  In  the  H 

plane.  Figure  15  , the  beams  show  the  .13  db  difference  at 

broadside,  and  smoothly  coalesce.  In  the  E plane.  Figure 

16  , there  are  sharp  jogs  in  the  curves.  For  the  rectangular 

grid,  the  jog  occurs  as  the  grating  lobe  pair  exits  from 

real  space.  For  the  triangular  grid,  the  sudden  power  loss 

at  sin0o  = .82  is  due  to  the  grating  lobe  pair  entering 

real  space.  The  main  beam  falloff  is  slightly  greater  than 
1/2 

cos  ' 00  for  both  grids  in  either  plane. 

2.3  Comparison  with  Limiting  Cases 

To  check  the  analysis,  and  in  particular  the  details 
of  the  computational  procedure,  several  limiting  cases  have 
been  examined.  For  these  cases,  the  slab  dielectric  con- 
stant is  set  to  ej.=l  and  the  LSE^^q  mode  is  independently 
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sin  0 


Figure  15.  Comparison  of  Main  Beam  Power  Levels  for  Rectangular  and 
Triangular  Grid  Configurations  of  WR137  Elements  - H Plane 


Figure  16.  Comparison  of  Main  Beam  Power  Levels  for  Rectangular  and 
Triangular  Grid  Configurations  of  WR137  Elements  - E Plane 


excited  in  the  upper  and  lower  aperture  regions. 

For  e =1,  the  LSE  modes  degenerate  to  the  TE  ^ modes 
r no  ^ no 

of  empty  rectangular  waveguide.  Therefore,  an  appropriate 
set  of  check  cases  include  H plane  scanned  thin  wall  rectangular 
grid  arrays  of  rectangular  elements  and  special  triangular 
grid  examples  which  have  appeared  in  the  literature.  Thick 
wall  rectangular  grid  cases  can  also  be  used  as  checks 
provided  the  wall  thickness  is  not  too  great. 

The  geometry  for  the  rectangular  grid  examples  is 
shown  in  Figure  17.  The  aperture  dimensions  are  A and  B, 
and  the  lattice  vectors  are 

(48)  Si  = d X 

' ' — ‘ x-o 

and 

(49)  S.2  = dyi^ 

For  A=d  and  B«d  , exact  solutions  for  E and  H plane 

X y 

element  patterns  have  been  obtained  using  function  theoretic 

techniques  to  construct  the  reflection  coefficient  of  the 

(8) 

driven  TE^^q  mode.  Figures  18  through  20  show  magnitude  and 

phase  of  active  reflection  coefficient,  F,  of  the  TE2.0 

square  waveguide  in  thin-wall  square  lattice  configuration  for  H 
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Figure  18.  Comparison  of  Exact  (8)  and  Approximate  Modal  Solutions 
for  Active  TEj^q  Reflection  Coefficients  - Thin  Walled 

Square  Elements  H-Plane  D .5714 


Figure  19.  Comparison  of  Exact and  Approximate  Modal  Solutions 
for  Active  TE^q  Reflection  Coefficient  - Thin  Walled 

Square  Elements  H-Plane  .6205 


Figure  20.  Comparison  of  Exact^®^  and  Approximate  Modal  Solutions 
For  Active  TE^q  Reflection  Coefficient  - Thin  Walled 

Square  Elements  H-Plane  DvA  = .6724 


plane  scan.  The  solid  curves  are  the  exact  solution^  as 

(8) 

obtained  by  Wu  and  Galindo.  The  dashed  curves  are  obtained 

from  the  current  formalism  using  the  first  five  TE  ^ modes 

no 

to  approximate  the  aperture  field  distrubution.  The  minor 

discrepancies  between  the  results  are  removable  by  including 

additional  higher  order  modes  in  the  modal  computations. 

Figures  21  through  23  show  magnitude  of  T as  a 

function  of  H plane  scan  angle  when  the  H plane  metalic  walls 

have  finite  thickness,  t = • ^“^x*  lattice  remains  square. 

The  solid  curves  were  obtained  by  Galindo  and  Wu  using 

30  feedguide  modes  to  represent  the  aperture  field.  The 

dashed  curves  are  from  the  present  analysis  using  the  first 

nine  TE„_  modes.  The  agreement  is  excellent, 
no 

i^erture  field  approximations  consisting  of  the  first 
few  TE^^  mode  functions  may  be  used  to  compute  H plane 
element  patterns  of  triangular  grid  arrays  of  rectangular 
apertures. 

Figures  24  and  25  show  H plane  element  pattern, 

, , .2  (10) 

Cl  - |r|  ]cos9,  as  computed  by  Amitay,  Galindo,  and  Wu 

and  using  the  approximate  limiting  case  of  the  present  analysis 


Figure  22. , Comparison  of  Modal  Solutions  for  Active  TE 

Reflection  Coefficient  r-  Thick  Walled  Square  Elements 
H-Plane  DxA  = .6205,  t = .1D„ 


Figure  23.  Comparison  of  Modal  Solutions  for  Active  TE 

Reflection  Coefficient  >-  Thick  Walled  Square  Elements 
H-Plane  Dx/\  = .6724,  t = ,1D. 


Figure  24.  Comparison  of  Modal  Solutions  for  H-Plane  Element 
Pattern  When  Aperture  Field  is  Approximated  by  TEj^q  Mode 
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The  lattice  is  45  triangular,  with  lattice  vectors 


Si  = i.ooaxx 

— * — o 


Sz  = .504X^  + .504X^ 


The  rectangular  apertures  are  .90BX  by  .4i(..  In  figure 
24,  the  aperture  field  is  approximated  using  only  the 
TEj^q  feedguide  mode.  For  Figure  25  , the  TE2Q  mode 
is  added.  Again,  the  comparison  is  quite  good.  The  resonance 
near  sin9=.6,  was  found  experimentally  by  Diamond^^^^  in  an 
investigation  of  the  central  element  pattern  of  a 95  element 
array . 

The  deviation  of  computed  results  in  Figures  24  and 
25  arise  from  two  sources.  The  first  is  an  inability 
to  read  the  published  curves  to  sufficient  accuracy.  The 
second  source  of  error  is  number  of  space  harmonics  used  to 
obtain  the  published  curves  and  the  present  results. 

2.4  Convergence  of  Numerical  Results 


From  the  results  of  the  previous  section,  it  is  event  that 
the  nxamerical  solution  implemented  here  converges  uniformly  as 
the  numer  of  aperture  modes,  space  modes -or  both  is  increased. 
However,  to  ensure  that  the  convergence  is  indeed  uniform,  the 
array  properties  of  the  WR187  element  were  examined  as  the  mode 


count  was  varied  in  the  aperture  and  free  space  regions. 

One  to  eight  aperture  modes  were  considered  and  the  circle 

of  convergence  was  varied  from  (=  2tt/X)  to  13. 9K^  at 

2.5GHZ  and  33.4k  at  6 GHZ. 

o 

For  maximum  circle  of  convergence,  the  solution 
converged  rapidly  as  the  number  of  aperture  modes  was 
increased.  At  both  frequencies,  four  or  five  aperture 
modes  were  found  to  be  sufficient. 

With  the  number  of  aperture  modes  held  fixed,  the 
circle  of  convergence  was  uniformly  increased.  Beyond 
roughly  5k^  at  2.5GHZ  and  lOk^  at  6GHZ,  the  solution  became 
stable. 


Except  for  the  usual  effect  of  truncating  the  modal 
series  at  an  inordinately  premature  point,  no  convergence 
anomalies  were  evident  in  the  computations. 


3,0  PROPAGATION  CHARACTERISTICS  CF  TWIN  DIELECTRIC 
SLAB  LOADED  RECTANGULAR  WAVEGUIDE 


The  analysis  of  the  radiation  properties  of  an  infinite 
array  of  dual  band  elements  requires  a complete  description 
of  wave  propagation  in  the  inhomogeneous ly  loaded  feedguide  . 
The  geometry  of  the  dual  band  element  is  shown  in  Figure  26. 
The  element  consists  of  a rectangular  waveguide  bifurcated 
in  the  E-plane  by  a septum  of  thickness  ,4  ,and  symmetrically 
loaded  with  full  height  slabs  of  lossless  dielectric  parallel 
to  the  guide  narrow  wall.  Both  slab  thickness  and  spacing 
are  arbitrary  in  the  ranges  0<  2&,/A<l,  0<2g/A<l,  and  c^,  the 
relative  dielectric  constant  may  take  on  any  real  value, 
e >1. 

^ (12-14) 

Several  investigators  have  studied  wave  propaga- 

tion in  similar  guiding  structures,  primarily  to  provide 

bases  for  perturbation  calculations  of  ferrite  phase  shifter 

(15) 

properties.  Collin  has  obtained  general  expressions  for 

mode  functions  and  modal  propagation  constants  for  the 

asymmetric  single  slab  case.  These  results  are  equivalent 

to  the  symmetric  twin  slab  results  for  short-circuit  sym- 

(12) 

me try  in  x.  Seckelmann  has  obtained  general  expressions 
for  LSE  „ (i.e.,  TE  ) mode  functions  and  propagation 


W*5W," 


•UJ^  H 


constants  of  the  symmetric  twin  slab  case. 

The  modal  fields  of  the  inhomo^eneously  loaded  rec- 
tangular waveguide  shown  in  Figure  26  obtained  in  a 

straightforward  manner.  Recognizing  that  E and  H modes  with 
respect  to  remain  decoupled  at  the  dielectric  interfaces, 
arbitrary  waveguide  fields  can  be  decomposed  into  E-type 
(LSM)  and  H-type  (LSE)  modes  with  respect  to  Thus,  an 

arbitrary  field  may  be  expressed  in  terms  of  either  complete 
node  set.  A one-to-one  correspondence  exists  between  the 
modes  of  each  set  (i.e.,  eigenvalues  of  the  E ind  H mode 
set  are  eigenvalues  of  the  type  mode  set  for  the  given 
boundary  value  problem) . The  LSE  and  LSM  modal  fields  are 
then  obtained  via  a component-by -component  comparison  of 
the  modal  fields  in  either  set  corresponding  to  a particular 
eigenvalue.  The  components  of  the  type  mode  functions  are 
then  proportional  to  either  the  voltage  or  current  distri- 
butions of  the  equivalent  circuit. 

The  modal  spectrum  for  the  structure  is  obtained  via 
a transverse  resonance  technique. 

3 . i Mode  Functions 

It  is  well  known  that  E and  H modes  with  respect  to 
surface  normals  remain  decoupled  at  planar  interfaces 


60 


Pljjij  ijij,  1^1.4 1-'  1 


(15) 

between  dielectrics.  Thus,  in  the  inhomogene ously 
(in  x)  loaded  guide  shown  in  Figure  26  , decoupled  modes  of 
the  structure  will  be  E and  H with  respect  to  x,  or, 
equivalently.  E-type  (LSM)  and  H-type  (LSE)  with  respect  to 

z . 

For  the  infinite  phased  array  of  dual  frequency  elements, 
the  aperture  fields  are  expressed  as  the  superposition  of 
transverse-to-z  mode  functions  which  satisfy  the  vector 
equations  (see  Appendix b )• 


( 52  ) 
(.53  ) 


- ^t'^t 

(x,y)  = o)c[I  + — 


V V 

h Y;Z"e".  (x,y)  = coy  [I  + -] 

-o  i.  yi.  ,2 


(z  X e ' (x,y) ) 

"*Q  ““i 

(hj|(x,y)  X z^) 


In  equations  ( 52  ) and  ( 53  ) / y is  the  longitudinal  (z 
directed)  wavenumbers;  Y'  and  Z"  are  modal  Immittances;  I 
is  the  unit  transverse-to-Z  diadic, 

(54  ) I = 2o2o 


and  is  the  tranverse-to-z  gradient  operator. 
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'55)  ’t=£o^ 

The  prime  (')  is  used  to  denote  LSM  modes  [for  which 
h (x,y)  s 0];  the  double  prime,  to  denote  LSE  modes 

X 

[for  which  e (x,y)  s 0] ; and  single  index  is  used  rather 

X 

than  a double  index. 

The  desired  modal  representation  of  transverse  fields 
is 

(56)  E.  (x,y,z)  = Z V!  (z)e!(x,y)  + Z V|!  ( z)e  V (x,y) 

^ i i -1  j J J 

(57)  H.(x,y,z)  = Z I!(z)h!(x,y)  + Z IV  (z)h!!  (x,y) 

“t  ^ i 1 j J 3 

where  (z)  and  l^(z)  (ot  = ',")  are  z dependent  modal 
voltages  and  currents  satisfying  the  transmission  line 
equations 

(58)  liVi(z)  = -j7iZiIi(z) 

(59)  |j  1.(2)  = -jY.Y.V.{2) 

Since  the  guide  is  uniform  (and  assumed  infinite)  in  z, 
the  z dependence  of  the  modal  voltages  and  currents  is 
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exp  [-jYj^z],  hence 


( 60  ) V^(z)  = 


(61  ) (z)  = I^e 


-jYz 


(i 

r 


i 


I 

A 


Equations  ( 52  ) and  ( 53  ) may  be  solved  to  obtain  re- 
lationships between  the  mode  components.  For  LSM  modes, 
the  X component  of  h’ (x,y)  is  taken  as  zero  (h^  = 0) , and 
equation  ( 52  ) results  in: 


{ 62  ) h{(x,y)  = Gj^Cx) 


( 63  ) e! (x,y)  = e’ . (x,y)x 

”“X 


— o 


k"e^(x)-ic?  (X) 


9x 


^ e;i(x,y)i:^ 


(64  ) 


Z! 

1 


ke^(x)-K?  (x) 

Y !o)e 
1 o 


where  e' . (x,y)  is  a solution  of  the  scalar  wave  equation 
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V 65  ) 
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(Vj 


kJi(x)) 


e;i(x,y)  = 


with  = k|(x)  + = k^e^(x)  - y|/  subject  to  the 

boundary  conditions  of  the  guide  cross-section,  ^^.(x)  is 
the  x-dependent  dielectric  constant  of  the  cross-section 
For  LSE  modes,  the  x component  of  e"{x,y)  is  zero 
(e"  s 0)  , and  solution  of  (53)  gives: 

( 66  ) eV(x,y)  = 


( 67  ) hV  (x,y)  = h^  (x,y)x^ 


( 68 


■c^(x) 


(X) 


^ h"  . (x,y)  V 
a x:  9y  XI  ' ' -k) 


Y 

X 


2 

k£  (X)-K?(X) 

r 1 

topY  . 


with  h". (x,y)  satisfying 
( 69  ) (V;  + k^  (X)  h". (x,y)  = 0 


over  the  cross-section.  The  wave  immittances  given  by 
equations  (64)  and  (6  8)  are  defined  such  that  the 

direct  proportionalities  of  equations  (62)  and  (66) 
are  obtained.  It  is  shown  in  Appendix  C that  the  type 
modes  possess  the  following  orthonormality  property  for 
the  bounded  cross-section  (CS)  of  Figure  26. 

(70)  • (hr^'  2o>  “ 'a6«nm 

CS 

where  a, 8 = (',"). 

Equations  (62)  through  (69)  r with  appropriate 
boundary  conditions,  are  the  complete  formal  solution  for 
the  mode  functions  of  the  symmetric  twin  dielectric  slab 
loaded  rectangular  guide  shown  in  Figure  26.  However, 
due  to  the  complexity  of  boundary  conditions  along  x,  the 
scalar  wave  equations  (65)  and  (69)  are  difficult  to 
solve.  If  the  transmission  line  direction  is  temporarily 
taken  along  x,  the  fields  in  the  guide  may  be  put  in  a 
representation  of  E and  H modes  with  respect  to  x.  Due  to 
the  degeneracy  of  the  rectangular  cross-section,  eigenvalues 
of  the  E(H)  mode  set  are  also  eigenvalues  of  the  LSM  (LSE) 
mode  set.  Thus,  corresponding  to  each  eigenvalue  of  the 


structure,  there  are  two  expressions  for  total  field. 
These  expressions  are  compared  component-by~component, 
resulting  in  particular  expressions  for  LSE  and  LSM 
mode  functions  in  the  symmetric  twin  dielectric  slab 
loaded  rectangular  waveguide. 

The  transverse-to-x  modes  of  the  twin  slab  loaded 
guide  are  obtained  in  standard  fashion  and  are  given  as: 
E modes  (h'  = 0) 

X 

4 ^t’^E 

(71)  e.  (y/2)  = 


Ae 

— 


-jY^z 


.mil 


cos 


Zo-JVj^sln 


2^  j 

b ~o 


ti 


(72)  hj^(y/Z)  = x^xe|(y,z) 

K^(x) 
we„e„(x) 


(73) 


z: 

1 


H Modes  (e"  e 0) 

X 


h!’(y,z)  = - 


i_  .in  ^ ^ . iViCos!^  2„I 


eV(y,z)  = hV(y,z)  x ^ 


'^i(x) 


Y!<  = ^ 

o))j 


In  the  above  equations,  is  used  to  indicate  results  in  the 
transverse-to-x  representation,  and  = /(^)  + y ? . 

The  modal  representation  of  tranverse-to-x  fields  is 


E^(x,y,z)  = LV!  (x)e.  (y  ,z)  + J:v::  (x)e;:  (y  ,z) 

t ^11  j j j 


H+.(x,y,z)  = Zl.  (x)n!  (y,  z)  + Zl!!  (x)  h"  . (y  ,z) 
t i 1 1 j J J 


where  Va(x) 
satisfy  the 


/\ 

and  Ia(x)  are  modal  voltages  an'  currents  which 
transmission  line  equations: 


(81) 


a; 


= -jK^(x)  Z^I?(x) 


(82) 


d rCt/  \ 
:s—  !•  (x) 
dx  1 ' ' 


■jK.  (x)  Y . V?(x) 
1 11 


The  single  mode  tranverse-to-z  magnetic  field  corres- 
ponding to  eigenvalue  of  the  LSM  and  E modal  subsets 

X 

must  be  equal.  Thus,  since  h' (x,y)  = h' (y,z)  s 0, 

(83)  I|(x)h^^  (y,z)y^  = I’  e"^'^^e^(x)e^^  (x,y)^ 


Using  (71)  in  (72)  , and  letting  A = ^l^ti^^i'^i 

results  in  an  expression  for  e'  (x,y)  in  terms  of  the  x 

xnm 

A 

dependent  modal  current  distribution  I^(x):* 


(84) 


e'  „(x,y) 
xnm' 


-3- 


N'  E (x) 
nm  r 


!• 

n 


(x) sin 


mlly 

b 


*in  equations  (84)  and  (86)  , the  double  subscrip b is 
used  to  explicitly  indicate  x and  y dependence. 


Note  that  for  m = 0,  the  LSM  mode  does  not  exist. 

Similarly,  the  single  mode  tranverse-to-z  electric  field 
corresponding  to  eigenvalue  Y of  the  LSE  and  H modal  sub- 
sets are  equal  giving 

— . j Y 2 

(85)  VJ(x)a"j(y,z)j(^  = -VJe  ‘ 

A 

Letting  B^=  and  using  (75)  in  (76)  results 

in 


(86) 


V" (x) cos 
n 


mny 

b 


The  coefficients  N^  and  appearing  in  equations  (84) 
and  (86)  are  normalization  constants  determined  by  appli- 
cation of  equation  (70)  . Complete  expressions  for  the 

normalizations  are  given  in  Appendix  C. 


Since  the  inhomogene ously  (in  x)  loaded  guide  is 
symmetric  about  the  midplane  (x  = 0),  the  modes  will  be 
either:  symmetric,  corresponding  to  an  open  circuit  plane 

at  X = 0 for  LSE  modes,  or  short  circuit  plane  for  LSM 
modes;  or  anti-  /mmetric,  for  the  converse.  Consider  the 
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equivalent  transmission  line  representation  of  wave  propa- 
gation in  X shown  in  Figure  27.  The  complete  current  and 
voltage  distributions  may  be  written  down  by  inspection  for 
either  of  the  indicated  terminations.  The  resultant 
distributions  are  then  appropriately  inserted  in  (84)  or 
(86)  to  obtain  the  cross-sectional  dependence  of  the 
X components  of  LSM  dlectric  and  LSE  magnetic  mode  functions. 
The  vector  mode  functions  are  summarized  in  Figures  28 
and  29  • 

In  normal  operation  as  a dual  frequency  phased  array 
element, the  excited  (propagating)  modes  of  the  symmetric 
twin  dielectric  slab  loaded  rectangular  waveguide  wxll  be 
the  LSEjiq  mode  in  tlie  low  frequency  band,  and  LSE^^q  and 
LSE2Q  modes  in  the  high  frequency  band.  At  either  band, 
the  tendency  will  be  for  the  field  strength  interior 
tc  the  dielectric  slabs  to  exceed  the  field  strength 
elsewhere.  This  characteristic  is  clearly  evident  in  the 
equations  of  Figure  29  for  slow  wave  propagation.  Exam- 
ination of  the  modal  voltage  expressions  shows  that  for 
<n  " proportional  to  cosh(|K^x|)  in 

"6<x<0,  and  to  sinh  Ik  I (x+a/2)  in  -a/2<x<-3-f.  Thus, 

' n ' 

for  greatly  slowed  waves,  i.e.,  Ypjj,  - the 

eyio(x,y)  is  nearly  exponential  in  the  air  regions,  and 
nearly  constant  in  the  dielectric.  Similar  characteristics 
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a.re  evident  for  the  LSE2Q  mode  (the  first  anti -symmetric  in 
X mode)  . 

Typical  e"^Q(x,y)  distributions  are  shown  in  Figures 
30  and  31  for  an  element  operating  at  2.5  GHz  and  6.0 
GHz,  respectively.  The  element  is  a WR187  guide  with  .250" 

(=6)  slabs  of  = 9 dielectric  located  .325"  (=a ) from 
either  narrow  wall.  At  low  frequency,  the  distribution  is 
roughly  uniform  between  slabs,  with  some  field  concentration 
in  the  vicinity  of  the  interior  air-dielectric  interfaces. 

In  the  region  |x|>6+6,  the  field  behaves  very  nearly  like 
cos{IIx/A).  At  6 GHz,  the  e'yj^Q(x,y)  distribution  is  entirely 
different,  showing  well  defined  field  concentration  about 
the  dielectric,  with  very  low  field  strength  in  the  air- 
filled  regions.  The  distributions  are  roughly  symmetric 
about  the  slabs,  with  non-zero  field  at  the  guide  center. 

As  the  dielectric  constant  is  increased  the  fields  become 
more  heavily  concentrated  in  the  dielectric,  and,  consequently, 
the  field  strength  at  the  guide  center  approaches  zero. 

At  6 GHz,  the  antisymmetric  LSE2Q  mode  is  also  propa- 
gating, and  all  other  modes  are  well  beyond  cut-off.  The 
®y20^^'^^  distribution  is  shown  overlayed  on  the  e^^Q(x,y) 
distribution  in  Figure  32.  This  comparison  shov/s  that  the 
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NORAAALIZED  FIELD 


NORMALIZED  FIELD 


Figure  31,  (x,y)  Distribution  for  WR187  Bifurcated 

Guide  With  .250"Thick  = 9 Loading.  6 GHz 
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LSE^q  and  LSE2Q  modal  field  distributions  are  quite  similar, 
differing  by  no  more  than  a few  percent  in  relative  magnitude 
for  |x/a|>.15/  but  having  opposite  symmetries.  It  is 
this  high  frequency  propagation  characteristic  and  the  fact 
that  /Yj^q-1-0  for  appropriately  chosen  dielectric 
constants  and  guide  geometries  which  provide  the  unique 
dual  frequency  array  element  potential  of  the  symmetric 
twin  dielectric  slab  loaded  rectangular  waveguide.  For, 
assuming  the  functions  are  exactly  identical  in  magnitude 
and  that  ‘r'2o/'^10  ” ^^^10 

magnitude  control  only,  may  be  excited  such  that  two  in- 
dependent phase  centers  are  located  at  roughly  the  positions 
of  the  slabs. 

The  similarity  of  the  e”j^Q(x,y)  and  dist- 

ributions, and  hence,  the  achieveable  high  frequency  phase 
center  independence,  is  directly  related  to  dielectric  con- 
stant for  fixed  geometry.  As  the  dielectric  constant  is 
decreased,  holding  cross-section  fixed,  the  LSE^^q  distri- 
bution approaches  the  TE^^q  distribution  of  empty  guide; 
and  the  LSE2q  approaches  the  empty  guide  TE2Q  distribution. 
These  trends  are  evident  in  Figure  33,  where  the  e^^Q(x,y) 
and  6^20  distributions  are  shown  overlayed  for  the 

WR187  guide  with  = 5 loading.  The  departure  in  magnitude 
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NORMALIZED  FIELD 


Figure  33.  ^yj^Q  ©^20  Distrubutions  for  WR187  Bifurcated 

Guide  With  .250"  Thick  = 5 Loading.  6 GHz 


between  the  distributions  is  significantly  greater  for  the 
lower  dielectric  constant  than  for  the  higher. 

While  the  above  discussion  shows  that  there  is  a strong 
influence  of  dielectric  constant  on  achieveable  high  freq- 
uency phase  center  independence,  it  should  not  be  construed 
that  large  dielectric  constant  is  generally  preferable  to 
low  dielectric  constant.  In  particular,  it  will  be  shown 
in  section  4-I  that  for  certain  geometries,  dielectric 
constants  on  the  order  of  = 9 may  lead  to  large  aperture 
reflections  which  are  difficult,  if  not  impossible,  to 
match  outJ 
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3 . 2 Mode  Spectrum 

The  modal  spectrum  of  the  symmetric  twin  dielectric 
slab  loaded  rectangular  waveguide  is  obtained  via  a 
transverse  resonance  procedure.  Representing  the  loaded 
guide  as  an  E or  H mode  transmission  line  in  x,  as  in 
Figure  27  » and  requiring  that  x = 0 be  either  an  open  or 
short  circuit  plane  results  in  four  dispersion  relations, 
in  X,  of  the  form. 


(87) 


where 


(881 


k"  = k"-  (S  - 


n 


'B 


nm 


(89) 


+k^ (e  -1) 
£n  n n r 


k = 2IT/X  is  the  free  space  wave  number,  mlT/B  is  the  y 

directed  wave  number,  and  y ™ is  the  z directed  wave- 

nm 

number.  One  dispersion  relation  is  obtained  for  each 
symmetry  condition,  in  x,  of  each  modal  subset.  The 
four  dispersion  relations  are  given  in  Table  2.  It 
shoiil.d  be  noted  that  the  forms  given  are  computationally 
unstable  due  to  both  the  various  tangent  evaluations. 
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and  the  indicated  divisions.  ''In^practice , the  expressions 
are  converted  to  combinations  of  sines  and.,  cosines , and  all 
divisions  are  removed.  ■ 


In  the  LSE  and  LSM  modal  subsets,  the  mode  indexing 
is  (n,m) . The  first  index,  n,  is  associated  with  the 
infinite  sequence  of  zeros  of  D(k,k^)  which  are  arrayed 
along  the  real  axis,  as  illustrated  in  Figure  34. 

The  first  zero  is  always  associated  with  the  even  symmetry 
solution  of  the  modal  subset  (i.e.  open  circuit  symmetry 
for  LSE  modes;  short  circuit  symmetry  for  LSM  modes) . 
Thereafter,  the  roots  of  the  even  and  odd  symmetry  disper- 
sion relations  are  interleaved.  Hence, for  even  symmetry, 
n is  always  odd,  and  conversely.  The  second  index,  m, 
appears  explicitly  in  the  auxilliary  dispersion  equation, 
(88),  and  is  directly  interpreted  as  the  order  of  the  y 
variation  of  the  modal  field. 

The  potential  of  the  symmetric  twin  dielectric  slab 
loaded  rectangular  waveguide  as  a dual  frequency  array 
element  arises,  in  part,  due  to  the  unique  migration  of  the 
roots  Yio  "<'20  frequency.  Figure  35  shows  an 

LSE  mode  dispersion  diagram  for  half  height  WR187  guide 
with  e = 9 loading.  Only  the  (n,o),  n = 1,2,3  and  (n,l), 
n = 1,2  roots  are  shown,  other  LSE  roots  being  considerably 
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Figure  35.  LSE  Dispersion  Diagram  for  Bifurcated 
WR187  Guide  With  250"  Thick  = 9 Loading 


more  evanescent.  ?or  kA/2  <1.40,  only  the  LSE^q  mode  is 
propagating.  As  frequency  is  increased,  the  1‘SE2q  mode 
begins  to  propagate,  and  Y20  ^^Pidly  approaches  Y]^q  ' 
ratio  Y2o/‘YiO  nearly  unity  for  ' I<A/2  > 2.5.  As 

frequency  is  further  increased,  the  LSEj^j^  and  LSE22^ 
modes  begin  to  propagate,  though  they  may,  of  course,  be 
pushed  further  out  by  decreasing  the  guide  height.  The 
result  is  that  over  the  range  2.5  < kA/2  < 3.0,  the 
propagating  modes  of  the  structure  have  virtually  identical 
dispersion.  Hence,  if  in  this  range,  the  LSE^^q  and  LSE2Q 
mode  functions  differs  only  in  symmetry,  as  in  Figure  32, 
the  half  height  WR187  guide  with  e^=9  slab  loading  will 
support  two  independent  phase  centers  over  an  18%  frequency 
band  centered  at  kA/2  = 2.75;  and  over  multiple  guide 
wavelengths,  2U/y^Q,  in  z. 

The  dashed  line,  y “ Figure  35  may  be  used 

to  estimate  the  mismatch  at  the  feedguide-free  space 
interface  for  broadside  excitation  (only  the  LSE^^q  is 
excited) . The  modal  admittance  of  the  LSE^^q  mode  is 
Yj^q/cjm,  and  the  modal  admittance  of  the  dominant  space 
harmonic  is  k/wp  at  broadside.  As  an  estimate,  then, 
the  reflection  coefficient  at  the  aperture  will  be  in  the 
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neighborhood  of 
(90) 

for  broadside  excitation.  For  the  WR187  guide  with  e^=9 

loading,  operating  at  KA/2  = 2.99,  equation  (4-39)  gives 
|r|  = .383,  whereas  the  exact  value  for  a rectangular  thin 
wall  array  of  these  guides  is  jr]  = .424.  For  the  same 
guide  operating  at  kA/2  = 1.25,  equation  (90)  gives 
|r[  = .110,  and  the  exact  value  is  |r|  = .285.  In  both 
cases,  the  implication  is  that  the  aperture  susceptance, 
which  is  ignored  in  (90)  , is  significant.  This  is  not 

entirely  surprising,  since  it  is  well  known*  that  equation 
(90)  is  exact  for  thin  walled  rectangular  grid  arrays  of 
empty  rectangular  grid  guide,  for  which  the  set  of  trans- 
verse wavemimbers  of  the  feedguide  is  the  set  of  transverse 
wavenumbers  for  the  unit  cell  guide. 

For  fixed  guide  wall  dimensions,  the  parameters  which 
most  strongly  influence  the  dispersion  curves  are  dielectric 
constant  and  slab  thickness.  The  location  of  the  slabs, 
denoted  by  the  ratio  a/3,  has  second  order  effects  in  the 
regions  .5<  a/3  £ 2.0.  a/3  ratios  outside  this  region  are 

not  of  interest  due  to  the  irregularly  spaced  high  freq- 

*See,  for  example,  Amitay,  Galindo,  and  Wu,(10)  pg  132ff. 
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uency  phase  centers  which  would  result. 

Figures  36  thorough  39  ^re  LSE  dispersion  diagrams 
for  half  height  WR137  guide  with  6 = .250"  and  a/3  = 1.00. 
The  curves  progress  from  loading  of  e^=3  to  £^.=9*  As  the 
loading  increases,  Y20  approaches  in  general,  and  the 

trend  is  toward  steeper  slopes. 

The  migration  of  the  LSE^q  cutoff  frequency  toward 
lower  frequency  produces  one  of  the  critical  trade-offs 
inherent  to  the  design  of  the  symmetric  twin  dielectric 
slab  rectangular  waveguide  dual  frequency  array  element. 

In  general,  the  overriding  array  design  criterion  will  be 
to  minimize  the  number  of  radiators  (or,  equivalently 
maximize  array  cell  size).  Thus,  for  dual  frequency 
operation,  the  high  frequency  band  center  operating  point 
will  be  in  the  vicinity  of  kA/2  = n.  In  this  region, 
the  higher  dielectric  constants  (e^  = 7,9)  provide  V20'^^10 
ratios  very  close  to  unity,  but  the  LSE^q  mode  is  propa- 
gating. Since  it  is  necessary  to  push  this  mode  out,  lower 
dielectric  constant  is  required. 

Similar  results  are  obtained  by  holding  fixed  and 
varying  only  slab  thickness,  6.  Figures  40  through  43 
are  LSE  dispersion  diagrams  for  four  slab  thicknesses  in 
half  height  WR137  guide.  The  dielectric  constant  is  5, 
and  a/3  = 1.00.  In  general,  the  behavior  with  thickness 
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Figure  39.  LSE  Dispersion  Diagram  for  Bifurcated 
WR137  Guide  With  .250"  Thick  e = 9 Loading 
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Figure  40.  LSE  Dispersion  Diagram  for  Bifurcated 
WR137  Guide  With  .063"  Thick  = 5 Loading 


93 


Figure  42.  LSE  Dispersion  Diagram  for  Bifurcated 
WR137  Guide  With  .188"  Thick  e =5  Loading 
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Figure  43.  LSE  Dispersion  Diagram  for  Bifurcated 
WR137  Guide  With  .375"  Thick  = 5 Loading 
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is  very  like  the  behavior  with  e^.  However,  comparison 
of  Figures  39  and  43  shows  that  high  dielectric 
constant  is  preferable  to  extreme  thickness.  In  the  two 
figures,  the  LSE^q  mode  enters  at  roughly  the  same  freq- 
uency, but  the  ratio  of  z directed  wave  numbers,  y20^\o' 
is  significantly  nearer  to  unity  for  e^=9  than  for 


‘1.0  ARRAY  APERTURE  DESIGN 


In  this  section,  the  trade-offs  leading  to  a 
practical  array  apertxjre  design  are  presented  for  the 
hypothetical  array  performance  given  in  Table  3.  Hie 
operating  bands  are  15%  centered  at  4 and  8 GHz.  Sixty 
degree  (60®)  principle  pleuie  scam  coverage  is  required. 
The  array  is  to  provide  30  db  gain  for  60®  principle 
plane  scan  at  4 GHz.  The  first  sidelobe  level  is  to 
be  below  20  db  auid  RMS  sidelobe  levels  are  to  be  below 
35  db.  Nominal  feed  losses  of  1.5  db  at  4 GHz  and  2.8 
db  at  8 GHz  are  assumed.  To  reduce  the  inherent 
difficulties  in  matching  out  the  aperture,  an  aperture 
mismatch  loss  at  broadside  of  no  greater  than  .6  db 
is  required. 

To  provide  the  required  performance,  a Taylor, 
n=3,  25  db  SLL  distribution  is  selected,  resulting  in 
an  aperture  gain  of  37  db  at  4 GHz. 

The  principle  result  of  this  section  is  the 
determination  of  an  element/grid  configuration  which 
minimizes  element  count  while  holding  high  frequency 
grating  lobe  contributions  to  levels  consistent  with 
the  sidelobe  requirements. 
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4.8  GHz  +7.5% 

+60°  in  either  principle 
plane 

30  db  0 4 GHz 

20  db 
35  db 

1.5  db  0 4 GHz 

2.8  db  0 8 GHz 

<.6  db  in  both  bands 


Frequency 
Scan  Coverage 

Apert\ire  Gain  0 60° 
Principle  Plane  Scan 

1st  Side lobe  Level 

RMS  Side  lobe  Level 

Feed  Losses 

Aperture  Mismatch 
Loss  0 0°  Scan 


Table  3 

Prescribed  Array  Performance 


4. 1 Aperture  Design  Trade-Offs 


It  is  well-known that  an  equilateral  triangle 
lattice  configuration  minimizes  element  count  for  a 
given  scan  requirement.  As  a basis  for  comparison, 
it  is  therefore,  convenient  to  first  consider  a lattice 
specifically  taylored  to  minimize  element  count  at 
low  frequency.  This  lattice  and  the  near-in  grating 
lobe  diagrams  for  4 and  8 GHz  are  shown  in  Figure  44. 
The  lattice  base  is  1.780"  (.602X  @ 4 GHz). 

It  is  immediately  apparent  from  Figure  44c  that 
the  six  near-in  high  frequency  grating  lobes  will 
migrate  well  inside  the  coverage  sector  for  all  scan 
directions  and  pose  a potential  limitation  on  achieve- 
able  peak  sidelobe  level  and  main  beam  gain.  This 
difficulty  may  be  alleiviated  somewhat  by  appropriate 
selection  of  the  element  configuration  and  high  freq- 
uency excitation  modifier,  R,  for  the  LSEzo  mode  as 
discussed  in  section  2.2.  However,  it  should  be  noted 
that  the  six  grating  lobe  locations  actually  represent 
twelve  independent  beams  (six  each  E and  H with  respect 
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Figure  44.  Equilateral  Triangle  Lattice  Configuration 
Which  Minimizes  Element  Count  @ 4 GHz 


to  the  array  normal) . Since  there  are  only  nine  design 
parameters  (half  aperture  dimensions  A and  B;  septum 
thickness;  4;  relative  permittivity,  e;  slab  thickness, 
6;  ratio  a/3;  lattice  spacings,  d and  d ; and  R)  , it 

X y 

is  clear  that  there  is  insufficient  control  for  the 
cancellation  of  all  beams. 

The  only  available  means  of  cancelling  the  beams 
residing  aJ.ong  the  v axis,  the  (0,1)  and  (0,-1)  lobes 
is  to  separate  the  half  apertures  by  one  half  the  y 
lattice  spacing  d^,  resulting  in  a scanning  sub-array 
pattern  with  nulls  coincident  with  the  grating  lobe. 
This  coincidence  is  maintained  for  H-plane  scanning, 
however  degrades  in  all  other  scan  planes  due  to  the 
imbalance  in  half -aperture  reflections  induced  by  the 
phase  taper.  In  the  circumstance  that  the  half  aper- 
tures are  spaced  at  less  than  0.5  d^,  a significant 
fraction  of  the  radiated  power  is  delivered  to  the 
E-mode  (0,1)  and  (0,-1)  beams  for  all  scan  directions, 
including  broadside.  In  calculations  with  the  half 
apertures  separated  by  0.26  d^,  it  was  found  that  as 
much  as  20%  of  the  power  was  delivered  to  each  of  the 


two  E-mode  beams  at  broadside  scan. 


Cancellation  of  the  remaining  eight  beams  is  a 
considerably  more  difficult  problem,  and  turns  out  to 
be  virtually  impossible  for  this  large  lattice  spacing. 
Again,  effective  reduction  of  radiated  power  into  the 
unwanted  lobes  must  be  obtained  by  locating  the  null 
of  a scanning  subarray  pattern  at  or  near  the  grating 
lobe  location.  One  subarray  spacing  has  already  been 
used  to  cancel  the  oi  axis  lobes;  the  remaining  degree 
of  freedom  is  therefore,  the  x separation  of  the  appar- 
ent high  frequency  phase  centers.  In  the  event  that 
these  separations  can  be  extended  to  0.5  d^^,  good  grating 
lobe  rejection  can  be  achieved.  However,  as  a practical 
matter, such  wide  x displacements  of  apparent  high  freq- 
uency phase  centers  are  not  possible. 

As  was  shown  in  section  3.2,  a reasonable  upper 
limit  on  aperture  x dimension  is  in  the  vicinity  of 
A=A  at  the  center  of  the  high  frequency  band  for  moder- 
ate relative  permittivity  loading  (£^.=5)  . Larger 
aperture  dimensions  can  be  achieved  by  reducing  the 
permittivity,  or  using  very  thin  slabs;  however,  such 
an  approach  is  generally  counter-productive  in  that 
bandwidth  is  reduced  and  H-plane  scan  loss  is  increased 
at  high  frequency  due  to  the  widening  dissimilarity 


between  LSEio  and  LSE20  nodes  which  occurs  for  parameter 
variations  in  these  directions.  Consequently,  in  order 
to  achieve  the  wide  apparent  phase  center  displacement 
for  this  grid,  it  is  necessary  to  place  the  slabs 
very  close  to  the  narrow  feedguide  walls,  giving  a/3 
ratios  of,  typically,  0.56  to  0.62.  This  is  not  a 
very  attractive  solution  since  it  produces  considerable 
aperture  field  asymmetry  about  the  apparent  phase 
centers,  thereby  invalidating  the  primary  assumption 
that  the  high  frequency  aperttare  distributions  are 
roughly  symmetric  about  these  phase  centers.  As  a 
practical  matter,  then,  it  is  impossible  to  kill  off 
the  eight  remaining  beams. 

The  best  result  obtained  for  the  1.78"  base 
equilateral  grid  are  summarized  in  Table  4.  The 
element  dimensions  are  given  in  the  table  caption. 
Clearly,  the  1.05  db  power  loss  to  the  grating  lobes 
is  intolerable,  and  it  is  necessary  to  reduce  the  grid 
size  significantly. 

As  mentioned  above,  cancellation  of  off-axis  high 
frequency  grating  lobes  requires  that  the  apparent 
phase  centers  be  separated  by  roughly  0,5  d . In 
addition,  to  maintain  reasonable  aperture  field  symmetry 
about  these  phase  centers,  ratios  a/3  should  lie  in 
the  approximate  range  0.8  < a/3  < 1.25.  Allowing  for 
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4 GHz 


8 GHz 


Broadside  Reflection  Loss  (db)  1.18 
*H-plane  Scan  Loss  at  60°  (db)  6.44 


Power  Loss  to  Grating  Lobes 
(broadside  scan) 


Broadside  Gain  Loss  (db)  1.18 


1.05 

6.80 


1.05 


2.43 


Table  4 

Performance  of  Dual  Frequency  Element  in  1.78"  Base 
Equilateral  Traingular  Grid  at  4 and  8 GHz.  A = 1. 
S = .412",  4 = .358",  d = 1.7  80",  d = 1.540",  a 

3 = .301",  6 = .198",  = 5.0  ^ 


*Includes  cos  6 beam  broadening  factor. 


a wall  thickness  of  0.062"/  a practical  set  of  element 

and  grid  x-dimensions  is  then  A = 1.376",  d = 1.500", 

X 

and  3+6/2  = .375",  resulting  in  the  specification  of 
slab  thickness,  6,  as  no  greater  than  .130",  or 
6/A<.094  and  a/3<.835.  Prom  Figures  36  through  43, 
it  is  seen  that  a relative  permittivity  of  5 will 
provide  the  desired  feedguide  characteristics  in  both 
frequency  bands  for  this  aperture  size  with  a/3  = 1 
euid  6/A  = .091.  Bringing  the  a/3  ratio  into  the 
required  range  produces  only  a small  perturbation 
on  the  dispersion  curves  shown  in  Figure  41.  Consequent- 
ly, the  selection  of  a=.247",  3=*309",  and  6=.132" 
will  provide  the  necessary  control  of  off-axis  grating 
lobes  while  keeping  higher  order  feedguide  modes  well 
below  cutoff.  By  requiring  that  the  LSEn  mode  be 
attenuated  by  8 db  per  wavelength  at  the  high  end  of 
the  8 GHz  band,  the  y aperture  dimensioned  is  obtained 
as  B = . 400" . 

The  aperture/grid  parameters  determined  so  far 
guarantee  some  cancellation  of  the  off-axis  grating 
lobes,  and  it  might  be  assumed  that  by  reducing  d^ 
for  the  original  grid  by  the  ratio  1.5/1.78  would 
result  in  an  acceptable  configuration,  giving  a 29% 
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VBmm 


increase  in  element  count  over  the  "optimum"  number. 
Unfortunately  this  spacing  is  still  on  the  large  size 
and  results  in  considerable  power  loss  to  the  grating 
lobes  at  relatively  small  scan  angles.  It  is  there- 
fore, necessary  to  further  reduce  the  cell  along  the 
E-plane  to  d^  = .960". 

The  final  configuration,  shown  in  Figure  45, 
with  its  low  and  high  frequency  grating  lobe  diagrams 
has  a cell  area  roughly  half  that  for  the  grid  optimized 
for  element  count  at  low  frequency,  and  results  in  an 
array  of  2444  elements.  With  five  phase  shifters  per 
element,  this  increase  seems  (at  first  glance)  rather 
unattractive.  However,  this  comparison  is  quite  mis- 
leading. The  twin  dielectric  slab  dual  frequency  array 
element  concept  provides  for  simultaneous  excitation 
by  two  entirely  independent  feed  systems,  and  con- 
sequently the  simultaneous  radiation  of  two  indepen- 
dent beams  at  two  widely  separated  frequency  bands. 

If  the  alternatives  for  multifrequency  operation  are 
considered,  the  dual  frequency  element  is  suddenly 
very  attractive. 

One  such  alternative  is  the  use  of  wideband 
elements  in  the  grid  depicted  in  Figure  44  to  obtain 
simultaneous  aperture  usage.  The  obvious  disadvantage 
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(o)  Eltm«nt/6rid  Configuration 


(b)  Low  Froquoncy  Grating  Lobo  (c) High  Fraqutncy  Grating  Lob* 
Diagram  1.96  \oy*3.07  Diogram-^^Dx  *.98, ^/0y=  1.54 


Figure  45.  Triangular  Grid  Configuration  for  Dual  Frequency 
Operation  over  16%  Bands  Centered  at  4 GHz  and  8 GHz 


of  this  scheme  is  that  it  requires  eight  phase  shifters 
and  four  dip  lexers  per  element  quartet  or  unit  cell. 

In  comparison  with  the  dual  frequency  element/  this 
configuration  requires  twice  the  number  of  low  freq- 
uency controls,  and  half  the  number  of  high  freq- 
uency controls  per  unit  area  of  the  array,  plus  di- 
plexers . 

A second  alternative,  in  a broad  sense,  is  to 
design  two  entirely  independent  single  frequency 
systems.  However,  it  is  clear  that  only  under  very 
special  circumstances  could  this  system  be  considered 
a viable  multi frequency  concept. 

Predicted  principle  scan  plane  performance  for 
the  configuration  in  Figure  45,  is  shown  in  Figures  46 
through  54,  at  the  end  and  midpoints  of  the  4 and  8 
GHz  bands.  In  these  figures,  the  performance  measure 
is  taken  as  the  power  transmission  coefficient  int"> 
the  radiating  beam(s)  . Mainbeam  scan  loss  is  deter- 
mined by  adding  10  logic  (cos0)  to  the  curve  for  the 
(0,0)  beam. 

Figures  46,  47,  and  48  show  E and  H plane  per- 
formance at  the  low  end,  middle,  and  high  end  of  the 
4 GHz  band,  respectively,  for  scan  out  to  sin 0=. 975. 
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Figure  49.  Propagating  Beam  Power  Levels  - H-Plane  Scan, 

f = 7.36  GHz 
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Figure  51.  Propagating  Beam  Power  Levels  - H-Plane  Scan, 

f = 8.64  GHz 
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Figure  52.  Propagating  Bean  Power  Levels  - E-Plane  Scan 

f = 7.36  GHz 
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Figure  53,  Propagating  Beam  Power  Levels  - E-Plane  Scan, 

f = 8.0  GHz 
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POWER  TRANSMISSION  FACTOR  (dB'^ 


In  this  range,  no  grating  lobes  enter  real  space.  Maxi- 
mum broadside  mismatch  loss  in  this  band  is  less  than 
.3  db  and  is  readily  matched  out  by  any  number  of  means. 
The  E-plane  fall-off  out  to  approximately  60°  is  some- 
what better  than  might  be  expected  for  the  subarray 
array  factor  at  these  spacings,  and  decreases  at  the 
upper  end  of  the  band.  The  H-plane  fall -off  is  typical 
of  a planar  array  of  rectangular  apertures. 

Figures  49,  50  and  51  show  the  H-plane  power  levels 
in  the  propagating  beams  in  the  8 GHz  band.  To  eval- 
uate the  array  characteristics  over  this  band,  the  LSEio 
and  LSEzo  feedguide  modes  are  assumed  to  be  generated 
2A  10  behind  the  aperture  plane,  where  X i o is  the 
LSEio  guide  wavelength  evaluated  at  8 GHz.  This  long 
feedguide  phase  length  results  in  considerable  excita- 
tion of  the  (-1,-1)  E-mode  at  the  band  extremes.  Evi- 
dently, this  path  must  be  shortened  for  such  wide  band 
operation . 

At  midband.  Figure  50,  the  (-1,-1)  and  (-1,0) 
grating  lobes  are  well  within  the  desired  range  of  rms 
sidelobe  level  and  will  not  result  in  any  significant 
perturbation  of  the  far-out  sidelobe  region.  This 


level  of  cancellation  is  achieved  using  the  modifier 
R = 1.2  exp  { “jl4°  )/  which  seems  to  also  result  in 
significant  reduction  of  the  H-mode  (-1,-1)  and 
(-1,0)  beams  at  the  band  edges  as  seen  in  Figures  49 
and  51.  However,  the  modifier  clearly  has  little 
effect  on  the  E-mode  beams,  and  tliese  levels  must  be 
controlled  by  a proper  choice  of  a feedguide  phase 
length. 

In  Figure  51,  it  is  seen  that  the  (-2,-1)  H- 
mode  beam  is  heavily  excited  at  the  upper  end  of  the 
band.  However,  this  occurs  only  at  the  H-plane 
extreme  of  the  scan  volume,  and  will  be  of  only  minor 
consequence . 

The  8 GHz  band  E -plane  performance  is  shown  in 
Figures  52  through  54.  In  general,  the  best  perform- 
ance occurs  at  the  lower  end  of  the  band.  Due  to  the 
choice  of  y lattice  spacing,  the  grating  lobes  remain 
outside  real  space  throughout  most  of  the  scan  volume 
at  this  end  of  the  band,  and  are  only  moderately 
excited  upon  entering.  From  the  results  at  8 and  8.64 
GHz,  it  is  clear  that  this  is  the  most  effective  means 
of  controlling  spurious  beam  levels  in  this  plane. 

A;  the  high  end  of  the  band,  the  (-1,-1)  E 
and  H moae  beams  are  very  heavily  excited  and  will 


have  considerable  impact  on  the  general  side  lobe  level 
for  scanning  beyond  25°.  However,  except  for  the  snail 
dips  occuring  near  (-1,-1)  lobe  incipience,  there  is 
little  effect  on  main  beam  gain.  Again,  as  for  the 
H-plane  results,  a general  improvement  in  grating  lobe 
levels  may  be  expected  for  shorter  feedguide  phase 
lengths . 


4.2  Experimental  Evaluation  of  the  Bifurcated  Twin 


Dielectric  Slab  Loaded  Rectangular  Waveguide 


Dual  Frequency  Element 


The  dual  frequency  element  design  shown  in 
Figure  45  was  built  and  tested  in  H-plane  waveguide 
simulators  over  an  8%  frequency  band  centered  at  4.6 
GHz  and  a 16%  band  centered  at  8 GHz.  In  general, 
the  experimental  results  were  in  excellent  agreement 
with  predictions. 

The  simulators  are  shown  in  Figure  55.  They 
are  constructed  of  brass  and  are  soft  soldered.  A 


single  element  section  and  parallel  wall  simulator 
section  is  used  for  both  bands,  resulting  in  near  broad- 
side simulation  in  the  upper  band,  and  wide  angle  scan 
simulation  in  the  lower  band.  The  parallel  wall 
section  has  2.250"  x .960"  cross-section  and  provides 
imaging  as  shown  in  Figure  56 . In  the  low  frequency 
band,  only  the  TEio  mode  propagates  in  the  simulator. 

In  the  high  frequency  band,  up  to  eight  modes  can  propagate 
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Figure  56.  Simulator  Imaging  of  Element  Section 
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The  element  section,  shown  in  Figure  56  with  the 

simulator  imaging,  consists  of  two  half  elements  and 

two  quarter  elements,  resulting  in  19. l'^  and  41° 

H-jJlane  scan  simulations  at  the  centers  of  the  high 

and  low  bands,  respectively.  The  dielectric  slabs 

are  stycast  HiK,  = 5.  The  element  section  is  7.2" 

long  (2A  at  4 GHz),  with  tapered  3.6"  300fi  card 

y 

loads  inserted  at  the  rear  of  the  section  on  either 
side  of  the  dielectric  slabs  and  the  feedguide  mid- 
plane. The  long  tapers  are  necessary  to  eliminate 
rearward  radiation  and  reduce  reflections  at  the  load 
discontinuities. 

Figure  57,  shows  measured  results  in  the  upper 
half  of  the  low  frequency  band.  Results  in  the  lower 
half  of  the  band  were  not  obtained  due  to  the  freq- 
uency limitation  of  the  network  analyzer.  To  ensure 
measurement  accuracy,  the  experimental  band  was 
sampled  discretely  in  80  MHz  increments,  and  a short 
circuit  reference  was  established  at  each  frequency 
step.  Measured  reflection  coefficient  magnitude  is 
in  the  ranee  .34  to  .41*  with  the  peak  and  minimum 
at  4.24  and  4.08  GHz,  respectively,  and  is  typical  of 
the  variation  in  H-plane  gain  loss  with  scan  observed 
in  many  broadside  matched  phased  arrays.  The  phase 
of  the  reflection  coefficient  is  nearly  constant. 

*For  H-plane  scan,  it  is  permissible  to  represent  feed- 
guide  port  pairs  (upper  and  lower  element  halves)  by 
a single  "effective  port".  Consequently,  the  reflec- 
tion coef  fici-''nt  magnitude  at  the  simulator  port  is  a 
factor  o’  /2  greater  than  that  in  either  half  element. 
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Measured  and  predicted  simulator  results  are 
shown  in  Figure  58.  The  calculated  resulcs  lie  well 
within  the  range  of  the  measurements.  Also  shown 
in  the  figure  is  a least  mean  square  straight  line 
fit  to  the  experimental  data  and  /(l-ZcosB)  , where 
9 is  the  simulator  angle.  These  two  curves,  in  com- 
parison with  the  predictions,  show  quite  clearly  that 
the  analytical  model  provides  an  excellent  descrip- 
tion of  the  array.  The  scatter  of  experimental 
data  about  the  theoretical  results  is  due  primarily 
to  errors  in  element  section  fabrication  which 
resulted  in  small  air  gaps  between  the  feedguide  broad- 
walls  and  the  dielectric  primarily  in  the  interior 
of  the  section.  By  assuming  a maximum  reflection 
coefficient  magnitude  of  0.035  at  the  internal  discon- 
tinuities, the  discrepancies  in  the  results  are  accounted 
for  throughout  the  band. 

In  the  8 GHz  band,  up  to  eight  waveguide  modes 
will  propagate  in  the  simulator.  Below  8.08  GHz, 
the  first  five  simulator  modes  will  propagate.  However, 
for  a properly  fabricated  element  section,  only  the  TEjo 
mode  is  excited.  Above  8.08  GHz,  three  of  the  eight 
modes  will  be  excited  by  the  array  interface.  These 
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FREQUENCY  (GHz) 


Figure  58.  Comparison  of  Measured  and  Predicted  Reflections 
Coefficient  Magnitude  at  the  Simulator  Port/  4.0  - 4.32  GHz 


are  the  TEio/  TE21,  and  TM21  modes,  corresponding  to 
the  (0,0);  H mode  (1,0)  and  (-1,-1);  and  E mode  (1,0) 
and  (-1,-1)  beams,  respectively. 

The  equivalent  network  representing  the  simulator 
discontinuity  is  a 3-port  below  8.08  GHz,  and  a 5- 
port  above  8.08  GHz.  Consequently,  below  the  (2,1) 
waveguide  mode  cut-off,  the  measured  reflection  coe- 
fficient at  the  simulator  port  is  given  as 

(91)  |Sii|^  = -1  +IS22I''  + |S33l"+  |S23|"y2/Y3+1S32|^  Y3/Y2 

where  the  ports  are  defined  in  Figure  59.  Above  the 
cut-off  frequency,  the  simulator  dominant  mode  self 
reflection  term  is  complicated  function  of  the  self 
and  cross  coupling  scattering  parameters  of  the  re- 
maining ports  in  the  network.  Since,  in  the  analysis 
presented  here,  the  interface  scattering  blocks 
Si  2 and  S22  are  unnecessary  for  the  determination  of 
element  performance,  they  are  not  calculated*,  an.d  it 


*The  calculation  o'  Si  2 and  S 2 2 requires  prohibitively 
large  amounts  of  computer  core.  Roughly  120K,  decimal, 
words  are  required  for  S 2 2 for  the  convergence  radii 
considered  here.  ~ 
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and  it  is  therefore,  not  possible  to  compare  theoretical 
and  measured  results  above  8.08  GHz. 

Measured  results  in  the  7.32  to  8.64  GHz  band  are 
shown  in  Figure  60.  As  for  the  4 GHz  band,  the  experi- 
mental band  was  sampled  discretely  to  ensure  measure- 
ment accuracy.  The  sampling  rate  is  roughly  every  160 
MHz,  with  an  exact  short  circuit  reference  established 
for  each  sample  point.  The  reflection  coefficient 
magnitude  is  in  the  range  .35  to  .46,  and  the  phase  is 
nearly  constant. 

In  the  measurements,  no  attempt  was  made  to 
load  terminate  the  higher  order  modes  of  the  simula- 
tor. This  leads  to  an  inherent  error  in  the  results 
which  is  associated  with  the  reactive  termination  of 
the  higher  order  modes  by  the  simulator  flare  transi- 
tion. This  error  should  be  insignificant  for  a reason- 
ably well  fabricated  element  section  since  the  higher 
order  beams  are  only  weakly  excited  above  8.08  GHz. 
However,  as  discussed  above,  some  irregularities  in 
element  section  fabrication  did  occur,  resulting  in 
week  excitation  of  the  TE20  simulator  mode. 
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Figure  60.  Measured  Simulator  Port  Impedance  7.32  - 8.64  GHz 

Sampled  in  160  MHz  Increments 


A comparison  of  measured  and  predicted  results  in  the 
7.32  to  S.OO  GHz  region  of  the  experimental  band  is  shown 
in  Figure  61.  Also  shown  in  the  figure  is  a straight 
line  least  square  fit  to  the  experimental  data.  As  in  the 
4 Ghz  band,  the  theoretical  data  falls  within  the  range  of 
the  experimental  results,  and  the  least  square  fit  has 
approximately  the  same  slope  and  magnitude  as  the  calcu- 
lated curve.  The  maximum  deviation  of  measured  reflec- 
tion coefficient  from  the  theoretical  value  is  .051  and 


occurs  at  8 GHz. 
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Figure  61.  Comparison  of  Measured  and  Predicted  Reflection 
Coefficient  Magnitude  at  the  Simulator  Port,  7.32  - 8.00  GHz 


5.0  ELEMENT  EXCITER  DESIGN 


Tne  design  of  an  exciter  for  the  slab  loaded  dual  fre- 
quency element  tjresents  a particularly  intriguing  problem. 

As  demonstrated  in  section  3,  the  dominant  mode  dispersion 
in  the  inhomogeneous ly  loaded  guide  is  roughly  linear  with 
frequency  for  practical  element  configurations.  However,  the 
slope  of  ^ /k  is,  in  general,  considerably  greater  than  unity. 
Consequently,  the  use  of  a bidirectional  exciter,  such  as  a 
stub  or  slot,  requires  load  termination  at  the  back  of  the 
feedguide  to  ensure  proper  aperture  excitation,  and  results 
in  a 3dB  power  loss. 

5 . 1 Exciter  Concept 

A unique  uni-directional  exciter  concept,  shown  in 
Figure  62,  has  been  developed  which  alleviates  this  dif f ir*”’' ty . Vhe 
exciter  consists  of  three  stripline  fed  flared  notch  antennas 

configured  to  provide  maxiiom  coupling  to  the  driven 
feedguide  modes  in  either  band.  The  center  probe  is  the  low 
frequency  exciter,  and  is  placed  well  in  advance  of  the  high 
frequency  exciters  (outer  two  probes)  to  maximize  the  low 
frequency  isolation.  The  high  frequency  exciters  butt  directly 
into  the  dielectric  slabs. 
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TOWARD  APERTURE  PLANE 


Figure  62.  Stripline  Fed  Notch  Exciter  for  Twin  Dielectric  Slab 
Loaded  Rectangular  Waveguide  Dual  Frequency  Array  Element 


The  basic  exciter  is  shown  in  Figure  6 3.  It  is  formed 
by  symmetrically  etching  the  outer  conductor  of  symmetric 
stripline  board  to  form  flared  notches  which  terminate  with 
maximum  aperture  at  the  board  edge  and  short  circuit  at  the 
notch  bottom.  The  stripline  center  conductor  is  configured 
to  cross  the  notch  regic-n  at  a right  angle  to  the  notch  center 
line,  and  terminates  in  an  open  circuit.  By  appropriately 
selecting  the  distance  X2»  from  the  center  conductor  center 
line  to  the  notch  short  circuit,  and  the  distance  y2>  ftom 
the  notch  edge  to  the  center  conductor  open  circuit,  the 
exciter  is  matched  in  the  feedguide  environment  over  the 
operating  band. 

A particularly  attractive  feature  of  the  flared  notch 
exciter  's  that  the  stripline  board  is  plugged  into  the  back 
of  the  element.  This  considerably  simplifies  feed  design. 

In  addition,  the  phase  shifter  and  exciter  may  be  integrated 
into  a single  unit. 

The  exciter  geometry  shown  in  Figure  62  results  in 
natural  isolation  between  the  low  and  high  frequency  probes 
in  the  low  frequency  band.  Since  the  stripline  outer 
conductor  is  etched  only  near  the  board  edge,  the  center  probe 
provides  a short  circuit  bifurcation  of  the  guide  in  the 


vicinity  of  the  high  frequency  probes.  Consequently,  the 
outer  probes  are  effectively  below  cut-off  guide  in  the  low 
frequency  band,  and  their  excitation  is  governed  by  the 
longitudinal  separation  between  center  probe  notch  bottom 
and  leading  edges  of  the  outer  probes. 

In  the  high  frequency  band,  the  outer  probes  couple  into 
the  LSEj^q  mode  of  the  loaded  half  width  guide  and  some  natural 
probe  isolation  is  achieved  due  to  the  low  field  strength 
along  the  outer  conductors  of  the  center  probe.  However, 
the  discontinuity  at  the  center  probe  termination  results 
in  scattering  back  into  the  low  frequency  port.  This  diffi- 
culty may  be  removed  by  introducing  an  appropriately  placed 
shorting  stub  along  the  low  frequency  probe  center  conductor, 
but  the  impact  of  this  technique  on  high  frequency  aperture 
field  distributions  has  not  been  determined  and  remains  a 
design  problem  for  future  consideration. 

5 . 2 Experimental  Investigation  of  the  Stripline  Fed  Flared 
Notch  Exciter 

Experimental  investigation  of  the  stripline  fed  flared 
notch  exciter  and  exciter  design  were  initiated  during  the 
contract  period.  In  general,  the  anticipated  results  were 
obtained,  and  demonstrate  the  viability  of  the  concept. 
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Notch  exciter  designs  were  developed  to  provide  better  than 
2:1  VSWR  looking  into  the  load  terminated  feedguide  over 
greater  than  10%  band  widths  / and  greater  than  50dB  probe 
isolation  in  the  low  frequency  band. 

The  baseline  exciter  design,  a scaled  version  of  a 
previously  designed  antenna  element,  is  shown  in  Figure  64. 
The  dielectric  is  Rexolite.  The  etched  notch  region  has 
uniform  width  of  .058”  for  a length  of  .154"  from  the  notch 
short  circuit,  and  then  flares  smoothly  to  a width  of  .400" 
at  the  board  edge.  The  centerline  of  the  strip line  center 
conductor  crosses  the  notch  .135"  from  the  short  circuit. 

In  the  vicinity  of  the  notch,  the  stripline  impedance  is  80 
ohms,  and  transistions  smoothly  to  50  ohms  at  the  connector 
junction.  For  the  experimental  investigation,  only  the 
center  conductor  open  circuit  location,  y2-,  was  varied. 

Five  exciter  elements,  differing  only  in  open  circuit 
location,  were  fabricated  and  tested  using  an  HP  Automatic 
Network  Analyzer  (ANA) . The  values  of  for  these  elements 
designated  PI  through  P5,  are  given  in  Table  5.  In  the  4GHZ 
band,  exciter  P5  produced  the  best  overall  match.  In  the  8 
GHZ  band,  PI  gave  best  results. 

Measured  VSWR  for  the  P5  exciter  is  shown  in  Figure  65. 


(b)  Excittr  Board 


: 1 
L , Figure  64.  Baseline  Notch  Exciter  Design  > 


Exciter 

^2 

PI 

0.070 

P2 

0.090 

P3 

0.110 

P4 

0.130 

P5 

0.150 

Table  5 

Open  Circuit  Stub  Lengths  for  Experimental 
Excite rs 
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The  measurement  band  is  4 to  6 GHZ.  From  4 to  4.7GHZ,  the 
exciter  is  poorly  matched,  and  is  operating  in  the  vicinity 
of  its  low  frequence  cutoff. From  4.7  GHZ  to  5.75  GHZ, 
the  mismatch  is  below  2:1.  The  low  frequency  cut-off 
phenomenon  is  known  to  occur  when  the  open  and  short  circuit 
stub  lengths  (x^  and  y2  in  Figure  64)  become  electrically 
short.  Consequently,  by  increasing  the  notch  depth  and  center 
conductor  stub  length,  the  exciter  operating  band  may  be 
readily  reduced  to  the  band  of  interest. 

The  ripple  in  the  well-matched  region  of  the  experimental 
band  may  also  be  largely  eliminated  by  judicious  selection 
of  the  stub  lengths.  In  the  study  of  isolated  notch  antennas, 
it  was  found  that  similar  ripple  in  reflection  coefficient 
occurred  when  the  electrical  lengths  of  the  stubs  were 
significantly  different.  Since  the  spectrum  of  guided  waves 
in  the  region  of  the  notch  short  circuit  may  be  determined 
in  a straight  forward  manner,  the  proper  stub  length  ratio 
is  obtainable  by  either  analytical  or  experimental  means. 

Measured  VSWR  for  the  PI  exciter  looking  into  load 
terminated  half  width  loaded  guide  is  shown  in  Figure  66. 

The  measurement  band  is  7.5  to  8.5  GHZ.  The  match  is  below 
2:1  throughout  the  band.  It  is  evident  from  the  figure  that 
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the  low  end  of  the  experimental  band  is  quite  close  to  the 
low  frequency  cutoff  for  the  exciter  and  that  some  stxib 
lengthening  is  required  to  extend  the  bandwidth  to  cover  the 
full  16%  operating  band.  The  discontinuity  at  8 GHZ  is  due 
to  the  automatic  head  change  on  the  ANA.  While  the  ripple 
shown  in  the  Figure  is  not  large,  some  improvement  can  be 
achieved  by  adjusting  the  stub  length  ratio. 

Probe  isolation  is  determined  by  exciting  a single  high 
frequency  probe  in  the  three  probe  load  terminated  configuration 
and  measuring  return  power  in  the  center  probe.  In  the  high 
frequency  band,  the  measured  result  corresponds,  roughly,  to 
equal  excitation  of  LSEj^O  LSE20  modes  in  the  probe  free 
region  of  the  test  device,  and  consequently  is  an  improper 
excitation  for  high  frequency  operation.  Below  the  loaded 
halfwidth  guide  cut-off  frequency,  the  measurement  approximates 
the  band  isolation  at  broadside. 

Measured  isolation  from  3.5  to  8.5  GHZ  is  shown  in 
Figure  67.  The  center  probe  is  the  P5  exciter  and  is  inserted 
into  the  guide  1.80"  beyond  the  PI  exciters.  The  high 
frequency  exciters  are  butted  against  the  dielectric  slabs 
and  arranged  such  that  the  slab  and  exciter  midplanes  are 
coplanar.  Below  5. 8 GHZ  (the  halfwidth  loaded  guide  cutoff 
frequency)  the  measured  isolation  is  greater  than  40dB,  and 
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exceeds  50dB  from  3.5  to  5.0  GHZ.  As  anticipated,  the 
isolation  in  dB  shows  a roughly  linear  decrease  from  5.5  to 
5.8  GHZ.  Above  5.8  GHZ,  the  isolation  rapidly  degenerates 
to  12dB  at  8.5  GHZ. 

The  weak  isolation  in  the  7.5  to  8.5  GHZ  band  is  due 
to  the  LSEj^q  mode  constituent  of  the  waveguide  field  in  the 
plane  of  the  center  probe  board  edge.  As  was  shown  in 
section  3.2,  the  LSEj^  relative  modal  electric  field  strength 
along  the  loaded  guide  centerline  v/ill  be  on  the  order  of 
.3  to  .4  for  a sleib  relative  permittivity  of  5.  Consequently, 
improved  high  frequency  isolation  can  be  obtained  by  either 
tapering  the  center  probe  board  thickness  near  the  edge,  or 
by  introducing  an  appropriately  located  shorting  stub  along 
the  low  frequency  probe  center  conductor.  In  general,  the 
first  alternative  seems  best  since  it  minimizes  aperture 
perturbations,  but  may  be  impractical.  The  second  approach 
has  several  difficulties  associated  with  perturbations  of 
the  field  at  the  radiating  aperture,  but  is  readily  implemented 
Clearly,  the  high  frequency  isolation  presents  a problem 
which  will  only  be  resolved  through  further  study. 
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6.0  CONCLUSIONS 


The  bifurcated  twin  dielectric  slab  loaded  rectangular 
wave  guide  has  been  shown  to  be  a viable  candidate  as  a 
dual  frequency  array  element  which  results  in  a considerable 
reduction  in  electronic  components  relative  to  other  multi- 
frequency aperture  techniques.  In  comparison  with  wide  band 
elements,  the  bifurcated  twin  dielectric  slab  loaded  rectangular 
waveguide  element  requires  21%  fewer  controls  per  unit  aperture 
area  for  equivalent  scan  and  gain  specifications. 

To  provide  the  multi -mode  aperture  control  required  for 
dual  frequency  operation,  a unidirectional  stripline  fed 
notch  exciter  concept  has  been  investigated  which  results  in 
greater  than  50dB  isolation  between  high  and  low  frequency 
probes  in  the  low  frequency  operating  band.  Preliminary 
exciter  designs  have  been  shown  to  remain  well  matched  over 
greater  than  10%  bandwidths. 

An  element  design  for  operation  over  16%  frequency 
bands  centered  at  4 and  0 GHZ  has  boon  fabricated  and  tested 
in  waveguide  simulators.  Measured  results  are  in  excellent 
agreement  with  theoretical  predictions  in  the  4.0  to  4.32 
GHZ  and  7.32  to  8.08  GHZ  bands.  Measured  results  outside 
these  bands  were  not  obtained  due  to  limitations  of  measurement 


equipment  (below  4 GHZ)  and  computational  practicallity 
(above  8.08  GHZ),  but  are  expected  to  show  similar  close 
agreement  with  theory  in  small  array  testing. 
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APPENDIX  A 

EVALUATION  OF  THE  Epgr,nm  f^OR  RECTANGULAR  LATTICE 


The  coefficients  E „ are  defined  by  the  integral 

pqr/nm 


B+L 1 A/ 2 


(A-1) 


W.nm  “ / dxe^(x,y)-e 

bi  -A/2 


*a 

pqr 


where  A is  the  x dimension  of  the  guide,  B is  the  y dimen- 
sion of  the  guide,  bi  is  the  half- thickness  of  the  septum, 
e^(x,y)  is  an  LSM^  (p  = ')  or  (p  = ")  feedguide 

electric  mode  function,  and  e (x,y)  is  a cell  guide  elec- 

-®pqr 

trie  mode  function  which  is  an  E mode  with  respect  to  the 
array  normal  (r  = 1)  or  an  H mode  (r  = 2) . By  appropriate 
choice  of  ordering  in  both  feedguide  and  cell  guide  regimes, 
single  subscripts  may  be  used  to  identify  the  modes.  These 
indices  are  taken  as  i for  the  feedguide  mode  modes,  and  a, 
for  the  cell  modes,  giving 


B+b 1 A/ 2 

^0,i  ^ f f • ®.*a 

bi  -A/2 


y) 


Expressions  for  e^(x,y)  are  given  in  chapter  3,  and 
expressions  for  e (x,y)  are  given  in  equation  8. 
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In  all,  there  are  eight  forms  of  the  integral 
corresponding  to  the  possible  inner  products  of  feedguide 
and  cell  guide  modes.  These  are: 


1.  LSM-Symmetric  Modes  with  E-Modes 

2.  LSM-Symmetric  Modes  with  H-Modes 

3.  LSM-Anti symmetric  Modes  with  E-Modes 

4 . LSM-Antisymmetric  Modes  with  H-Modes 

5.  LSE-Symmetric  Modes  with  E-Modes 

6.  LSE-Symmetric  Modes  with  H-Modes 

7.  LSE- Antisymmetric  Modes  with  E-Modes 

8.  LSE-Antisymmetric  Modes  with  H-Modes 


Since  the  integrals  possess  many  terms  in  common,  these 
terms  are  defined  first  and  will  be  used  as  simple  variables 
to  shorten  the  expressions.  These  terms  are: 


g sin(kj^+iCj^)  3 
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2.  LSM  Synunetric  Modes  with  H-Modes . 


(A-6) 
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{k  _I1  - k „I!) 
ya  x0 


where  II  and  I2  are  defined  in  equations  (A-4)  and  (A-5) , 
respectively. 

3.  LSM-Antisynunetric  Modes  with  E-Modes 
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(A-7) 
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4.  LSM-Anti symmetric  Modes  with  H-Modes. 
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where  I3  and  I!,  are  defined  in  equations  (A  -8)  and  (A  -9) 
respective ly . 
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5.  LSi;-SyiT\metr.i  c Modes  with  E-Modes 
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6.  LSL-Syinmetric  Modes  with  H-Modes 


(A-13)  E 
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where  1 5 is  defined  in  equation  (A  -12) 


7.  LSE-Antisymmetric  Modes  with  E-Modes 
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8.  LSE-Antisymme trie  Modes  with  H-Modes 


(A-16) 
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where  le  is  defined  in  equation  (A  -15) 
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APPEIfflIX  B 

DERIVATIOij  OF  DIFFEREHTIAL  EQUATION  RELATING  e AND  H 


The  homogeneous  Maxwell  field  equations  are 


(B-1)  V X E (x)  = -jo)pH(x) 


(b-2)  V X H (x)  = jcoeE  (x) 


To  separate  out  longitudinal  components  (z  directed)  take 
vector  and  scalar  products  of  (B-1)  and  (B-2)  with 

I-.-Q 

Thus 


(B-3) 


jwyH(x)xz^  = ^x(VxE(x))  = 


E(x)  + VE^(x) 
d Z’  — — z — 


fi  (Vz'2  * ^‘2))  * <V2o 


- f;  S (5) 


(B“4) 


^ • (VxE(x))  = -V^  • (Z^xE(x)  ) 

where  e(x)  is  the  tranverse  to  z electric  field. 

Simi larly , 

(B-5)  X E(x)  = V^H^(x)  - 1^  h(x) 

(B-6)  jcoeE  (x)  = V * (HxZ  ) 

Substituting  for  E (x)  in  (B-3)  from  (B-6)  gives 

2 

j,  A V . A 

(®-7)  - ^ e(x)  = ja)u(-I  + ] • (h(x)  x 

Recognizing  that 

e(x)  = e ^^^e(x,y) 

/\  < 

h(x)  = e ^^^h(x,y) 


for  uniform  (in  z ) media,  (b-7)  reduces  to 


(B-8)  ^ e(x,y)  = t^  + -^T^l  • (h(x,y)  x ^} 

And  from  C^4)  and  (B-5) 

V.  V 

(B-9)  Yh(x,y)  = we  [*I  + ] • (^xe(x,y)) 


For  the  inhomogeneously  filled,  uniform  in  ^ wave- 
guide, uncoupled  modes  will  be  either  LSE  (e^  = 0)  or  LSM 
(h^  = 0) . For  LSE  modes,  solution  of  (B-8)  gives 


(x,y)  = 


- h”(x,y) 

k ^-kf.  ^ 


e^  is  related  to  h^  by  an  impedance.  It  is  convenient, 
therefore,  to  define  a modal  admittance  Y’’,  such  that 


e^(x,y)  = -h"(x,y) 


and 


//< 


dA(h.xz  ) 
—1  -o 


• e . 
-3 


In  this  manner,  (b-8)  may  be  rewritten  as 


(B-10)  YZ"e"(x,y)  = wp  U + -^]  • (h"(x,y)xz^) 

Similarly 

V V 

(B-11)  YY'h'(x,y)  = wefi  + -^]  • (z^xe'(x,y)) 
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APPENDIX  C 

ORTHOilORHALIZATION  OF  FEEDGUIDE  I10BE  FUiKTIOUS 

Let  L , li  and  E ,H  be  linearly  independent  character- 
jstic  relations  of  Maxwell's  source  free  equations  in  the  cyl-- 
indrical  guide  shown  in  Figure  D-1.  The  guide  is  uniform  in 
z,  resulting  in  z dependencies  given  by 

- j Y z ~ j Y z 

e , e 

Assuming  a lossless  region,  then  the  conjugates  of  the 
characteristic  solutions  are  also  solutions  of  Ma>well's 
Equations.  Consider  the  curl  equations 

(C-Ia,b)  VxE^  = 

(C“lc,d)  VxH^  = j'^'^-n 

where  * denotes  conjugation.  Taking  scalar  products  of 
(C-la,b)  with  H*  cind  respectively,  and  adding,  cjives , 

(c-2)  H*  . VxE  + II  . VxE*  = 0 

-m  — n -n  --m 

.Similarly,  from  (c-lc,d) 
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-rr. 
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(C-3) 


E* 
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VxH*  = 0 

-TO 


Adding  (C-2)  and  (C-3)  and  manipulating  results  in 


(C-4) 


0 = . (E*xH  + E^xH*) 

t -TO  — n — n —TO 


(E*xH„  + E^xH*) 

—TO  — n — n —TO 


where  is  the  transverse  gradient  operator,  and  the  longi- 
tudinal gradient  has  been  specifically  shown.  Since  the  entire 
z dependence  is  embodied  in  the  exponentials,  (C-4)  becomes 


(C-5) 


V.  . (E*xH„  + E xH*) 
t — m — n ~n  — m 


= • (^t  ^ »t  ^t  ’"^^t  ^ 

m n n m 


where  the  transverse  field  is  indicated  by  the  subscript  t. 
In  equation  (05)  the  longitudinal  components  of  the  fields 
have  been  ignored  on  the  right  hand  side  due  to  the  Zo  . 
operator.  Let 


( 06) 


(x,y) 


z 


and  similarly  for  the  other  explicitly  transverse  field 
quantities  (subscript  t) . Then,  application  of  the  diverg- 
ence theorem  (z-dependent  integrals  cancel  out)  results  in 


(C-7)  (Y. 


n’^m>  /A*  • 


(e*  X h + e X h*)da 
-m  -n  -n  -m 


= 0 


where  S is  the  cylindrical  cross-section  and  da  is  the 
differential  unit  of  transverse  area. 


The  fields  E , H / E , and  H are  assumed  to  be  linearly 
— n — m — in 

independent.  Hence  Yj^  ^ Yj^/  and  (c~7)  may  be  rewritten  as 


(C-8) 


/A.. 


(e*  X h + e X h*)da  - 0 
— m — n — n — m 


Since  the  direction  of  propagation  (±z)  should  not  effect  the 
result  (c-8) , consider  the  z dependencies 


“DY  z +1Y  z 

n ^ m 
6 / 0 


Then,  by  entirely  equivalent  steps, 


(C-9) 


(Y  +Y  ) 
'n  'm 


IL 


(e*  X h„ 

—in  — n 


e X h*)da 
-n  -m 


0 


or,  for  Y 7^-Yn,/ 
n n' 


(c-10) 


lU 


(e*  X h - e X h*)da  = 0 
— m — n — n — m 


Adding  and  subtracting  (C-8)  and  (09)  results  in 
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fC-ll) 

//-. 

(h 

-m 

-n 

(c-12) 

f.k' 

— m 

£o)da 


0 


0 


Equations  (C-11)  and  (C-12)  are  the  desired  mode  orthogonality 
relations . 

In  the  instance  that  then,  assuming  the  eigen- 

values are  not  degenerate,  equation  (C-7)  is  satisfied  indepen- 
dent of  the  value  of  the  integral.  From  equation  (C-lO),  with 
m n 


(C-13) 

e*  X h 

X h* 

— n — n 

— n 

— n 

Substituting  (C-13) 

into 

(C-14) 

//-■ 

(e* 

X h + e 

X h*)da  = constant 

V - 

— n 

— n -n 

— n 

yields 

the  desired 

mode 

function 

normalization  integral 

(C-15) 

IL  . 

iK 

X ^o)do  = 

cons  tant 

If  the  constant  in  equation  (C-15)  is  taken  as  1,  then 


f 

! 

J 


I 


in  the  field  representations 


E 

H 


Z V.e. 
1—1 


Z I.h. 
: 1—1 


has  units  of  power. 

Using  the  expressions  for  and  h^  (r=',")  from  section 

4 in  equation  (c-15)  with  unit  constant,  and  rearranging 

results  in  the  following  integrals  for  the  normalization 
r 

constant  N. : 


1.  LSM  symmetric  modes  (r=') 

A/z  b 


where 


(N 


, S) 

nm 


dx  / dy 
-A/2  'O  r 


1 


t (>:) 


• ? m i',  y I , s , . 
!in‘'-r;^-i  I'  U) 
D n 


= bf:;?  (2-c'  •■)  (2k'  ) 

1 2 I C ^ 1 Z ) 

+ i-|8'  |‘yS  {2r'A)n(>  ') 

' ? Z 7 t r 


-t  (>:' ) 8' 8’ (K't)r(K') 
t 1 ? •)  f. 


+ |C  I (2'.'’  u)| 


cos.: ' 8 
r r ' 

= — ^1 — sint  ' 

f 

i un  ‘ 


(”  ~ — ; ; — cosK  rS  t rcot.  , sim 

s.in.  £ f ! k'  ( 

r 
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r-w  ...  ..i— 


and 


t(x)  = 


o(x)  = 


Si(x)  = 


SaCx)  » 1- 


1, 

x=|xl 

j» 

x=-‘j  1 X 

1, 

x=lx| 

-1, 

x=- j 1 X 

sirx 

X 

+1 

, sinx 

X 

t (x) 

sin^x 

X 

2.  LSM  antisymmetric  modes  (r=M 

A/2  b 


(N'^) 

nm 


-A/2  0 ^ " 

= b||S2(2K'6)  + |-|E'i|  ' 

+i_lEJ|2|s2(2K;.6'a(K;.) 

■Pi.l5.y-F5S,  (K'  6) 

C-  0 

r 

+ |F'  |2|Si(2ic'a)J 


E;  = sinK'g 


G_K 

E'z  — cosk'6 


where 


N)|0. 


1^ 

pi  = 22S£li[cosK' 6 ^-^-tanK*  8siriK^$] 


sxriK  a 


F = 


I Re{  ElEk*},  k'  = 


Im{ElEi*},  <1  = -3\<' 


3.  LSE  symmetric  modes  (r=") 
A/2  b 


(n:®) 

nm 


dx  jdycos^™|V"®(x)  1 2 
-A/2  0 

rmb{|Si(2K"e)  + S'l’^fSi  (2k^6) 

+ |B'l|^fs2(2K"6)a(K';) 

-t(K")6l'8’^6S3('c"6)a(K") 

G 0 t 

+ |C"  1 2|s2(2K”a)a(K")| 


where 


r = Neumann  factor  = 
m 


1,  m/0 

2 , m=0 


B"  = cosk"6 


Bil  = •^sinK"6 
C" 


2?^{cosk" 
sinK  a E 


- ^tan<  Bsinx  6] 
ic  e 

e 
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4.  LSE  antisymmetric  modes  (r=") 
A/2  b 


(N"5 

nm 


-A/ 2 


dx  J dycos^S^I  V^^(x)  | ^ 


rmb|fs2(2K"6)a(ic")  + 

+ |E2‘|"|s2(2!C^6)a(icp 


- F6S3 (k^6) 


+ |f"'  rf5  2(2K*'a)| 


El'  = -sinK"6 


F" 


^osk"6 

e 

SinK"gr  _ .ir  . k"  . ..  ^ ■ 

+ pTCOtK”8sinK^6 

e 


F = 


EiE'^, 

Re{E'iEH*}, 


k"  =k" 


-3\K' 


i(2k^6) 


j0.5{El'-Ei'*}E'k 


k"  = -j  1 k" I , 


APPENDIX  D 
PROGRAM  LISTINGS 

This  Appendix  gives  listings  of  all  programs,  and 
subprograms  required  to  reproduce  the  numerical  results 
presented  in  this  report.  In  general,  the  listings  are  self- 
explanatory.  The  language  is  FORTRAN  (extended)  and  the 
programs  are  designed  to  run  on  CDC  6600,  6700  and  CYBER  73 
series  computers. 

The  Appendix  has  two  subsections.  The  first  gives 
listings  of  main  programs.  The  second  section  gives  listing 
of  subroutines  and  function  subprograms  required  for  execution 
of  the  main  programs. 
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ooo  oooo  non  o o rs  o o o o o n 000000 


PROGRAM  PlTPAT  CINPUTfOUTPUT* TAPE5»INPUTf TAPE6suUTPUTfTAPt«) 

computation  of  BIFURCATEO  T^I'^  niELECTRIC  SLAB  LOADED  RECTANGULAK 
wAVEGOlDE  ARRAY  LLEMENT  RAdIATJON  CHARACIEHISTICS, 

REUUiRtS  CALCOMP  library  AnO  SUBROUTINES  GIVEN  IN  APPENDIX  0,2, 

dimension  powdOfbiitSTisn ,sii(2) tS2i(2) tMORocao) tipao) tiunoi , 
UR(I0),JJ(10) 
dimension  IbUFUOOO) 
real  KtKEfL 

Integer  Pi»QifTu(ioj»yo(io),Ro(io)fSiGono)isiG(2o),siGio(io) 

complex  Y,YA»L10»C20»C1  »C2,^ J,C.3,R1  .R2fCa 

Complex  sii  (?ot2oj  ,s2H25o»20)»vo 

common  /ARRAY/  AL » BL f OL t P I ♦ SEPTl , TPI ,EPSf SI (2) , S2C2 ) 

Common  /cnSrv/  y(20) , yac2S0) 

common  /MODES/  K(20) ,^E (20) ,GAMwa(20) fMODEl (20) » ISYM(20) »NNC20) ♦ 
IMm(20) »MOOORD(20) 
common  /PO/  Pl»ul 

data  TPI*C/b.2S31853071  7R6»i  I ,ft(j28526/ 
data  AJ/(0,«1  .,)/ 

CALL  plots  (IBUFf lOOOtA) 
knt*o 

INPUT>  FREQUENCY  IN  GHZ,  fOA  RELATIVE  SLAB  PERMITTIVITY,  EPS# 
FEEDGUlOE  HEIGHT,  RD,  IN  UjCHESa  A(ALPHA),  B(HETA),  AND  D(DELTA), 
IN  iNChESA  and  aperture  PLaN  Tu-REFFRENCE  PLANE  SEPARATION  IN 

feeoguide  wavelengths,  ALEn, 

READ  (S*0OO)  FO*EPS,BD,A,B,DtSFP,ALEN 
IF  (EoF(5),NE,0)  go  to  160 
Al«2,A(A+BfD) 

write  (6,900)  F0,EPS*A1  ,BD,a,B,D,SEP,ALE.N 
INPUT>  limits  OF  sampled  grating  lobe  space,  pi  AND  01 
REAL)  (B,02O)  PI, 01 

INPUT>  lattice  vectors  Si  aND  S?,  IN  INCHES, 

read  (S,800)  Sl,S2 
DO  10b  I«l»2 
SlI(I)*Sl (I) 

82l(I)«S2(I) 

CONTINUE 

igrd«i 

SsSl(l)A82(l)+Sl(2)AS2(2) 

IF  (AbS(S) ,LT,1,E-10)  IGRD*? 

0X*S1 (1 ) 

DY»82(2) 

INPU7>  number  of  modes  to  establish  ordering,  NMODE'1,  #MAX  of 
20AA  and  frequency  hand  DEsI&MatiON,  LOHIAsLO  or  HIA, 

READ  (b,810)  NMUDE1,L0HI 
IhLO*LPHI 

IF  (8EP,G7,1,E-10)  IHL0S2HL0 
INPUT>  SINE  SPACE  SCAN  RANGE  AND  INCREMENTS, 


READ  (5*800)  3THS»8THfc*8THI,8PH8*8PHE*8PHl 
NTH«INT((8THE*81H8)/8THI+0,5)+1 
NpHilNT((8PHE-8PH8)/8PHl+0,5)+l 
C 

C INPUT>  number  op  FEEOCUIOe  MODE8  TO  BE  USED  FOR  APERTURE  FIELD 

c approximation*  nmooesp  and  number  of  beams  to  be  plotted*  lobes. 

c 

READ  (5*620)  NM00ES*LUBES 
NMvNMODEStl 
DO  106  laltNMOOES 
C 

C JNPUT>  NM0DE8  FEEDGUlDE  MqdE  DESIGNATIONS  PE.G.*  LSEOJOO* 

C tSMOlOl*  L8E0201P. 

C 

READ  (5*830)  MORD(I) 

106  CONTINUE 

IF  (LOBES.EU.O)  60  TO  108 
IF  (LOBES, GT, 10)  STOP  'LOBEsaIO* 

C 

C INPUT>  PLOBEBP  beam  DESIGNATIONS  IN  FORM  P*U*R 

C 

DO  107  Ial*L0BES 

READ  (5*820)  TO(I)*UO(I)«HO(I) 

Sl60(I)«00(I)+Ql  + U(T0(I)^Pl)»(2*Ql  + l)  + (R0(I)-l)»(2APl  + l)* 
l(2«Gltl) 

107  CONTINUE 
C 

C INPUT>  complex  LSE0200  mOdE  VuLTAGE  MuDIFIEK*  Rl#  AnO  COMPLEX 

c upper  element  half  VOUAGF.  modifier,  R2, 
c 

108  READ  (5*800)  Rt*R2 
write  (6*900)  R1,R2 
write  (6*910)  SlfS2 
write  (6*915)  P1*UI*NM0UES 
L«FO/C 

al«a*l 

bl«bfl 

ol»d»l 

BlBbOFL 

SePTL«8£P9L 

C 

c compute  dispersion  relation  at  center  band. 

C 

CALL  L8MLSE  (NMOOEl) 
alen«alen/(l*gamma( j ) ) 

IF  (ALEN,GT,1,E-05)  AlEMsAlEW 
AlENbaLENI 

C10«CEXP(AJATP1*GAMMAIU*AI.eN*U 

C20«CEXP(AJATPI*GAMNA(2)*ALfN»L) 

c 

c IMPUT>  half  bandwidth  (INTEGER)  IN  PEHCENI*  iBWp  AND  PERCENT 

c bandwidth  STEP*  ISTEP, 

C 

READ  (5*820)  iBwflSTEP 
N8w«2*IBWt1 
C 

C COMPUTE  performance  over  BaNL, 

c 

00  155  J8w«lfNBw*I8TEP 

0W»1.tO,OIF(J8W-1*IBwJ 
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F»8W*F0 

IF  (iB^tEU.U)  60  TO  1U4 

write  (6t<»05}  8y<*F»EPS«Al*BD«A»Bf0f8EPtALEN 
109  L«F/C 
AL«A91. 

BL«BH 

dl«o»l 

BiabDAL 

septi«sep*l 

Bll«0,5»(Bl+8EPTU) 

IF  (IBW.EQ«0)  60  TO  1091 

Call  lsmlse  cnmoded 

1091  Cl«C10FCEXP(-AJ»TPI*6APMA(l)FALEN*L) 
C2»C209CEXP(-AJ*TPI*GAWMAt2)FAUfcN*L) 

DO  MO  1«1*2 
81(I)«S1I(I)«L 
82(n«S2I(n«L 
MO  CoVTIr^UE 
JiO 

Do  M6  iBltNMOObS 
MS  JaJfl 

IF  CM0D0RD(J).NE,M0R0(I))  60  TO  H5 
K(I)«K(J) 

KEa)«K£(J) 

6AMMA(I)aCAMMA(J} 

NN(l)aNN(J) 

MM(l)a^»M(J) 

MOOEl(I)aMODEKJ) 

ISYM(I)btSYM(J) 

MQDORdM)aMORD(I) 

M6  CONTINUE 
iWRITeai 

CALL  NORM  (NM00E8I IWRITE) 

8PH»8PH8*8PHI 
DO  150  IPHbIiNPH 
KNTaKNTtl 
00  M60  iBltlO 
DO  M63  ITHal,5l 
PQWUf  ITH)B»100#00 
M63  CONTINUE 
M64  CONTINUE 

8PHaSpHf8PHl 

CPHb80RT(1,»8PHAF2) 

STHB8TH8-8THI 

IF  aHL0,EQ,2HL0)  WRITE  (bf920)  SPH 
IF  (IHL0.EQ.2HHI)  WRITE  (6»q30)  8PH 
Jl«0 
KJlaO 
C 

C TAKE  theta  CUTS. 

C 

DO  HO  ITHBl.NTh 
STHB8TH+8THI 
8TCITh)«8TM 
CTHb80RT(1.-8TH»»2) 

Ua8T;i9CPH 

Vb8TH»8PH 

C3«CEXP(AJ»TPHV*BM) 

C4«R2/C3 

C3«R2*C3 


or»r»r»r»  •-»  oooooo  o o o 


EVALUATE  SCATTERING  BLOCKS  811  A^0  S2l 

CALL  SCTMAT  (XHLOtNMODEStUf VtSl 1 «S21 tN^OOf ISIG) 

JaJl 

KJBKJl 

00  117  laltlSlG 

IF  CREALCYAtm.Lt.O,)  GO  To  117 
IF  (Jl.EQ.O)  GO  TO  1162 
DO  1161  Tlslfjl 
IF  (I.EQ.SI&am  GO  TO  117 

161  CONTINUE 

162  J«J7l 

IF  (J.LE.IO)  GO  TO  1168 
KjaKJfl 
SI610(KJ)>I 
165  S1G(J1>1 

17  CONTINUE 
Jl>J 
J2*J 

IF  (J2.GT.10)  J2alO 
KjlaKJ 

IF  (L0HI.EQ.2HHl)  GO  TO  120 
Pa2,FCABS(Y(l)) 

PoaO. 

COMPUTE  POWER  TRANSMISSION  COEFFICIENTS  EOR  FIRST  10  PROPAGATING 
beams  in  GRATING  LOBE  SEQUENCE*  aT  LOW  FREQUENCY  AND  POwtR 
REFLECTION  COEFFICIENTS.  paINPUT  POwER#  PO»SUM  OF  REAL  PUwER  IN 

all  beams  and  fcedguioe  modes. 

00  US  lalfjl 
iSGaSlG(I) 

PTa(CABS(S21(lSG*l)fS21 (ISG,NM))a*2) AREAL (YA(ISG))/P 

SaREAL(YA(IS6))-1.0 

iPTal 

IF  (ABS(S).LT.1,E-10,AND,ABs(SPh*CPH).LT.1.E-10)  IPTaiPTAl 
IF  (ABS(SPH-1,),LT.i.E-10.AnD.I3G,GT.ISIG/2)  IPT«IPT-2 
IF  (ABS(CPH«1,),LT.1.E-10.An0.ISG.LE.ISIG/2)  IPTalPT-2 
IpTalABS(IPT) 

PTalPTFPT 

POapO+PT 

IF  (PT.LT.l.E-iy)  PTal.E-10 
IF  (1.6T.10)  GO  TO  lie 
POW(IfITHjalO,AAL06lO(PT) 

18  continue 

PRUa(CABS(Sll  CUntSlKl  tNM))*A2)ACA83(Y(l))/P 
POapOtPRU 

IF  (PHU.LT.l.t-lO)  PRUal.E-iO 
PRUalO.*AL0610(PRU) 

PRLa(C^BS(Sll(NM«i)tSll CNM.NM))AA2)acA0S(Y(l))/P 
POapOtPRL 

check  CONSERVATION  OF  tNERQY,  JPO  AND  IM  AR£  MINUS  THE  NUMBER  OF 
digits  to  WHICH  CONSERVATION  OF  ENERGY  IS  APPROXIMATED  BY 
SOLUTION. 


P0aALOG10(ABS(P0«1.0)7l.E«60) 

IPOaPo 

IF  (IP0.LT.-99)  IPOa-99 
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IF  (PRU.LT.l.E-lO)  PRL«UE*10 
PRL^lO.^ALOGlOCPHU) 

CALL  CONSHV  (SUiS^ltNMOOf  tst&»lM) 

write  (6*9«01  SlH»PRU»P«L»TPft»lN'»(PUW(I»nH)»l«l,J2) 

GO  TO  130 

compute  power  TRANSMiSIOiVi  aMD  rffuection  cueepjcients  at  HIG 
freouencv, 

V2»1An(0,25*TPIAS1(1)*u) 

V0*RUV2 

V2«CAUS(V0) 

P«CABS(Y(l))ACAbS(V(2))AV2**2 
IF  (IhLU.NE.LOHl)  60  TO  12^ 

l>  INFIN1TE88IMALIY  thin  SfPTUM  ANU  RECTANGULAR  GRID, 

P0«0*0 

no  121  i>i«ji 

I8G«8X6(I} 

PT«HEAL(YA(I8G) J*CCAS8(82l(l8Gf nTS2l (I86»2)*V0)*»2)/P 

8>REAL(YA(I86))*1,0 

1PT»1 

IF  CAB8t8J.LT,l,E-lO.ANU,AB8(8PH*cPM),LT«l.b-lO)  IPT«IPT+1 
IF  lAB8l8PH*l,).LT.Ut-10,ANO.lSG,GT, 1816/2)  lPTalPT-2 
IF  (Afl8tCPM*l,) ,LT.1«£-10,AnD,ISG,LE.ISIG/2)  IPTslPT-2 
IPT«IA88(IPT) 

pt»ipt*pt 

PO«PO+PT 

IF  (PT,LT,1,E-10)  PT*l.flO 
IF  CI.GT.IO)  GO  TO  l2l 
POWClf ITMJ«lO,FALUGlO(PT) 

CONTINUE 

PRlaCA0S(Y(n*(8ll  (1 1 DtSll  (1  »2)*vO)»Y2)/P 
po»potpki 

IF  CPRl.LT.l.E-lO)  PRl»l.e-lO 
PRl«l0.*ALOGl0(PRl) 

PR2«CAB3(YC2)»(SU(2»1)+Sn  (2f2j*vO)*A2)/P 
P0«P0+PR2 

CHECK  C0N8ERVAT1UN  UP  EN’ERgV. 

P0*ALOG10(AB8(PO»l,0)tl.E-80) 

IPO«PO 

IP  (IP0.LT.-99)  IP0*-99 

IF  (PH2.lt. I.E-IO)  PR2al,E-lO 

Pr2«10.»ALOG10(PR2) 

CALL  C0N8RV  (SUi82l»NMOD»l8l6flM) 

write  (6f9«0)  8TH»Pr1  fPH2,IpO»lM,(PUw(IfITH)fI«i.J2) 

GO  TO  130 

CASE  2>  Thick  septum  or  triangular  grid, 

Pb2.09P 

P0«0.0 

DO  125  IsltJl 
I8G«81G(I) 

PT»KEAL(YA{lSG))AtAeS(((S2l(ISG»n»Ca  + S21(l8G»NM)*C3)FCl  + 
l(82l(IS6»2)*C4tS2l (lSGfNM+i)*c3)AC2*vO)A*2)/P 
8«REAl(YA(ISG))-1.0 
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IPT»1 

IF  (AB8(8).LT,i.E*lO,A^O,A»sC8PH*cPH).LT.U£«lO)  IPT«IPT  + l 
IF  CABS(SPH•1,),LT.1•E•10,A^D,I8G. 67,1816/2)  |P7»IPT*2 
IF  (AbSCCPH*!,), IT. Ufc-lO, AMD, lS6,Lt, 1816/2)  IPT«lPT-2 
1P7«IAB8{1PT) 

PT*1PTPPT 

PO»PO^P7 

IF  (PT.LT,l.r-10)  PT«l.E-10 
IF  CI.GT.iO)  GO  TU  J25 
POWlI»ITH)»lO,*AL06lO(PT) 

125  CONTINUE 

PRu*CABS(Y(i)PC(sn  (i»i)*c^+sn  (1  ♦nm)*c3)*ci+(su  (ii2)*c« 
i+sll  (1»NM+U*C3)*C2«V0)*P2)/P 

PRUsPKUtCA8S(Y(2)P((Sn  (2fl)»C4+SlU2»N(^)*C3)»Cl*t8lK2»2)*C« 
UsnC2»NM+n»C3)*C2»V0)**2)/P 
P0«PO+PRU 

IF  (PRU,LT, 1,6*10)  PPU«l,fc*iO 
PRU«lO.AAL06lO(PRU) 

PRt«CAB8(Ya)4((8ll  (NM,n*Ci»+8U  (mm,nm)*C3)»Ci  + (SU  (NM*2)»Ca 
H8U  (NM*NM4-1)»C3)»C2*V0)*»2)/P 
PRL*PRL  + CABS(Y(2)*(  (SlKNM+i,  1 )*C«*S11  (NM+1»NM)*C3)*CI  + 
l(Sll(NMtlf2)*C4f811(NM+ltNM+l)»C3)*C2*V0)**2)/P 

P0»P0+PKL 

C 

C check  C0N8ERVATIUN  OF  ENERGY. 

C 

P0*ALO6l0(A88(P0*l,0)tl.E*B0) 

IPO«PO 

IF  (IP0,LT.*P9)  IP0«-99 

IF  (PRL.LT,l,e*10)  PRL*l.E-lO 

PRL«10.*ALOG10(PRU 

Call  C0N8HV  (SU»S2l*NM0Dfl8l6»lM) 

write  (6»9a0)  8TH»PRU»PRLfIp0flM,{PUW(I»lTH)»I»l»j2) 

130  CONTINUE 
190  continue 
C 

C PRINT  beam  designations. 

V 

lQl«2>YUlf  1 

IS6SISI6/2 

DO  ;9l  I*lfJ2 

DO  1«10  Il«l»' OBES 

IF  (816(1), NF.,SIG0(IU)  GO  TO  1910 

Jv)(ll)8l 

GO  ro  1411 

1410  continue 

1411  lR(I)sl 

IF  (8IG(I),6T,ISG)  IR(I)»2 

IP(I)«(SIG(I)-l-aR(I)-l)Y‘ISG)/10l*Pl 

Ioa)*Sl6(I)-l-Ul-(lRCI)-l)*ISG-(lP(I)+Pl)»iai 

141  continue 

WRITE  (6»950)  (IP(I)fIQ(I)flRa),I«l»J2) 

IF  (KJ1,EG.0)  go  to  143 
DO  142  IclfKJl 
IR(I)*1 

IF  (8IG10(I),GT,ISG)  1R(I)«2 

IP(I)*{8I610(I)*1-<IR(1)*1)*I8G)/101*P1 

IO(I)«SIG10(I)-l-01*(IRn)-i)PlS6*(IP(I)+Pl)*I01 

142  CONTINUE 

WRITE  (6.960)  UP(I).IQ(I).lR{I)»I«l.KJi) 
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143  IF  (LOBbS.EQ.O)  CO  TO  ISO 

IF  (KNT.CQ.I)  call  plot  (0o,S««3) 

C 

C plot  PLOBE8F  BEAMS 
C 

call  plot  (ti.f.2f2) 

ISTRTtl 

L0BE«MIN0(4,t08tS) 

CALL  PLTCAL  (A«B«OfF*EPS*Al,BOtD)(»OY,8EPfPOMfSlGOyLOBE*lSTRTtJJt 
18TfSPH*lGRDfNTH) 

IF  (LOBES. LE. 4)  GO  TO  150 
CALL  PLOT  (0.*.2t2) 

L0BE«MIN0(4yL0BtS«4) 

iSTRTiS 

CALL  PLTCAL  (AfBfUtFfEPS«Al.BD*i)XtOTtSEP*POIf«fSlGOtLQBEtlSTRT(JJt 
l8T»8PHflGR0fNTH) 

IF  (LOBES. LE. 6)  60  TO  150 
CALL  PLOT  (0.».2t2) 

LOB£*MINO(2»LOBES-a) 

I8TRT«9 

CALL  PLTCAL  ( A f 8»0yF  yEPSt  A 1 ,BO*DX yOY  »SEPt POi^ tSlCO »LOBE ♦ ISTRT  ♦ JJ ♦ 
iSTtSPMiIGROfNTH) 

150  Continue 

155  CONTINUE 

Go  TO  100 

160  Call  plot  (o.f.2«2) 

Call  plot  (io.»o.t<>99) 

ENDFILE  4 
CALL  EXIT 
C 

800  FORMAT  (8F10.0) 

810  FORMAT  (I5»A2) 

820  FORMAT  (1615) 

830  format  (A7) 

900  format  (lHlf46Xy37HDUAL  FREQUENCY  ARRAY  ELEMENT  PATTERNS »/// i 

l5Xf 13MELEMENT  DATA&y/i 1 OX ySHFO  s ,F5.2»5Xy6MEPS  » »F5.2»5X»4HA  s y 
2F5#3*5Xy«H8  « yP5y3y/y lOXySHACPHA  * yF5,3y5Xy7HHETA  » yP5,3y5Xy 
38H0ELTA  a yF5.3y5X y9H8EPTuM  a yF5.3y5Xy 1«MALEn/CAMRDA  s yF5,3y/) 

904  format  (5Xyl8HV0LTAGE  MOn IFjERS&y / y I OX y 5HR 1 s y 2F 7 , « y 1 H J . b*  i 
15hR2  a y2F7.4ylHjy/) 

905  format  (IHI  ylOXy7HCASt&  y F<j .2 * 2HF0 y // y 1 OX y 4Hf  s yF5.2ybXy 
16HEPS  a yFS.2y5XyAHA  ■ yF5.2f5XyaMB  a y F5.3 y / y 1 OX y BHAlPHA  * y 
2F5t3t5Xy7HBETA  « yF5.3ySXy8HDELT A « yF5,3y5Xy9H5EPTUM  a ,F5.3y 
35Xy7HAL£N  a yF5.3y/) 

910  format  (SXyllHAHRAY  DAT AA y / , 1 OX y anSl  a y F 6 . 3 y I M y , F5 , 3 f 5X y AhS2  sy 
lF6.3ylHyyF5.3) 

915  Format  (lOXyAMPl  ayl3y5Xy4MQl  8 , j 3 y 5X y BHNMOUES  ayHy/ylOX, 

12AHALL  DIMENSIONS  IN  INCHES,/) 

920  Format  (IHlySXyUHSlNCPHI)  • yF5.2y///ylXy7HSIN(Th)  ySXySHPRUySX, 
13HPRLy7Xy3HIPOy2Xy2HlMy3Xy3lHPO^ER  IN  EXCITED  BEAMS  (OB)  —Ay//) 
930  FORMAT  (IHl y5XyllH8lN(PrtI)  S yF5. 2 y /// y I X y 7hSlN ( TH ) y Sx y 3HPK  I y 5X y 
13HPH2y7Xy3HIPOy2Xy2HlMy3Xy3)HPO^EH  IN  EXCITED  BEAMS  (D8)  —Ay//) 
990  format  (2XyF5,2y4Xy2(iXyF6.2ylX),3Xy2ClXyI3ylX)yl0(lXyF6,2ylX)) 

950  format  (/ySXySHPyQyR  a y 26X , 1 0 ( 1 X y 12 y 1 H y y I2y iHy y I 1 ) ) 

960  Format  C//y5Xy2bHSPACE  MOOES  NOT  PRINTED  ARE4  y / y 5X  y 8HP  y Gl  y k a y 
ll0(2(l2ylHy)yIly3X)) 

END 

•• 
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n 


PROGRAM  OESMOO  (INPUT»0UTPUt*TAPF5»INPUT»1 APEb«0UTPUTiTAPE7) 

computation  of  first  four  routs  of  symmetric  and  ANTI-SYMMtTRiC 
LSE  and  L8M  MODE  DISPERSION  RELATIONS  WitH  M«0  FOR  TwIN  DIELECTRIC 
SLAB  LOADED  RECTANGULAR  aAvEGUIOE  IN  RANGE  I .LE #K«A/2,AND. 
K*A/2,LE.a, 

CREATES  DATA  FILE  INPUT  FOR  PROGRAM  UESGN. 

DIMENSION  BBU6J«G(16) 

REAL  K(l6)fKA2 

COMMON  /WAV60/  AtbtOfBl fTPI,EP8 
data  CPI fTPI/3. 756964667 *6. 2831 65308/ 

INPUT>  A(ALPHA)»  B(BETA)»  AND  D(D£UTA)»  IN  INCHES#  GUIDE  HElGHTf 
Blf  IN  INCHES#  AND  RELATIVE  SLAB  PERMlTTIVITVt  EPS. 

00  read  (StSOO)  A«BfD«Bl«EPS 
IF  (EOF(5l.NE.Ol  call  EXIT 
AAs2«#(AfBfO) 
write  (6t9l0) 

write  (6*900)  AfBfDiBltEPS 
write  (7*900)  A*B*D*B1*EPS 
EPSU«SQRT(EPS-1.0) 

DO  IlO  Isl*16 
BB(I)«-£PSQ*UE«'0 

no  continue 

Ka2».97 

DO  I30  IslflOl 
Ka2*Ka2>0.03 
FSCPI9KA2/AA 
CALL  F0URMD8  (F*B8) 

Do  120  J«IfI6 
K(J)«KA298B(J) 

TaKA2#«2-K(J)«ABS(K(J)) 

G(J)>8QRT(AB8(T)} 

IF  (T.LT.O.)  G(J)«-G(J) 

BB(J)>BB(J)*1.0 

IF  (BB(J).LT,-EP80)  BB(J)»-ePS0^1 ,E-10 
120  CONTINUE 

WRITE  (7*900)  KA2 
WRITE  (6*900)  Ka2 
WRITE  (7*900)  (K(J),J«1*16) 
write  (6*900)  (G(J)*J>1*16) 

UO  CONTINUE 

GO  TO  100 

c 

600  FORMAT  (6F10.0) 

900  format  (SFIO.6) 

910  format  (IHl) 

END 

$$ 
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PROCRAM  DESGM  ( input  t out PUT,TAPt5*INRUTfTARt6sfJUTFUTfTAPhl0j 


COMPUTATION  OF  BROaOSIDE 

array  of  bifurcatfd  twin 

WAVEGUIDES  VS.  K«A/2. 
which  CREAIES  TAPEIO. 


Scan  element  i«ibMATcH  fur  Infinite 
OiFLECTRIC  slab  LUAuFD  HECTanGI'LAH 
use  In  CONJUNCTiOiv  WITH  PRUG«A*^  OLSmOU 


dimension  NM(10I  »NE(10)*T1  (16)  f T (6«)  »IM0Dt(2)  fPT  Ci> )♦  P TDb  f 2 » 1 0 l ) i 
lPRU2)»PRlOB(2)fSSl  (2J»SS2(2) 

REAL  KtKEtKAUOn  iMgS 

integer  PlfUl»SVMC2j 

COMPLEX  YiYAtSU  (20,20)  .821  (250, ?o) 

COMMON  /ARRAY/  AL , BL »UL ,«1L , StPI L , TP  I , EPS ♦ SI (2) , S2 (2 ) 
common  /CNSRV/  Y(20) ,YA(250) 

common  /MODES/  K (20)  ,KE  (2u)  .gamma  (2u) , MODE  1 (20  ) . ISYM  (20 ) ♦ NM20  J , 
IMm(20) ,M0D0R0(20) 
common  /PCiI/  Pl.Ul 
DATA  TPI/6.2831B5i0717<)6/ 
data  IM0DE,SYM/3HLSM,iHL8t« IHA, 1H8/ 

iwrite«o 


1NPUT>  SEPTUM  THICKNESS,  sEPT,  JN  INCHES#  AMD  LATTICE  VECTORS 
SSI  AND  3S2  IN  Inches. 


READ  (5*600)  SEPT,SS1*SS2 

INPUT>  number  of  modes  for  APERTURE  FIELD  APPROXIMATION#  AND 

limits  of  sampled  (.rating  lobe  space*  pi  and  Ql, 

read  (5*6)0)  NMUUES«Pl,(>il 
mlo«i 

IF  (8EPI.lt, l.E-10)  ML0*2 
NOaNMUOtSf 1 
NmOOEsNMOOES 

READ  (10,600)  A*S,0,B1,EPS 
IF  (EOF(10),N£,0)  CO  TO  190 
AA>2.#(A>BtD) 

write  (6,900)  AA, Bl, A, B,D, Sept, EPS*NmUDE3, SSI, 8S2, Pl.Ul 
AP»2./(AA*TPI) 

write  (6,990) 

00  ISO  IFal.lOl 
read  (10,800)  KA(IF) 

APl«KA(IF)#AP 

al«a*api 

8LaB#APl 

ol«dfapi 

BIL»B1#AP1 
M2Ba(0.5/BlL)#F2 
8EPTL*8EPT*AP1 
Sl(l)»APl#SSl (1) 

81(2)«AP1#S81(2) 

S2(l)«APi#S82(l) 

82(2)sAP1#SS2(2) 

S’^^aSl  (2) 

S..'>S2(2) 

READ  (10,600)  (T1(I)*I>1*16) 

I0«0 

DO  120  I«l*16 
TKl)»Tl(l)/KA(IF) 

M0«(I*l)/9tl 


1 ao 


Ms*^i0/3 
Do  UO  J*l»a 

I0*i0  + l 
MgMfl 

S*1  .-M2B»M^»2-T1  (lUABSCTKi)) 

GtSURTCABSCSn 
If  (8,LT,0,)  G«*G 
T(IO)«G 

110  CONTlNUt 
UO  continue 
J*1 

130  X«-1.E*25 

DO  HO  lalfbA 

IF  (X,GT.T(1)J  GO  TO  HO 

IASH(I-1)/16 

IRalHl-(lA-n»l6J/«+l 

IC»I-(IA«1)*16*(IB*1)*« 

loal 

XaTCn 

HO  CONTINUE 

GAMMA(J)ax 

T(l0)a-l,E+30 

IE«<**(IA*1HIB 

K(J)arUIE) 

GaEP8-l,+K(J)fABS(K(J)) 

8aSUPT(ABS(G)) 

IF  (G.LT.O.)  Sa-8 

KECJlaS 

M0DE«IA/3+l 

MOOtl (JialMOOE(MOUE) 

M0D£8M0U(IA»2)tl 

lSYM(j)a3YMCM0Dt) 

NMdCjaNMdC)  tl*lA/3 

NEdCjawEdC)tIA/3 

N2ad-IA/3)*NA^dCH(IA/3)*NpdCi 

lClaIC-IA/3 

NN  C J) aN2 

IF  dA.NE.3,OH,ICl.Nh,0)  GO  Tu  130 
Encode  do»9io»NiouORD(j j ) wodei  ( j)  ♦n2»  ici 

JaJtl 

IF  U,LE,NMU0E8J  GO  TU  130 
CALL  NORM  (NMODESf IwRITE) 

UaO, 

VsO, 

lSGa{3APHl)»(2*QHl)tOH-l 
Do  ITO  ILO*'^LOiMEU 
L0HH2HL0 

IF  {IlO.EG.2)  LUMlagHHl 
IF  CILO.EG.2)  31 (2)aO,5»Sl2 
IF  (IlO.EG.2)  S2(2)aO.S+S2? 

CALL  SCTmaT  (LOHl»NN-OUH'),VtSU  »S2t  ♦■'•■'•ObflSiG) 
IF  (IL0.EQ.23  GU  lu  18(' 

PT(IL0)*CABS(S21  nSG»  ntS21  ( I SG  » frO  3 J 
PT(lL03a«EAL(yAdSG)/V(n3*FT(Un) 

IF  tPTdL0),LT,l,t-10i  PTdiOJsi.E-10 
PTDb(lLO»IF)*lO,*AU)GH(PH  iL'.i)  J 
PRl  dL0)a0.b*CAb8(S11  ( t » n tsn  ( 1 .iiO) ) 

IF  (PRl  dLOJ.LT.l.E.lO)  PRl  fiLtUsi  .t-10 
PRlDBClLOjalC.^ALOGlOlPRl d| 0)) 
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PR2*0,*)*CABS(SU  (NOt  U+Sl  \ (MOt^V)  )*»2 
IF  (PH2.lt. I, t-iOj  pR2*i.flO 
PR2UBslU,»ALOGlO(PP2) 

GO  TO  IbO 

IbO  PT(UO)sRtAL(YAUSG)/VU  J)*CAtjS(S21  (ISGf  U)**2 
PTat0)*2.»Aes(PT(lL0)) 

IF  (PTULOl.LT.l.t-lOJ  PT(IlO)»1,E-10 
PTOb(lLUfIF)«iO.FALOGlO(PT(iLO)) 

PRl  ClLOJ«CAbS(SUntU)**2 
IF  (PRKILOJ.lT.I.E-Io)  PRl(ILO)*l.t«10 
PRIUB(ILO)«10,FALOGIO(PH1  (UOJ) 

IPO  CALL  COMSRV  (SU»S2l  *N'^00»IsIG»|M) 

IF  (ILO.EU.I)  WHITE  (b»950)  KA ( IF ) » LOmI i P 1 ( 1 LU ) * PTDB ( ILO  ♦ IF ) . 
IPHlOBdLO)  tPR2Dbf  IM 

IF  (ILO.FU.2)  WHITE  (o»9feO)  Ka (J F ) ♦ LOHI »PT ( iLU ) . P7DB ( ILU» IF ) t 
IPR1U0(ILO)»IM 
170  CONTlNUt 
ibo  Continue 
GO  TO  100 
190  Call  exit 
c 

eoo  Format  (bfio.o) 

810  format  (IblE) 

900  Format  (lMlf««Xf42riL>UAL  frequency  array  element  design  curves. 

!////» lOX,  13HELF.MENT  OATAA./,  15X.3HA  « » F b . 3 » bX . 3hB  s,Fo,3tbX» 
27hALPmA  a»Fo.3«bX»6MBETA  *.F6.3i5X » /WOELTA  s.Fb,3»5X» 

38HSEPTUM  s.F6,3»/. 

AlSX.bhERS  s.F5.2»bX.I2»bH  MoOES ♦ // . I 0 X » U HAHRA Y DATAA./. 
SlbX.AhSl  a.F6,3* IH» .Fb.3.bX,«H82  s»Fb.3» Ih. »FE.3»5X.«hPI  «»I3» 

bSx.RHui  J.U./.15X.29MALL  Dimensions  in  inches) 

910  Format  (a3*2I2.2) 

9i0  Format  (2X.A7.3x. Ab.2x»3F10,5»I«) 

940  FQRMAJ  UHl .AbX»38HPOAEH  TRaNSMISSTuN  FACTOR  AT  BROADSIDE./. 

159X. 12HVERSJS  k*A/2.////.10h  kfa/2  . IOhEXCITaTION. lOM  PT 
2.10H  PTCOB)  .lOh  PHllDB)  .lOh  PR2(Db)  » lOh  !M  ./» 
310X.10H  type  »/) 

9b0  format  l5x.F5,2»bx.A2.«X.«(2X.F7,3,lxJ.l7J 

9b0  format  (3X.F5.2»bX.A2.ax.3C2X.F7,3.lx)»l0X.i7j 

F nD 
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SLlBROuTlNfe  LSMLSt  (N^OntS) 

C 

C COMPUTE  CHARAClEKiSTir:  HumS  '-»■  r.TSPt'^STU'''  ►'tLATji,  . p.’ik  bY,  -'tf-’lL 

C INHOMOGENEUUSLY  inAOfcn  HECTA-.f^'l  aR  GUIuF. 

c 


c 

DEFIMTIONS& 

c 

B 

m 

GuIUE  MAPKOft  OP>'EvSinr-,  OivTDEt.  oY 

c 

Free  space  /.Avf.LFifCiiH 

c 

EPS 

m 

REIATIVE  OlELfcClRlC  CC'MSTA’.T  i>h 

c 

1 uagjmG 

c 

aL 

m 

DISTAMCE  FRQM  GuIUF  ‘-'AIL  TO  MtAKtST 

c 

Fl)Gt  OF  L'-IAHIi^G  si  Ah  fwAVtLtf-GTHSje 

c 

BL 

m 

ntSlAMCt  FRuM  SLAo  E r>Gfc  To  bYr->'tTKV 

c 

PlAwF  <»..  AvFLFr-G  I hs** 

c 

OL 

- 

SlAh  THiD'^tSS  •'•■^AVtLLtJGTnS?: 

c 

septl 

m 

SEPIUN  IhICKLlSS  aITm  RtSPEtT  Tl;  F«EF 

c 

SRAlF  a'AV/ELENGTH 

c 

TPI 

m 

6. 

c 

81 

m 

FIRST  Lattice  ''tCTow 

c 

82 

- 

SECi''  n LAI  TiCE  vEC  h R 

c 

K 

m 

X-I‘1PFCIEL  / T‘i  Alh  hFGlO''‘  OF 

c 

n-jll-F  JITH  HtSPtCI  in  FRPt  SR/-CE 

c 

v'fA  Vt^  L'*oFR 

c 

kE 

m 

X-l'|P(ClFu  V A vEkF  U'vF.h  !(.■  OIFlElIkIL 

c 

RFf.ir*.  UF  GuIUfc  i-ITh  RtSPECI  TO  FREL 

c 

SRACF  .'/.v'F'-u^bE> 

c 

gamma 

m 

LI"- GTU!'1;I  vAL  AAvEiv"!'Rfc.R  N;C'»MALl/tU  10 

c 

FmLl  ?pacE  AAvtnowhtR  ORPpwfcn  Ry 

c 

I .Ct-FASr.G  LUI*uFF  FRblOoE'N.CY 

c 

MOUEl 

• 

tyhf  ( I ,fc  , f LSfc  OR  1 S") 

c 

ISYM 

m 

SY‘"  F I Rit  (I  .E  , • S HR  A ) 

c 

NN 

m 

X ■''•I  F 1 iG  RAH  Ar‘F  TE  w 

c 

MM 

m 

Y rVf'F  RT  G RAkAi'iF.  IFR 

c 

MOOOHD 

m 

VtClPR  Cf'-f.TAliMi-G  ‘•’CDt  OFoIGGaTIOwS  -- 

c 

ilRi'FRfci/  HY  lNCPtASU.G  CuT-OFF 

c 

FiXF-vtif-CY 

c 

C,*^^$i$t^:$$if$tHi*$*^ilt$$****^^'^******ili**$^*^**1f***$****1^***^***$$i^*1f*$****** 

dimension  gam  cat  lO  1 1 1 ) f I^•fJPf  f A?1 ) ,su  (Uf  10»  11 ) ?SE  (Of  1U»  1 n »Nr^(  1 n » 

INECIl} 

REAL  KAP♦KAPEfM28*K,KE♦KT^C 

COMMON  /array/  ALiBLft'LfflfSfPILf  Tpl  fPHSfSl  (t>)  fS2(2) 
common  /modes/  K(aO)  iKE(20)  a(?0)  fMOOEl  (2u)  f ISYM(20)  f^.M(20) . 

IMM(20) iMOOORO(20) 

Al«TPl»AL 
Bl«TPl*BL 
Ol«TPl*DL 
eP8l«EP8-l, 
tP80»8QRT(EPSl) 

M2B«(0.5/B)P«2 
C 

C COMPUTE  K FOR  M»0 

c 

DO  ISO  MODEsUa 
MdBlHS 

IF  ((M00E/2)*2.E0.Mo0E)  MoejHA 
mi»}Ml8m 

IF  (M00E,CT,2)  M1.3HLSE 
M2«i0 

IF  (M00E.GT.2)  M2«ti 


186 


M3»M0Ufe/J 
KAP»*EP&Q^1 »E»10 
MO»i 
McM0-M3 
DO  160  Nslill 
J«l/^ 

Sst.Ol 

IF  (KaP.LT.U.)  8*0,99 
KaP*J*KAP+( l-J)*S»kAP 

KINC*0.319  a 

kap*kap»kinc  I 

1*0  3 

KaP*KaP*KINC  ? 

DIFF«DIFFI  1 

I«Ifl  1 

KAPfe«tPSlfKAP*AtfS(KflP) 

SsSUPI (ABS(KAPEJ)  j 

IF  (KaPE.LT.O.)  S»-8 

KAP£*8 

OlFM*niSP(MOOE»KAP,t<MPfc  »a1  »R1  »U  f EPS)  I 

IF  (I.EGi.l)  GO  TO  120  i 

IF  {KINC.LT.I ,E-lu)  GU  TO  1 30  i 

IF  (Ad8(DIFF*OlFFl),UE.l.E»lO)  GO  Tu  130 
IF  (01FF»01FFI)  110»130fl00 
KaP*»'AP«KINC 
KlNC*0,5»KlraC 
OlFFlsOIFF 
60  10  ion 
0IFF«DIPFI 

60  to  100 

IF  tAbStOlFF  ),Ll.AB8(0IPFn)  GO  TO  l«0 
60  TO  ISO 
KAP*KAP»KlNt 
KAPfc*EPSl tKAP^AttSCKAP) 

SaSyHTCABSTKAPF.)) 

IF  (KaPE.LT.O.)  S»-S 

kape«s 

80(MODE»NfMO)«KAP 
8E (mOuE ♦N»Mu)«KAP£ 

GbI ,-M?B*M«AH8(kAP)»KAP 
SaSUPTCABSCG)) 

IF  IG.LT.O,)  8*-8 
GA^IC^OOE  (Nf  MO)bS 
CONT INUE 

continue 

INCLUDE  T CI,E,»  dependence 

do  183  MnoE«l t« 

M3BM0DE/3 
Mgsi n+Mi 
DO  182  NBlflO 
KAP»S0(MODEfN»l) 

KAPE*SECMUDEfN«l) 

DO  I8l  mob2»m2 
80(MOdE»N»M0)*KAP 
8E(MODE»N»mo)*KAPE 
fis(mo«m3)F«2 

Bml*W2B*m 

G2»-KAP*AB8CKAP)+l,0-bML 


1 R7 


GAVCN'ODfcfNf  v()  jsS^ikT  ( At  S(f5(?) ) 

If  (S2*UT*0,)  GA'-^CfOOt  '<UDt»NfNO  J 

cijf^riNUE 

COk I iNUt 

COMIiSlUfe 

WHIIE  (6,«<0U) 

nn  iks 

A<  ( 1 ) 8 

ME (1)80 

COM! lUUt 

Isl 

GA^’i'iAd  js-1  ,fe  +30 
DO  iidu  ‘'"DoFsltU 

►"ASlhs 

IF  { Ci-’DDE/2J  ) M«4S,ha 

Vis3Mi,Sf^ 

IF  (f'0nfc,GT,2)  MlsJHLSfc 

M28in 

IF  (^!UDE,GT,2)  M2ett 

M38Iv,OoE/3 

DO  21u  M08l»f^2 

DO  20u  Nsl«10 

IF  (GAMMA{13,GT,GAM(Muot«?.fNon  r.u  In  2yo 

GA'-^mA  ( I 3sGAm(.v.0UE  »w,yo  ) 

Encode  Uo»viu«lMur>En  ))  ■'•i ♦''a*' tf 

I ASivtOuE 
IHSN 

IC*MO 
COf"!  iMiit 
COM) lUUE 
COMfiNUE 

QROEk  modes 

GAMUAf  IBf  K)8-1.L  + 3« 

KdJssriciflf  IR,  TO 
KF.d)sSECIA,lBfiO 
siHLSw 

IF  dA,GT,2j  •UsimsE 
IF  (lA.UT.S)  lir'UIOsMr'dO  + i 
IF  dA,GT,2)  r;F(R3sMe  dC)  + i 
N2siMM(It  ) 

IF  dA,GT,23  iJ2S'\E(]C3 

lCUlC-lA/3 

MODE  I (T3sM 

NMd)*M2 

MM(I3sICI 

ISYMC 1 )sih8 

If  dlA/23*2,ED,  lA)  TSv^f  1 )sl»-A 
EhCUDE  ( 1 ('tVoOiMOUOkOd  3 ) >-'i»'v2fTC) 

(ft  *020  3 Ir-.DUE  ( 1 3 *'''ni'ripnd  ) ,c,Ai-.MAC  ] ) , I 

Tal  + l 

IF  d ,LE.M*’'00ES3  00  ru  IV(. 

SEC*SECOiit  ( <3 
write  (ft*930)  SEC 
reIohn 


F0«MAI  (1H13 

Format  (a3»ai *212.23 


188 


function  DI9P(MfKfKEfA«BtO«C>*) 

COMPUTE  DISPERSION  RELATION  0(K,kE)  FOR  ARBITRARY  K ANQ  kE 


OEFINITIONS> 

M 

K 


KE 


A 

B 

0 

ER 


MODE  FUNCTION  DESIGNATION  - SEE  CODE 
X-DIRECTEO  wavenumber  IN  AIR  REGION 
WRT  KO 

x-oirected  wavenumber  in  dielectric 

REGION  WRT  KO 

alpha/lamboao 

BETA/lAmBDAO 

OELTA/LAMBDAO 

RELATIVE  PERMITTIVITY  OF  SLABS 


real  KfKE 
9A>SINC(KAA) 

CA>COBC(KPA) 

SBBSINCCKPB) 

CB«COSC(K«B) 

SD«82NC(KEP0} 

C0«C08C(KEP0) 

SkA>SXOX(K*A) AA 

SkB«8X0X(K*B)T8 

8kE»SX0X(KE*0)FU 

81«S1GN(U«K) 

S2>3lGN(I,tKE} 

GO  TO  (lOUf  U0(l20f  I30)t  M 


C 

C l8m  S^^MMETRIC  mooes 
C 

100  OISF«ERPK«SA*S1P(CB«CO*ERaK«SB*S1«SKE) 

0l8P«0lSPTCA*(ER*K»8B*8;*CD+KE*8r>»32*CB) 

return 

c 

C l3M  aNTI-SYMMEIRIC  modes 

c 

no  DI8P»ER»SA*SI*(SB+C0+ER*K»SkE»CB)4CA*(KEP5D»82*SKB-ER*CL)*C8) 

Return 

c 

C LSE  symmetric  modes 

c 

120  Ol8P«SKA»(KE*80*82*CBtK»S8»Sl*CD)tCA*tK*SB*81*SKE»C8*CD) 

return 

c 

C L8E  ANTI»8YMMETRIC  MODES 

c 

lio  OlSPiSKAiMCdPCD-KEPSKBTSOASjTtCAACSKBACO^SKETCR) 


SuBKOuTIN'E  imOHM  (NMuDtSt  U«iTE) 

8UBR0U1INE  NURi^i  COMPUTES  ThE  FtFO^iUlUE  ■''UDt 
values  for  The  COEFFITIFNiTs  APPFARINI,  liM  THE 
AMO  IHAT  •FOR  WHICH  SEE  OOcUMFMA T lOM* . 

NOlEi  THE  coefficients  b?P*ElR«F2R»FP»BiOH* 
REAL  OK  imaginary  #fvOT  COKpLEX^.  ALL  UThEH 
ALaA>S  real,  ALSOi  riOKMALjZAT JOMS  CARRY  AN 

oehenoence  which  Cancels  in  scattering  hatri 
MOT  IN  MUOE  FUNCTION  COMPutaTUIOnS, 


DEFINITICMS>  see  SUpRuUTJmE  LS^LSfc 


NOHMAUZATimiSf  AMU 
FXPRES8I0NS  f-UR  VHAl 


CDP»E1UP»E2DP  hay  be 
coefficients  are 

AODniONAL  1/l.AMBOA 

X compuTations»  HUI 


PRIN'E  (L8M  MODES) 

OOObLE  PRT^E  (LSE  MootS) 
SYMf.FT'UC 
anti-symmetric 


C 

C 

c 

c 

c 

c 


l?)*SS2(2) 

lO),CPClO)tCDPUO)» 

iO)»Nps(io)tNPACio)» 


SUFF1XES> 

P 

DP 

S 
A 

REAL  KAP,KAPEfK»KEfNP5iNPA»M0PS,NDPA 
Common  /array/  ALfRLiDLfRlL,SEPTl.,TPl»EPS»SSl 
common  /COEFS/  BlP(lO)fB2P(lO)tblDPC10)»B20P( 
lElPUO)iF2P(10)»EiOP(  10)  »E2oP(lL)  ,FP(10)  »pnp( 

2N0P8UUJ  ♦NOPA(lO)  ♦K0UNTC20) 
common  /MOOES/  RAP(20)iKAPE(20) tRAMlpO) fMuOEl (20) ♦ISYmI (20) » 
INn(20j ♦MM(20) tMUD0RD(20) 

A«TPI»AL 

B«TPI*BL 

0«TPI$OC 

IF  (IwRITE.EQ.n  WRITE  CbtRoOj 

IL8m8»0 

IlSMA.O 

IlSESsO 

IlSEAsU 

DO  IRO  lalfNMOOES 
K«KAP(I) 

KE«KAPE  (n 
akbkpa 

BK«R*B 

dke«d»ke 

MOOE«MOoei(I) 

I8YM*ISYM1(1) 

M«MM(I) 

IF  (M0DE,EQ,3hL3E) 


GO  TO  Ho 

lsm  modes 

IF  (ISYM.tO.lHA)  GO  TO  100 

symmetric 


iLSMSslLSMSfl 

K0UnT(I)»IL8MS 

8lPaLSM8)«C08ClBK) 

B2PCILSMS)«EPS*R»SInC(BK)/Ke 
B2P(IlSMS)«B2P(1LSMs)*S16M(j ,»K)fsi6M(l,fKE) 
CpaL8MS)»-SlNCCBR)*(C08C(0KE)+Cn8C(8K)f8INC(0K£)/B2PaLSMS)) 
CP(IL8M8)«CP(ILSMS)/SImc(AK) 
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ANOKMsO,5»8R*(bL<'Sl  (<i,»BK)4fDL*Sl  Cd  .♦DkE  ) *b  1 F C U SMS)  ♦♦a  + l)L<‘Sa  (2,  » 
lDKE)*Sn»Ntl.»KE)*H2P(lLSMS)**2)/FpS 

2»2.0*fcaPULSMS)»B2P(Iu8MS)*nL*S3(DKt)»SlGlv(l.iKe  )/FP8 
S+AL+SI  C2^*AK)$CP{11  SM{j)»»2) 

NpS(lLSM8)»S0RT(AbS(AiM0KM) ) 

GO  fO  liO 

c 

C A^TUSYMMEIRIC 

C 

iOO  iLSMA.iLSMAf I 
K0UNT(I)*IL8MA 
BL»BL*S1RN{i,,KJ 
ElPdLSMAjsbTi'lCCBK) 

E2P(lLSMA)«-FPS»K*COSC(bK)/KF 

FP(ILSMA'«CUSC(BK)*(COSC(DKE)-Kf  ♦tai'  TiBn)^Sli'JC(l;Kt)»Slf;'^(  1 , tKE)/ 

1 (EPS*K)  )»SIG*l(U«K)/SRiC(AK) 

FaEjP(lLSMA)^fc2P(lUSMA) 

AN0HM«0,S»BlL»(bL*S2(^,*6K  )+OL»l81  (2,<'DKt)*h  1 P U LS'*A ) **2  + 

I (82  C2,*0Kt)  ♦EaPULSMA  )♦>>?•?, *pFS3rnKE) ) *81  Givti  .♦Kfc  ) )/tPS  + 

2Al*S1 t^•«AK)♦FP(ILS^  A)A»?j 
NpAnLSMA)*SrjRT(ABS(  AijnpM) ) 

BL*BL*SlGfN(l.fKj 
GO  10  130 
C 

c lse  modes 

c 

UO  Rm«1*0 

IF  (M.EO.O)  RMS2.0 
IF  (ISYM.tO.lHA)  GO  TO  120 
C 

C symmetric 

c 

iLStSsILSES+l 

KOUnT(I)*iuSES 

8lDP(ILSES)sCUSC(HK) 

B2DP(IL5ES)«K*81MC  (BK)»SIGM(1  .»'')*8IGN(1  . *Kfc)/Kfc 

CDP(ILSES)*C08C(B^)♦(L^SC(0K^•K♦TA^(T(BK)♦Sl^.'C(0Kt)*SIGN(  1 . fK  j/kE) 
l^Sl6'^n.»t')/SINC(AK) 

AnOKMsO,5'I>«M<‘B1L»(PU*S1  (2 , *bk  ) tUL*  ( S J ( 2 . FDKt ) fh  1 [)p  ( RSES ) A*2  + 
18lGNCl.iME)F{S2(2.*0«t)*82np(lLSFS)F»2-2.»BlDP(U  sf S)* 
2B20P(lLSES)*83(DKfc)) )+AL*82{2.Fak)*SIGNM ,»k)ACDP( TLSfcS)»»2J 
MDPS(IL8E5)«SUH1 (ARS(ANURM)) 

GO  TO  130 
C 

C ANII-SYMMETRIC 

c 

120  IcSEAsILSEAf I 

AL»AL*8IGfY(l,»K) 

BL>BIfS1GN(1.»K) 

K0UNT(I)*ILSEA 

ElOP(IL8EA)«-SIivC(BK) 

E2DP(ILSEA)*KAC0SC(BK)/kE 

F0P(lLSEA)«-8INC(bK)*(C0SC  (D»<E)+K*SlNC(DKt)/'  (i'fc<‘TANT(BK) ) )/ST.-jC(Ak 
1) 

F*EIDP(1LSEA)FE2DP(ILSEA) 

ANOHMaO,5*RM*BlL*(8LFS2(2.»RK)tOl  *(Sl (2,*DKE)»FlDP(ILSEA)#*2+ 
tSiGNd .♦KE)*((82(2.F0KE)AE2DPCTLStA)*02-2,0»F*S3(DKE) ) ) ) + 
2AL*32{2*faK)»FDP(1LSEA)fa2) 

N0PA(ILSEA)«SUR1 (ABS(AMORmI) 

AL«AL*S1GN( 1 , »K) 
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c 

c 

c 

c 

c 

c 

c 


SuBROuTINt  SCTMAT  (LOhI  »NyODEStUO»VOiSU  »S2l  fiMMOOfSIGn 

SCI^AT  computes  the.  fFtOGUlDE  • FREE  SPACE  SCATTERING  MATRIX 
blocks  SU  and  821  FOR  An  iNFIMTE  RECTANGULAR  GRTD  ARRAY  OF 
T^IN  DIELCTRIC  slab  loaded  RECTANGULAR  YjAVEGUIOES. 

lattice  vectors  si  and  S2 

LOHI  SPECIHES  a'hICH  FREQUENCY  RaND  ►L0HI«2HLn  FOR  LOW  FREQUENCY 


c 

r 

band* 

L0hIs2mH1  for 

MiGh 

Frequency  «anup. 

w 

c 

notes 

IF  the  septum 

thickness  is  not  equal  to  The  wall  thickness 

c 

L0hI»2MLUt  and  TmE  high 

and  lo-a  frequency  unit  cells  are 

c 

C 

IDENTICAL. 

c 

c 

DEF1nITION3> 

c 

LOHI 

M 

«2hlo,  for  trianguiar  grid  or  Thick 

c 

SEPTUMF  a2HHI»  FOR  THIN  SEPTUM*  HJCH 

c 

FREQUENCY  BAND*  AnO  RECTANGULAR  GRID 

c 

nmuoes 

m 

numbfr  of  fefoguide  modes  used  to 

c 

APPROXIMATE  APERTURE  FIELD 

c 

uo 

• 

81N(THETA)FCDS(PHn 

c 

VO 

• 

S1N(ThETA)*SIN(PH1) 

c 

Sll 

m 

FEEDGUIOF  self  REFLECTION  SCATTERING 

c 

BLOCK 

c 

821 

• 

FEEDGUlOE  To  SPACE  mOuE  VOLTAgE 

c 

TRANSMISSION  coefficient 

c 

NMUO 

• 

Numbfr  of  feeoguiof  mopes  in  unit  cell 

c 

SIGI 

m 

number  of  space  modes  • 2FPIF01 

c 

Y 

m 

FEFOGUlOE  MODE  ADMITTANCE 

c 

YA 

m 

SPACE  mode  admittance 

c 

***  FOR 

other  definitions  see  sub  lsmlse  *** 

C 

C SUFFiXES>  SE^  SUBkOuTINE  NORM 

DIMENSION  AM(20J 

REAL  KAPfKAPE»K»Kt»MPS»'vPA*NOPStMDPAfKZl  »KTC250)  ♦KX(250)  »kY(250) 
complex  KZ♦AJ»CUE^M(^bO)  fF.SN(25u*20)  ,SU  (20»20)  lYf  YAfEXPBi  (250)  f 
1S21 (250»20)»FXPb(250)»YRIP(pOt2U^ ,IWTGRL 

Integer  p»RtPi *ui »sir,tsiGi isige 

common  /array/  ALf8LfUL»HIL,SEP1l.  ,TPItEPS*3l(2),S2(2) 
common  /CN3HV/  Y(20)  # YA(250) 

common  /C0EF3/  blPf lOi»fa2PCtO)tBlDP(10)»82DP(lO)fCP(lU)tCDP(lO)» 
lElP(lO)  ♦F2PO  0)  fEJDPt  IP)  fE2rP(tO),FPf  iO)  tFOPClO)  ♦NPSUO)  fNPA(lO) , 
2NdPS(  10)  tNOPAnO)  »K0"NT(2o) 

common  /COSSIN/  8kAPAL(20)»CKAPAL(20)»8kAP8(20)»CKAP8(20)» 
1SKAPD(20)»CKAPD120) ♦SkxaL(250)»CKXAL(250) fSKXA(250) »CKXA(250) f 
2SkXB(250) tCKX8(250),SKXD(25o)»CHXD(25U) 
common  /MODES/  KAP(20)  fKAPE(20)  fr,AM(2U)  fM0DE(2«)  »ISYM(20)  »NN(20)  ♦ 
IMMC20) tMODORDCPO) 

Common  /pq/  p»Q 
Data  aj/(0,»i.)/ 

100  AbTPUAL 
8*TPIABL 
ObTPIAOL 
BB«IRlFblL 

sep*tpifsepjl 

FPSU«1./TPJAF2 

Pl*2FP4-l 

Qls2AQfl 
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o r»  o o r»  o 


8ICI«2*P1*01 

8lG2«3IGl/2 

NmOO«NMODL8 

IF  (L0HI.EQ.2Ht0)  NmOO«2*NMoOE8 

DO  UO  NilfNMOOES 

KE«KAPE(N)*U 

MsKAP(N)4:A 

XbKAP(N)»B 

AM(N)»O.5*HH(N)/01L 

T«GAM(N)/(l,-KAP(N)*Ae8(KAP(N))) 

Y(N)«CMPlX(T»0.) 

IF  (6AM(N).LT,0.)  Y(N)«CMPLx(0.»T) 

IF  IH0DE(N),£0,JHLSE)  Y (N)«i ,/Y (N) 

8kAPAl(N)«3INC(W) 

CkAPAlCN)«C08C(w) 

8KAp8{Nj«8lt\JC(X) 

CKAP0CN)«COSC(X) 

8KAP0(i\)s8INC(KE) 

CKAPO(N)*COSC(KE) 

IF  cLOHi.EO.gHHi)  CO  10  no 

Y(NYNMODES)«V(fO 

no  comhnue 

compute  inverse  lattice 

T*8l(n»S2(2)-82(l)Y8l(2) 

CELLA»1,/SQRT(A0S(T)) 

T«l./T 
TlX«T»S2(2) 

T2X»»TF81 (2) 

TlY«»T*82(n 
T2Y»T$8l(l) 

81G»0 

compute  free  space  <^fAVE  nUvRERS  A^D  nAvE  AUMITTAnCES 

no  130  L«i»2 
Ll»2-L 
RL1»L1 
Ll«L-l 
RL2*Ll 

00  130  Jl*ifPl 
J»Jl»P«l 
UoUU+J*T 1 X 
vsvof j*Ti  y 
DO  l3o  KltlfOl 
8IG*SIG+1 
KbK1«u*1 
Kx(SlG)sUtK»T2X 
KyC8IG)sv+K$T2Y 
K21s1,«kx(SIG)Y»2-ky(SIG)**p 
AA*SQRT  (ABSCKZn  ) 

KZ*CWPLX(AAf 0.) 

IF  (KZ1.L7.0.)  KZsC«PLX(0.,,AA) 
KT(SlG)*8QRT(A8S(l,-KZn) 

IF  (Aa.lt, l.E-lU)  Go  TO  120 
YA(SIG)»PL1/XZ+HL2*kZ 
GO  ro  130 
120  YACSIUlPPLa-XKZ 
130  cominue 
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compute  sines  ANU  cosines  of  KXf(EUEMENT  UIMbNSTUMS) 

DO  Hu  SIGBltSlOl 

TaKX(SlG)«A 

8kXAUS1G)«SIN(1} 

CKXAL(8IG)bC0S(T) 

TbKX(SIG)*U 

8KX0tSlG)*81NCT) 

CkXU(SIG)«C08(T) 

TsKX(SIG)*B 

8KXb(SIG)*81N(T) 

CKXb(SlG)sCUS(n 

TBKX(SlG)*(AtB>U) 

SKXA(sIG)bS1N(T) 

CKXA(S16}>CUS(T) 

TBKY(SIG)«eb 

EXPbl  (81R)*CEXP(Aj*n 

TbKY{SI6)»(B9+SEP) 

EXPb(8lU)«CEXP(Aj»T) 

continue 

COMPUTE  coupling  coefficients#  ESN(SICfN) 

DO  18o  NBlfNMQOES 

TbAm(n)*B2 

T2bAM(N) 

If  (Mm(N),EU.O)  T2B1.U 
DO  HO  SI6sI«SIGl 
TlBAbSO'YCSlG)} 

IF  lA8S(T-n»*2).LT,UE*lO)  60  TO  ISO 
C0EFMC3lG)»n ,«EXP81 (SIG) *)*PMM(NJ )YT2/CT*TiB$2) 

IF  (MUOE(N) ,EU,3ML8E)  COEFM(SI6)bCOEFM(SIG)FKY TSIG)/T2 
GO  TO  1S5 

CO£FM(S1G)bo,5*AJa8B 

IF  C^M(N),EU,0)  CUEFM(8IG)Bp,BC0FFM(Sl6) 

IF  (MOOE{N).EU.3HLSM)  COEF  M(SI6)bcOEFM  (SIG)  »3IGf,  ( 1 , * K y (STg)  ) 
CoEFM(SIG)»COEFm(8IGI •CELIA 

continue 

DO  170  SIGBlfSIGl 
LR«1 

IF  (SlG.6T.SIG2)  LPbS 

ESNCSlG#N)*INT6HL(Kx(SIG) »Ky(SIG) tKTCSiG) »LHiSlG#N) 

E8N'(SIGiN)bfps0»C0EFM(S16)*eSN(SI6#N) 

IF  aOHl.EOtZHHl)  GO  TO  170 

E8N(SIG»N+nmOOES)«CONJG(FxPr(SIG) )*ESN{SIGf N) 

continue 

continue 

form  scattering  matrix  Blocks  su  anu  sai 


form  matrix  triple  product 

MAT(TRIP)*mAT(C0nJ6(E8n) )*MAT ( YA)YMATCE3n) 

and  matrix  FOR  inversion. 

DIAG(Y)tMAT(TRlP) 

DO  210  lABl.NMOU 
DO  200  XBBlfNMOD 
TRIPUA,IB)b(0.#O.) 


SlHUtlBJKO.tO.} 

IF  (lA.tQ.IB)  Sn(lA»lR)«2.*YtlA) 

DO  19u  SIOltSiGl 

TRIPCiAf  TB)«THIP(IA»IB)^C0NjG(ESN(5lG»lA))<'tA(8lGj*t8M(Siri»lB) 

190  CONTlNlJt 

200  continue 

210  CUNTiNUb 

DO  220  lAsltNMOD  * 

TRlP(IAtIA)«Y(IA)+TRIP(IA,lA) 

220  continue 
C 

C CSIMEQ  returns  MAT{SU+DEL(I»J)1  ‘VMERE  OELCltJ)  IS  ThE  KrOnECKEH 

C delta  FUNCTION 

C 

Call  csimeo  (TRiP»NMOo»sii»NMun,KS) 

c 

C SOLVE 

C MAT(32l)«MAT(ESN)*'<AT(311fDEL(ltJ)) 

C 

DO  2SU  SIGaifSiUl 

Do  29u  IBslfNMOD 

d21(SlGtl6)s(0«f0,) 

DO  230  lAiltMMOU 

82l(8lGf  IB)»82USIG»Ib)^ESN(8lG.TA)*SlUlA»iB) 

230  CONIInUE 
2ft0  CONTINUE 

2b0  continue 
C 

C SOLVE 

C MAT(S11)bMAT(31U0EL(I»J)) 

C 

DO  260  XAslfNMOO 
Sll(lA*lA)BSn(lAflA)«i,0 
260  CONTINUE 

return 

END 

#6 


o o n o o o o o o r>  o o r>  r>  o o o o o o o o 


complex  FUNCTlOiv  IMTCiHL  (^  X ,KY  t k T,LK i SIG  t Im) 

function  intgrl  computes  The  t^tegmal  (n«  x)  portion  he  the 
COUPLING  coefficients*  ESN(S16*M) 

n01E>  mnemonics  AHE  CHOSFv  to  cOiNClOE  <AlTH  MOTaTIO.,'  in  report 


defimtions> 

KX 

KY 

kt 

LR 

SIC 

N 

♦♦♦  FOR 


- X*niREC7ED  wAVENUMBF R OF  lOhF  «RT  ko 

• Y-PIPECTFU  aAVENOMHeR  of  LH6F  *^RT  KO 

• transverse  wavenumhem  of  uObE  rmt  KO 

• s}«  FOR  E»M00ES»  s?»  fur  M-MOUES 

m Nui^'bF  R OF  GNATII'iC}  LURE  IN  InTERi'^AL 
ORhEPlNG 

• NuARFR  UF  aperture  MOoE  In  INTERNAL 
OROtPlMG 

UThfR  pEFlNTTlUNS  SFt  SUB  LS^■LSE 


SUFFiXES>  SEE  SUflROtiTINt  NORM 


complex  Rl «R?fR3*R«fRS*R6*R7*KR»PP*Kio,AJtZl ♦^?*Z3*ZR»Z5*Zb*TT*TU* 
1TV*TU1*TV1»11P»12P 

real  KX*KY  f KT tK*Kt*NPS«NPA*MDPS*NDPA,KAP»KARE 

Integer  sig 

common  /array/  AL*BL*UL*BU fSEPTL*TPT*EPS»Sl(2).Saf2) 

common  /clefs/  blP(  lOj  *fa2P(  1 ft)  *H1DPU0)  fBeOPT  lOj  fCPllO)  »CDPC10)  ♦ 

lElPUO)»E2P(10)*nDP(U))»E2nPllO),FP(lO)*FDR(  lO)tNPS(10).NPA(l0)» 

2N0PS(10)  iNDPA(U')  *KOUi'iT(20) 

COMMON  /COSSIN/  8HApAL(2ft) »cRAPm| (2u) »3nAPrt(20) *CkAPB(20) * 
1SkAPO(20) fC^APO(2o) *SkXAL (?50) *CR>4C(2bO) *S^XA C250) *CkXA(25v) * 
2SKXB(2*jO)  ♦CRXB(2'jU)  *8KXU(?Sft)  »Cr  X0(2*50) 
common  /MOOES/  RAP(20)  *kAPF(?uT  *r,AMC?u)  fMuOtl  (20) » TSYmI  (20)  ♦ 
1NN(20)  »MM(20  ) ,M(JDORO(20) 

EQUIVai  ENCF  (k1  tAI)  » (R2*A2),  (r3.A3)*  (Ra,Aa)  * (kS,AS)  t (RF»*Ab)* 

1 (R7*A7)  , (W8,  A3) , (K<),AV)  , (RlofAlu)  ♦ (Z1  »B1) » (Z2*B?)  ♦ (Z3*B3)  ♦ (Z'lfBR)  » 
2(Z5fB5) ♦ (Z6*B6) 

DATA  aJ/( 0. » 1 . )/ 

DATA  SOR2/1  .«IR213'5623730R/ 


initialize  temporary  STOrAgF 


R1*(o, 

*0.) 

R2«(0. 

.0,) 

R3*(0. 

* 0 • ) 

R«»10, 

♦ 0.) 

R5*(0, 

*0.) 

P6*(0, 

♦ 0.) 

R7»(0. 

*0.) 

R8»(0, 

*0.) 

RR»(0, 

♦ 0 , ) 

RIOb(o 

.♦().) 

Z1»(0, 

♦ 0.) 

Z2K0. 

♦ 0 , ) 

Z3»(0. 

f 0.) 

Z««(0, 

*0.) 

Z5*(0, 

♦ 0,) 

Z6«(0, 

*0.) 

MbKU'JnT  (N) 

KaKAP(M) 

KE«KAPE (N) 
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MODt«MOUFl  Cn)  i 

ISYM*1SYMI  (N)  i 

ALH»Urt  J 

alR«i,5-al«  i 

ikt»o  I 

IF  (KT.UT,l.F*lO)  iKTsl  1 

C 

C CfJMPuTfc  TERf^S  COMMON  TO  ALl  INTFQRALS  ^ 

c 1 

C KAP  imaginary  i 

c I 

IF  CK.QT.O.J  GO  TO  100  ^ 

Tu*Aj»K/(U«K»AbS(K))  I 

Tul*Aj  ^ 

RIsISkXB(SIG)*CKAPB(N j*AJ*SKAP8tM)»CKX8(8l6) )/(KX+AJ*K)  s 

R2a(SKXB(SlG)fCKAPB(W)-Aj*SKAPB(N)»CKXB(8lG) )/(KX-AJ*K)  ' 

R7«(CKXALC8IG)»tKAPAL(N)+Aj»SKXAL(8IG)*SKAPAUN))/(KX-AJ*K)  ; 

RB«(CKXALl8IGJ*CKAPAL(N).Aj*SKXAL(8lG)*8KAPAL(N))/(KXtAJ*K) 
R9»(8kXAL(SIG)*CKAPaL(N)«AJ*8KAPAl(N)*CkXAUSIG)  )/(KX-AJFK)  < 

R10«(SKXAL(SIG)»CKAPAL(N)+Aj*SKAPaL(N)*CKXAL(8IG) )/(KX+AJ*K)  J 

C I 

c kape  terms  1 

c I 

100  TaABS(KX)  \ 

TV*nE/(EPS-kE**2)  I 

Tvl*(1.0fO.)  I 

IF  (AB8(T-Kt),Gl .l.E-151  GO  TU  llO  t 

R3«U,Et?7»0.)  J 

IF  (KX.LT.O*)  A3a(CKX0(SlG)*CKAPD{N)+SKXD(8IG)*8KAPD(N))/(KX-KEJ  J 

R«all,E+27*0,)  ] 

IF  CKX.GT.O.)  AaaiCKX0(SlG)fCKAPr)(N)-8KXD(SlG)*SKAPU(N))/(KXtKt)  j 

A5aIPl*UL  \ 

IF  (KX.LT.O.)  A5alSKXU(SlG)aCKAPn(N)-CKXD(SlG)*SKAPD(N)  )/(KX»Kf.)  | 

AfeBlPl*DL  1 

IF  (KX.GT.O.)  A6sCSkXD(SIG)*CkAPD(N)^CKXD(8IG)«8KAPI)(N)  J/(KXtKEJ  i 

GO  10  120  i 

no  A3a(CKXU(SIG)*CKAP0(N)  + SKxDfSJG)*8KAPU(N)  )/(KX-kE) 
Aa*(CKXU(SlG)ACKAPDlN)-SKxO(SIG)*SKAP()(N)  )/(KX  + KEj 
A58(8KXD(ST6)>t‘CKAP0(N)-CKX0(Si6)*SKAP0(N)  J/(KX-KE) 

A6a(8KX0(8lG) ACKAP0(N)+CKxnfSIG)*8KAPD(N))/(KX+KE)  1 

c : 

c kap  weal  ^ 

c ; 

UO  IF  (K.LT.O.)  go  TO  lau  ? 

Tusk/(1,-K»»2) 

Tu1=(1.0.0.)  i 

IF  lAbS(T-K),GT.l.E-18)  GO  TO  130 

AialPl^BL 

IF  (Kx.GT.O.)  Ala(SKXb(SIG)*CKAPp(N)+3KAPb(N)»CKXB(SlG))/(KX+K)  ^ 

A2alPl$bL 

IF  (KX.GT.O.)  A2a(St<Xb(SlG)*CKAPR(NJ-8KAPbCN)ACKXb(8lG))/(KX*K)  I 

R7S(l.E+27fO.)  ^ 

IF  (KX.LT.O.)  A7a(CKXAL(SIG)*CKAPAL(N)4SKXAL(SlG)*SKAPAL(N))/(KX-K  J 

I ) i 

R8a(l.fc+27.U.)  I 

IF  (KX.GT.O.)  A8t(CKXAL(SIG)ACKAPAU(N)-SKXAL(SlG)aSKApAL(N))/(KXfK 

1)  I 

A9*IPIAAL  1 

IF  (KX.LT.O.)  A9a(SKXAL(SlG)fCKAPAL(N)»CKXAL(SlG)*SKAPAL(N))/(KX"K  | 


AlO«TPI*AL 

IF  (KX.6T.0*)  A10a(3KXAL(8lG)»CKApALCN)<^CKXAU8I6)*SKAPAUN)}/(KXt 
iK) 

&0  TO  140 

Al«(8KXB(816)PCKAP8(N)^8KAPB(N)*CKXB(SIG))/(KXfK) 

A2S(SKX8(SXG)»CKAP6(M]-8KAPB(N)fCKX8(SI6))/(KX»K) 

A7«(CkXAL(8IC)*CKAPaI-(N)+8KxAL(8IC)*8KAPAL{N))/(KX-K) 

A8»(CKXAL(SIG)*CKAPAUM)-8KxAL(8lC)PSKAPAt(N))/(KXtK) 

A9«(8KXAt.(3IG)*CKAPAL(N)-CKxAL(SlG)*8KAPAL(N))/(KX-K) 

A10b(SKXAI(S][6)PCKAPAL(N)4CkXAUSIG}P8KAPAL(N))/(kX4K} 

All«2,*KX/(KXF*2»KEPAB8(KE)) 

Al2»2,*KE/(KXP*2-KE*AB8(Kfc)) 

Ai3«2,*KX/(KXP*2-K*ABS(K)) 

A14s2,«K/(KX>I'F2*K»AB8(K}  ) 

IF  (^*00E,E0,3ML8E)  CO  TO  32o 
IF  (ISYM.eq.IhA)  go  to  230 

t8M  • symmetric  and  LSe  • antisymmetric 

IF  CA3.GT.l,Ef2S.0P,A4,GT,l,E425)  CO  TO  150 
Zl«(R3tH4-AlUP8KXS(S16} 

Z2«(Tv1»A12-R3tR4)*cKXB(SI6) 

Z3»(R3-R4-TVl*Al2)Y8KXfl(SlG) 

Z4*(All«R3*H4)*CKXR(81G) 

GO  TO  180 

IF  (A3*A4,GI,1,ET50)  GO  TO  l7() 

IF  (Aa.GT.l.EtZSl  GO  TO  IbO 

Bl«*SKAPO{N)*72/Kfe 

B2«bl 

B1«BU8KXB(SIG) 

82aB2«CKx6(SIG) 

R3«-8i 

R4««B2 

GO  TO  180 

B1«SKaPD(N)**2/KE 

B2*«B1 

B1«B1«SKXB(SIG} 

B2*B2«CKXb(SlG} 

83»B1 
B4«B2 
GO  TO  180 
Blau. 

B2au. 

B3*0. 

BtiiO. 

IF  (A7.GT,l,ET25.0R.A8.GT,l,Et25)  GO  TO  190 

ZSa(Al3-R7»K8}*SKxA(3lGJ 

Z6«(TUlFAlfl-R7TR8)*8KXA(SlG) 

60  TO  220 

IF  lA7*A8.GT.UE  + 50)  GO  Tu  ?10 
IF  CA8.GT.1 ,E+25)  60  TO  200 

bsaskxaloigi^s^apalcniyaz/k 

B6«»B5 

GO  TO  220 

B5«-8kXal(SIG)A3KAPaL(n)A*2/K 

BbBBS 

60  TO  220 

B5«0. 

Bb*Q. 

IF  (MOOE.E0.3MLSE)  GO  TO  34o 


TT»n.»0.) 

IlPtRlfR2+(blP(M)*(21f(K5fPfe)*CKX&(SIG)J 
UTT*B2P(Mj*(Z2fCR5-R6)*SKxfl(SIRj  n/tPS 

2*CPt^")*tZ5-KR9tRlO)*CKXA(Sl6)) 

l2P«Tu*(R2-Rnt?V*(BlPfM)*(73^(R5»R6)»CKXb(SlG)) 

UB2P(M)^(Z«t(R5tR6)»SKXB(8lG)))*TU*CP(MJ*(Z6+(R9-Rt0)*CKXA(SlG>) 

LSM  - SY^f^tTRlC 


WITH  E-MOOES 
lNT6«t«KX*IlP«l2P*Ky9»2 
WITH  H-HUOES 

IF  aH*E0,2)  INTGRL»KY*I1P+kX*KV*i2P 
IF  (IKT.FQ.I)  INTGRL«ItP/S0R2 
IF  (IKT.EU.O)  IN(TGRL«1MTCRL/KT 
INTGRL■1^^TGKL/^JP8(^') 

return 

l8m  - antisymmetric  and  use  • symmetric 

IF  (A3.6T,l,E+2S.0R,A«.GT.l,Ff25)  GO  TO  290 
ZJ«(AU-R3-R9JFCKX8(8IG) 
22«CTVl9Al2-R3tRa)*SKXP(8lG) 
23«(Tvl9Ai2-R3tR9) *0^X81810) 
2<l«(R3^Rtt«AU}9SKXH(SIG) 

60  TO  270 

IF  (A3AA9.GI.I.E*50)  go  TO  260 
IF  lAa,GT,l,E+2b)  GO  TO  250 
Bl«8KAPD(N)F*2/KE 
B2»B1 

B1bBI«CKXB(SIG) 

B2>B2«8KXB(SIG) 

B3«"BI 
B9Bb2 
CO  TO  270 

Bi»»8kAPD(NJ**2/KE 

B2*bl 

B1>B1«CKXB(S1G) 

B2Bb2«SKXB(SlG) 

B3BBI 
Ba>«B2 
CO  TO  270 
BI«0. 

B2«0. 

83*0. 

B9*0« 

IF  lA7.CT.I.Et2&.0R,A6.6T,t,Et2&)  GO  TO  260 

Z3s(R7tRe*Al3)FCKXA(8IG} 

Z6«(R7-R8-TulAAia)»CKxA(8IG) 

CO  TO  310 

IF  lA7*A8,CT,l.Ef50)  GO  To  3OO 
IF  (AB,GT.l.Ef25}  GO  TO  290 
B5«»CKXA(SI6)*8KAPAl(M9»2/K 
B6B-BS 
GO  TO  310 

B5«CKxA(8IG)F8KAPAL(N}**2/k 


BbaeS 
GO  TO  310 
300  B5>0* 

BbBO* 

310  IF  (MODE.LQ.SHLSe)  GO  TO  33o 

llP«R2»Hl+TUl»(felP(M)»{21+(R5+R61*SKXB(8IG)) 
UTVl»E2P(M)Fa2*(H5-Rb)*CKXfl(8lGT)*8lGNCl.fKEn/EPS 
2tTUl*FP(M)*(ZST(R9tRlO)*SKXA(Sl6)) 
l2P»TU*(RlTR2)iTV»TUl*(tlP(M)*(nf(R5-Rb)*SKXb(8lG)) 
l+TVl*E2P<M)*(Z<t+(R5^R6)*CKXB(SlG))*8lGNU,»KE)) 
2+TU*TUl*FP(M)*(26+(«9-Rl0)*8XXA(Sl6)) 

C 

C L8M  • antisymmetric 

c 

C WITH  E-MODES 

C 

INTgRL«AJT(-KX9IIP4I2P*KY»*2) 

C 

C mITH  H-mooES 

c 

IF  (LRtEQ.2)  IN1CRL«AJ*(-kY*IIP.kx*KY*I2P)  ■ i 

IF  (IKT.EO.I)  lNT6RE«-AJFItp/80K?  1 

IF  ClKT.eo.O)  INTCRl«1NT6RL/K1  I 

intgRl«intgrl/npa(m)  i 

return  I 

320  IF  (ISYM.EQ.IHS)  GO  TO  230  | 

IF  CI8YM,|Q,1HA)  GO  TO  l«5 

330  nP«Rl+R2+TVlFB20P(M)*(23t(R5-R6)»8Kx8(8I6))  1 

UBl0P(M)*(2<l  + (R5+R6)FCKXB(Sl6))tTUl»C0P(M)*(26t(R9-Rl0)*8KXA(SlG)}  J 
C I 

C USE  • symmetric  1 


C WITH  E-M00E8 

c 

intgru«aj*i;pfky 

c 

C WITH  H-MODES 

C 

IF  (UR. EG, 2)  INTGRLi-AJFKXFilP 
IF  llKT.ECi.n  INTGRU«8QR29Aj*ALKFllP 
IF  (IRT.EB.O)  INT6RU*INT6RL/KT 
INT6RL»INTGRL/N0P8(M) 

return 

3R0  IlP«R2-RltTUlF{TVl*E20P(M)F(Z3+(R5-R6)*CKX8CSIG))*8I6Ntl,,KE) 

1-EIDP(M)*(2R+(R5+H6)*8KXB(8IC)))+FOP(M)»(26-KR9-RIO)*CKXA(8IG)) 

C 

C USE  • ANTISYMMETRIC 

C 

C 

C WITH  E-MQOES 

C 

1NT6RU»-I1P*KY 

c 

C WITH  H-MOOES 

C 

IF  aR.EQ.2J  INT6RUsKX*HP 
IF  (IKT.EO.I)  INTGRU»-SGR2FaLR*11P 
IF  CIkT.EO.O)  INT6RU"INTGRL/KT 
INTCRu»IN1GRU/NDPA(M)  202 
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subroutine  C0N8KV  (SUf32l«NMOOtl8l6tlM) 
C 

C CONSERVATION  OF  ENERGY  CHCcK 

C 


C N0TE>  does  NOT  guarantee  THAT  ANSWER  18  CORRECT,  ONLY  GUARANTEES 


c 

e 

SELF»C0N8I8TENCV. 

If 

C 

OEFINITION8> 

cc 

811 

m 

FEEOGUlDE  SELF-REFLECTION  SCATTERING 

C 

BLOCK 

C 

821 

m 

FEEOGUlDE  TO  SPACE  MODE  VOLTAGE 

C 

TRANSMISSION  COEFFICIENT 

c 

NMOD 

• 

NUMBFR  OF  FEEOCUIOE  MOOES  IN  uNiT  CELL 

c 

IM 

m 

MINUS  the  number  UF  DIGITS  TO  I^HICH 

c 

CONSERVATION  OF  ENERGY  IS  APPROXIMATED 

C 

BY  solution 

C 

«««««*$♦««««*»«««««*«$««««««««*« A***************** 

Integer  sic 

complex  SU(20f20i«82l(250f20j«8(20f20)f YfYA 
common  /CN8RV/  V(20)«YA(250) 

IM«*10000 
DO  iSo  lAslfNMOO 
DO  140  IRbUNmOO 
8(1A»IB)>(0,»0,) 

00  100  SIGBltlSlG 

8(IA«lB)«S(IA(Ib)Y821(81GtlA)*C0NjG(YA(SIG}AS21(S]6*lB)) 

100  CONTINUE 

00  llO  NNBifNMOD 
OELRBO, 

oelsbo, 

IF  (NN.EQ.IA)  OELRbI. 

IF  (NN.EQ.IB)  0EL8B1. 

8(IAtIB)«8(IAfIB)*(0ELR+8n (NN.IA))AC0NJG(Y(NN)A(0ELS-Sn (NN»I8))) 

no  continue 

AbCABS(S(1Ai IB)) 

IF  (A,GT,1,E-20)  CO  TU  120 
lB»1000 
60  TO  130 
120  laALOGlO(A) 

130  iMBMAXOdf  IM) 

140  CONTINUE 
130  CONTINUE 

return 

End 

#F 
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8U8H0UTINC  PLTCAL  (A*B*0*F#eP8»A1 ♦Bl »DX,OY »SEP.PT,80f LObEt I8THT» 
lJjt8T,8PH#ICR0*NTH) 


C 

c calcomp  pattern  plotting  Routine 

C max  of  to  CURVLSt  51  POINTS  EACH. 

c 


c 

0EFINITI0N8> 

c 

A 

• alpha 

c 

B 

• beta 

c 

0 

• delta 

c 

F 

- FREQUENCY  (GHZ) 

c 

EPS 

• RELATIVE  PERMITTIVITY  OF  SLABS 

c 

Ai 

• x-dimenston  op  guide 

c 

B1 

- Y-DIMENSION  OP  GUIDE 

c 

OX 

• X 6KI0  SPACING 

c 

DY 

• Y GRID  SPACING 

c 

SEP 

• septum  Thickness 

c 

PT 

• POWER  transmission  coefficient 

c 

SO 

• LOBE  number  in  internal  ORDERING 

c 

LOSE 

• numbfr  of  bears  To  be  plotted 

c 

I8TRT 

• ARRAY  PICKUP  VALUE  FOR  CURRENT  PACE 

c 

JJ 

• LOBE  selection  VECTuR 

c 

ST 

• SIN(ThETA)  array 

c 

SPh 

• 8IN(Ph1)  of  plot 

c 

IGKD 

• GRID  TYPF. 

c 

NTM 

• NUMBFR  OF  POINTS  TO  BE  PLOTTED 

i 


i 

i 


-I 

I 


i 

n 

i 

t 


‘S 

'i 


t 

} 


i 


DIMENSION  PTnO*5l),JjaO)isT(5n,STlC5l)»P(5n 
integer  SOUO) 


L«1 

IF  (NTH.CE.Zl)  L»2  i 

Call  axis  **ifS,f o.*o.,o«2) 

call  symbol  (R,3«2.6f .lOtJHslNtO.tSl 

call  GREEK  (a.6t2,55t.lSfO.,8;  ^ 

CALL  symbol  (a,7»2.5t0,07f 1hOiO,,i ) 1 

Call  IAXIS  C2,0»J,0,30hPO*vER  transmission  factor  (DB)»30#^.|90.♦ 
1-30,«5.I  ] 

IF  UOHO.EO.I)  call  symbol  (3,75f9.6f ,lf ISHTRIANGULaR  GRId»0,»15) 

IF  (16RD,E0.2)  call  SYMbOL  (3,75i9.bf , 1 t IPHRECTANGULAR  GRlD,0.»l6)  1 
CALL  PLOT  (2.#3.»-3) 

13»I8TRT-1 

DO  100  l2slfL0BE  t 

I3*13tl 

14>«JJ(I3) 

ni«nth 

NaO 


Il»i 

90  00  91  lallfNTH 

IF  CPKUf IJ.GT,«32,R999)  Go  TO  92 
Nsl 


91 

92 


93 

94 


continue 


NaNtl 

IF  (N.GT.NTM)  go  to  9b 
00  93  laNfNTH 

IF  (PT(I«»n.LE.-32,«99999)  GO  TO  99 
Nial 

continue 


N2aNt»Nf i 
00  95  IsNtNi 
P(I)«PT(I«»I) 
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STl(I)sST(n 

CONTiNUt 

P(Nltn*-30,0 

P(Nl*2)*5,0 
8T1  (Nl  + nsO, 

STUNl+<i)sO,2 

CALL  Llf'iE  ISTl  (N)  »P(NJ  ,N?,1  ,L  * It'+i  J 

IF  (Ml. to, NTH)  (.0  TO  V6 

IiaNl+t 

GO  10  90 

CONI iNUfc 

CONI  I:\iUt 


CALL 

PLOT  ( 

•2, ♦ 0 * ♦•3) 

CALL 

symbol 

(2.5f",7bf,lf  aHi-  s 

♦ 0, 

♦ ‘'I 

Cali 

MUN  her 

(999,t999,.,l .F,0.. 

2) 

Call 

symbol 

(999. f999, , . 1 .an  OH2.0 

,*U) 

Call 

symbol 

(5.bf»,7bf,l.yHA  = 

,0, 

♦ «) 

call 

MJMBfcR 

(999,f999,«,l.Ai.O. 

♦3) 

call 

symbol 

(999, f999, 1 IN 

• .0 

• ♦ A ) 

Call 

GREEK 

(2,5t-l,05t.l5#0..1 ) 

call 

SYMBOL 

( 2 . 75 ♦ • 1 • 0 ♦ , 1 . ?H=  *0,. 

3) 

call 

number 

(999,t999,.,l ,A,0,, 

3) 

CALL 

symbol 

(999,.999,,.l,yH  I'-’ 

• « () 

, ♦ A ) 

Call 

symbol 

(b.bt"l,f,lt9HB  = ♦ 

A) 

Call 

number 

(999.f999,^.l,Pi,0, 

♦ 3) 

Call 

symbol 

(999,t999,f.l tan  IN 

. ♦() 

. ♦ A ) 

Call 

GPtiK 

(2,5fl.30f,15f0.f2) 

CALL 

symbol 

(2«75f«l ,25f , 1 »?H* 

♦ 0, 

♦ i) 

Call 

number 

(999, #999. ♦,! fR.O. ♦ 

3) 

CALL 

SYMBOL 

(999,f999..,l,aH  IN 

, . 0 

. ♦ A) 

Call 

symbol 

(5,bt»l,25f,I.tH0»0 

• ♦ n 

CALL 

symbol 

(999.,-l,30f.u7,lHX 

♦ 0. 

♦ n 

CALL 

symbol 

(999, ♦•! ,25f , 1 ♦ 3H  s 

,0 

.♦3) 

Call 

number 

(999,i999,,.lf0xtU. 

♦ 3) 

CALL 

symbol 

(999, ♦999,,,!, Am  IN 

• f 0 

. fA) 

CALL 

GREEK 

(2,Sf-l,55f,lS.0,,a) 

Call 

symbol 

(2.7bf»l,b0f,l f?Hs 

’ (^ . 

♦ 3) 

call 

NUMiBLR 

(999.,999,..l ,0,0,. 

3) 

CALL 

symbol 

(999, i999, , , 1 ,Ah  U 

. ,0 

• ♦ A ) 

Call 

symbol 

(5,bt*l,5f,l.lHp,o, 

♦ 1 ) 

CALL 

symbol 

{9?9,f-l,b5t,07,lHY 

»0, 

♦ 1 ) 

Call 

symbol 

(999, f-l ,s» , 1 ,5h  = 

♦ n. 

.3) 

Call 

number 

(999.f999,,.l,DYf0. 

• 3) 

call 

symbol 

(999, f 999, ,.l .Ah  IM 

, » 0 

. ♦A) 

CALL 

GREEK 

{2,5f-l.80f.l5,0.t5) 

Call 

SYMBOL 

(2.75f“l,75.,l»^H5 

♦ f' « 

♦ 3) 

CALL  NU^^.fltR  (999. ♦999, 
CALL  SYMBOL  (5,bf-l.75 
Call  numblp  (999, 1999, 
CALL  symbol  (999, .999, 
call  symbol  {3.bt-2.f, 


call  symbol  {3.bt-2.f,lfl9HcOl  p 
CALL  GPtFK  (a,9f-2.05i,r..  ,0,.<?n 
CALL  symbol  (a.95f-i’,o.,lf3f.  = ♦ 
Call  NONBER  (999,f999...1,SpH,0. 
CALL  symbol  ( 1 ,bf2,?bf  , 1 .6Ml  LGt 
call  symbol  (1 ,75t-2.B5f.l 

CALL  GPttK  (2.2bi»2.5..l5.0,,lBj 
CALL  symbol  (2.35f-2.S., 1 .5h  = ♦ 
AS»SO(ISTwT) 

call  fjUMBfcW  f999,f999,,,  5 .ASfJ'vf 
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,1  .EpSfO,,p) 

, 1 * SuStR  a ♦ 0 , ♦ b ) 

• 1 f SgR I 0 , , 3 ) 

,1.9h  1N,,o.»4) 

♦ 1 AHcO  I Pi  AN't 

r.. , 0 , , «?  n 

,lf3f‘  s t0,t3) 

. 1 f SPH.O, ,3) 

, IfbMl 


S 1 N ♦ 0 , . I y ) 


0.  ♦3) 

‘3) 

’'1 0 ♦ I'  ♦ ♦ B ) 
-1  ) 


0 . ♦ 3 ) 
0) 


TF  (Ll'Hfc,Lf,2J  bU  Tl.  no 
CALL  SY.’MUL 

CALL  bf»tFK  (5,b»-2.S».1Si(',,lH) 

CALL  SY^i^tiL  (S.e*-2,5»  * 1 f = *0,,^; 
ASSSU(IST(^T  + i ' 

Call  •jUMhE.W  ,As*0,*A) 

If  anMfc.LT.3J  bU  T()  Ho 
CALL  SYf'MuL  n.75t-2.7n,.l,u,u.,-5  ) 
Call  b^fct^  (2«25t»?,75t , ist n. 1 1 hj 
Call  SYN'HOL  (2.35t-2.75t,l»XH  = ,(),.3J 
ASSS('(  ISTWT+2) 

CALL  'v.Ur'MfcR  (99S(,,999,,.1  , Ac,i,.  ,0) 

IF  aiiMt.LT. a)  bn  TO  no 

Call  bV-'HUL  (5.*-2.7».l 

Call  GwfcfcK  (5.Mf2.7S*  .is.o,,  18) 

Call  SY^'■^UL  (5.e»-2.75»  . l . 5m  s .0,.i) 
A3  = S0( 15TW  f + 3J 

Call  (999.»999...1  .At;.0,»A) 

Call  MU'!  (m.M.-3.o,-3) 


,s 


! 


i> 

ii 


SuHKOUTl.'Jf:  FOuRmOS  (F»HB) 

C 

t SPbClAL  OlSPfcWSIO'')  RtlATin^  SOLVER  FOR  RAPID  COMPinATTON  OF 

C ELt^'twr  mismatch  at  bPCAuSlOfc, 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


deficit IUMS> 


F 

• FWEuut'wCY  Iiw  GHZ 

t»8 

• X-D1RFCTF,L  rtAVENUNtfltR  VECTOR 

A 

- ALPHA 

H 

• beta 

0 

• DELTA 

dBb 

• G'lUE  height 

IPI 

- 6.?»Sl*.,,,, 

EPS 

• relative  PERHiTTIvnr  OF  SLABS 

OlMfcf'i{>IUW  BH(l) 
RfcAL  K»KF»t<I.’MC 


CO^'''iUiM  /aAVOD/  Af  rJ»U*FBH*  TP]  .fcPS 

Data  L/n  .Hugascb/ 

fPlL*1Pi*F/C 

FpSlstP5«l.y 

AJSlPll  *4 

BJSIPII ♦B 

Dl  = lPU 

LsO 


Di)  17i;  |»|||QFS1  »a 
DO  lb«)  f'Sl»4 
LsLtl 


sisl/fj 
Ssl «01 

If  iK.u.o.i  s»v.v9 

Ki^Csy.JlA 
KsK-k  ii'ic 
TsO 

loo  KaK  + ts]NC 

DlFFsuIFFl 
IsI  + 1 

KESfcPSl tKAAbSCKj 
SsSwMT  (ABS(»vE  ) ) 

IF  (Kfc.LT.O.)  Ss-b 

kf:*s 

DJFF  lsr>iSFMM(jpF,«,Kt»Al  »8l#tU  »FPS) 

IF  lI.EU.n  GO  10  1?0 
IF  IKlliC.LE.l  ,t-lO)  GO  10  I3C 
If  lABSCniFFAi.aFFn  ,L£  , I ,£•!  0)  GO  Tt  130 
If  IDlFT  »l)1FF1  ) n0,130fU»n 
UO  K*K«Kjmc 

KINC»U,5*KINC 
DIFF IxDIFF 
Go  10  loo 
UO  DlFF«olFFl 
GO  TO  lijO 

liO  IF  lAbS(01FF ) ,LT  .ARStOIFFj  ) ) GO  To  140 
BbCL)*K 
Go  to  loo 
l«0  KeK-kimc; 

Hb(U*k 

loo  COMllNUt 


■I 


I i 


i 

s 
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SuBKOuTlNt  KfltT  Cw^^0^tS) 


I 


SPtClAU  SUtJPuUTIwF  FuR  CREATING  LSk(NtO)  AixO  LSEC^'*!)  K-HhTA 
DlARRAy.S 


DFf  Ii'Jl  I IONS> 
iviNUOtS 
HG 
KA 

IMUUt 

A 

6 

0 

BB 
TPi 
tPi 


nu'-^!>Y  akqu'^E'NI 

ng«!  At  TZFD  LtiNGnuniNAI.  !/> /WtNt"‘‘t'.Ek 
KFA/p 

MUDk  pk&TG'MATlONS 
AlPhA/LAmbDAO 
BhTA/L ANHDAo 
HkLl A/LAMhOAO 
((•"IDf  HFIP^T  J/LAhHOAn 
ft. ...... 

«kl.  AT]vt  PFK^iTlIvTTV  OF  SI.ARb 

OlMtNSIUN  RaM(2,2,2) ,so(?,2,?).SF{2*?.^) 
real  K0»KA.K«KE»w2R»Klf.'C 

CONMIJN  /KhETA/  BOlbtlun  .^AduD.r'lUPtCB) 
common  /aAVGO/  A.b.OtCtPtTPT.FPS 
Data  c/n  ,802B52b/ 


AASAfH+0 
F0sC/( rPTFAA) 
DFS0.03FF0 
NpsiOl 

FPSIsePS-1 .0 

FpSUSSOWTCFPSl) 


F xFO«UF 

wWllt  (b.ROU)  (iMUDt(l)»I»l»bJ 
nn  200  IFsl  ,NF 


FsF+OF 
K()X1HI»F /c 
KA(IF  JxkOFAA 
AIskOXA 
B l»KU»B 
DJXKOFD 

N28ato,b*C/(HB*F ) J*f2 
DO  180  '^OUtsl.2 
mi«MC0F. +2 
K«»tPSO+ I »t«l 0 
Mo*l 

00  IbO  NSI ti 

Jsl/N 

Ssl .01 

IF  (k.Ll.U.J  SS0.V9 

KsJ*K-f{l»J)FSFK 

KlNCsO.il  A 

KsK-K  i^■C 

IsO 

loo  KsKtKlNC 
DIFFsoiffi 
Isitl 

KEStPSl fKFABSCK ) 

S5SuRHAUS(^t ) ) 

IF  IKt.LT.O.)  Ss-S 
KE  = 5 

DIFF  l=f>IS<'’(Ml  .KfKt.Al  .R1  »01  ,FPSJ 
IF  (I.FL.l)  GO  TO  120 
IF  (KI-NiC.lT.I  ,E-10)  go  in  130 


i 

■i 


1 


IP  CABS(r)IFP»Dlf-Fn.LE.l,fe«lO)  (,n  fu  i50 
IP  (01PP*DIPFI)  UOtUOtlOO 
no  KsK-KlMC 

KINC«0.5PKIiMC 

dipfuoifp 

60  TO  100 
120  DIPF*DIPFI 
GO  TO  100 

liO  IP  (Ab8CDIFP),U  .ABSTDIFFD)  60  TO  1«0 
60  TO  150 
l«0  K«K»K1NC 

Ga&P8lfK«ABS(K) 

8k8UFT (ABS(G) ) 

IP  (G.LT.O.)  Sa-S 
KE«S 

150  80(^^OOEf^ifM0)*K 

SE{f^ODE»N|Mo)*Kt 

6«1 tO^KPABBCK) 

SsSuRr(ARS(6)} 

IP  (G.LT.O.)  S»*S 

GAP(WOOE»NfM0)sS*KA(IF ) 

UO  CONTINUE 
180  continue 
II»0 

00  190  ll*l*2 
00  1^0  I»l»2 
Iialltl 

IP  (II.EO.A)  GO  TO  190 
B6CII»IP)a6AM(If Iltl) 

K«SO(Itll»U 
G2««K»ABS(K)fl 
GaSQRT(ABS(G2) )»KA(IF) 

IP  (G2.LT.0,)  Gs-6 
BG(n  + 5»TF)*G 
190  continue 

write  (b»910)  KA(IFJ»(B6{I»lFm«l»6) 

200  CONTINUE 

C 

900  format  (1H1»5X»10H  »vA/2  * 1 0HXGAMMA»A/2  » / 1 SX  » 6 C 2 * A 7 t 1 X ) ♦ / / J 

910  Format  (8XfF5.2f2Xf6(2X,P6,3»2X) > 

END 


90 


8UfiK0uTi^Jt  I SfeMuO  f 

c 

C COMPuTt  ust  mHuE  FuMLtlO'\  ^ViX) 

c 

C OEP  liMll  IUI'<S>  SEi  aiJF.Ri  U1  T.iE  LSMuSt 

0IMt^8IU^  Vx(iM) 

REAL  KAH*KAPtt‘'»Pbi^PA»^t/PS*.\.nPA 

CcwMOfy;  /AkRaY/  AL*Pt*l'L*BU  ♦SbPtL  ,TPi  ,tP8»Sl«) 

/CUE>  3/  blB(  10)  *b?Mf  1 0)  »h1  CPU  u)  fB«il)P(  10)  »CP(  lo)  fCOP(lO) » 
iFiPCUO  tFdp(io)iEiPpn())»k?hP(jo)  ,FHno)»Enp(io).kPsuo)»^PA(iO)» 
iM)PS(  If'j  , U‘)  »K0U«.T(?1)) 

Cuf^^Ci'v  /vi(.')kS/  KAPtiJO)  tftAPF  f ?J)  *r,AV(?u  ) *'^Ul>b  (^0)  * iSY'^'tJO)  f C,\(  ^0)  ♦ 
iMMCiO)  .i*ncOHir)(?o) 

M*MM(\l) 

Asrpl»(ftL+BL+OL) 

xntL»o • 0?* A 

Pa-lPpt!L 

Oa-IPl*(H|,YL>L) 

NISKOU^I  (I'.) 

Bl s-8o,»B/A 
ni s»bo.*o/A 

If  (ISYC(c)  ,F'0,1ha)  C,l  10  luO 

Xs-1  *(12*A 

On  I3u  Uitbl 

Xsx+xuEL 

Y*AoS(X) 

IF  (X.GT.U)  GO  TO  no 

vxcnscoPC'vii  )*sr^t(KAp  (I,)  ♦(<♦«)  )*siGcn . »RAPtm  i/kdpsc^o 
VX(10i-I)*VX(I) 

GO  10  130 

no  If  lX.Gt.B)  GU  10  1?') 

TsKAPE(OU(b-X) 

vx(i)3(bioP(ci  )*cuS(.n  i-p^.)p(  vi  )*sT^c(T)  )/i'‘uP&(M ) 

VXU02-1)»VX(1) 

GO  TO  130 

1^0  Vx(I)«COSC(KAP(i\)YY)/f,r.pSC-  1 ) 

YX(lOi-I)*Vx{I) 

130  CONTIwUt 
GO  TO  IbO 
IRO  Xs-;,02»A 

no  i7o  i«i»bi 

XaXtXDEL 

YsAbS(x) 

IF  (X.GT.O)  GO  TO  150 

VXU)aFbP(M  )^SiNC(KAP(  \)»f  x + A)  )/^')PACfil ) 

VxU02-i)s-vx(I) 

60  TO  170 

150  IF  lX,G1*b)  GO  H)  IfeO 
T*KAPEf^)A(b-X) 

Vxa)s(tlUP(M  )^COSC(  T )-E20p(M  ) OS  I >'C  ( T ) ) / "luP A f M ) 
VX(102-I)S-VX(I) 

60  TO  WO 

160  Vx(I)»-SlNC(KAPCN)»Y)/M.)PA(!n  ) 

VX(102-I)»-VX(IJ 
170  continue 

c 

C call  PnINT/PUOT  KOUTiMt 

c 

160  CALL  PRiyT  (VXtjHLSFfB!  fOU 


SUBKOUTiNt  LSMMOO  (Ni) 

c 

C COMPUTt  uSM  mode  function  hVOO 

c 

C UEFINITIUNS>  SEE  SUBROUTINE  LSMtSE 

REAL  KAP.KAPEfNPSfMPA»NDP8*N0PA*TXtl0l) 
common  /ARRAY/  AL»BL»UU»81L,SEPTL»TPIfEPSfSl«) 

Common  /cuefs/  BiPtio) .a2PdO)»BiDPCio) iB20PUO) tCPCio) »cdpuo)» 

1EIPU0)»E2PC10)  tElOpnO)  »E2oPnOUFP(10J»FDP(10j*NPS(10)fNPA(lO)# 
2NOP8(10)tNDPAUO)»KOUM  (20) 

common  /MOOES/  KAP(20) ♦KAPE(20)»6AMC20)fMuDEC20)»ISYM(20) »NN(20)f 
1Mm(20) »MODORO(20) 

MsMM(N) 

A«TPU(AL  + 0L>DU 
X0EL«0.02*A 

8«-TPl*bL 
Oa-TPl*(BL+OU 
N1«K0uM  (N) 

Bl«-bO.»H/A 

Ol*-bO,*n/A 

IF  (I8YM(n),EU,1ma)  go  10  laO 

Xc«i .U2*A 

DO  Hu  Isltbl 

XsX^XoEL 

YsAbS(XJ 

IF  (X.GI.O)  GO  rO  110 

lx  (l)sCP(M)  40080  IKAP  (f' J ♦ (X^A)  )/NpS  INII  ) 
ixno2-n«ix(i) 

GO  TO  130 

no  u (X.GT.b)  GO  TO  UO 
T«KAPt(i04(b»XJ 

IX(I)a(blP(M)4C08C(T)-b2P{MU4STNr.iT)  )/NPS(Ni ) 

lxn02-l)5lA(U 

GO  TO  130 

UO  lx(i)sC0SC(^AP(N)4Y)/r-P8(M  ) 

IxUOg-nalxd) 

130  continue 
60  TO  160 
UO  Xa*UU24A 

DO  l7o  lalfbl 

XaXfXoEL 

YaAbS(XJ 

IF  (X.GT.O)  GO  TO  ISO 

IX(l)aFP(M)4C0SC(KAP(NJ4(X*A)  )/NpAlNl  ) 

IxU02-l)a-lx(I) 

CO  TO  170 

li>0  IF  (X.GT.b)  GO  TO  UO 
TaKAPt (N)4(b-X) 

lxCI)a(|tP(Nl ;4L0SC(T)-t2P(M  J»STNC(T) )/NPA(ni  ) 

IX(102-I)8«1X(I) 

GO  TO  170 

UO  IX(1)sSINCCkAP(w)4Y)/nPA(N1 ) 

IXU02-1)8-1X(U 
170  CONilNUt 
C 

c CALL  PKINT/PuOT  KOuTl^fc 

c 

UO  Call  h»nt  (iXf 3r’LS^’»8i .01 ) 
return 
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I? 

I' 


SuBHOuTINE  PRNT  (VXflMtBlfDj ) 

c 

C PRJNT/PIOT  ROU1INE  FUR  MoOf  FUNCTIONS 

c 

C DEFINITIONS> 


c 

vx 

m 

NUDE  amplitude 

C 

IM 

m 

M(iOE  TYPE 

C 

bl 

m 

NURMAlIZED  DIELECTRIC 

ROUnOARY 

POINT 

C 

01 

m 

NURMAlIZFD  DIELECTRIC 

80U*v'DARY 

POINT 

DIMENSION  VX(l)»IXU01)  ♦IP(i02)»TC(5«) 

integer  ouTtaLAwKt Plus* zero. slash 

DATA  OOT*BLANK*PLuSt'^iNUS»2rPO*SLA8H/lH,,lH  * IH+ ♦ iH-t  IHO  ♦ IH\/ 

IF»3He-Y 

IF  (IM,E0,3ML8MI  1F*3hH-y 

iblsbi 

lOl«Dl 

I8R*52+IBL 

IBL>82-1BL 

IDR«52fIOl 

IDL«S2-IDL 

XMAX*0« 

DO  100  1*1*101 
T«AB8(VX(D) 

XmAX«AMAXI (XMAX*T) 

100  continue 

WRITE  Cb*ROO)  IM,IF 
DO  no  1*1*21 

write  (6*910)  (JJ»5l*VX(JJ),JJ«l,jOU2l) 

no  continue 

Do  120  2*1*101 

T*50,fVX(I)/XMAX+0,25»SlGN(i,*VX(l)) 

1XU)«T 

120  continue 

write  (S*920)  XmAX 
IP(l)«DOT 
DO  130  1*2*102 
iPCnsBLANK 
130  CONTINUE 

IpaOL)«8LA8H 

IP(I8l)«SLA8h 

IP(1BR)*SLA3H 

IP{IDR)*SLA8H 

DO  180  1*1*51 

I1«I*1 

Ifi*10o-2AI1 

B*0,01«IB 

NaO 

Do  140  JJ*1*101 
IT»51-IABSCIX(JJ)) 

IF  (lT,NE,n  GO  TU  IRO 

NaNfl 

IC(N)sJJ<f1 

lP(JJfl)*PLU8 

IF  (IX(JJ),LT,0)  lP(JJfl)*MtNu8 
IF  UX(JJJ.EQ,0J  1P(JJ  + 1)*ZeRu 

190  continue 

IF  ((ll/10)*l0.E0tH)  GO  TO  150 
WRITE  (6*930)  (IP(JJ) *JJ*1 *i02) 

60  TO  160 
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IbO  wRiTt  (bfRaiO 
160  IF  (N.eu.O)  GO  TO  180 
DO  170  JJslfiJ 
lR>lC(jg) 

IP(lB)*bLANK 
170  CONMNUt 

IP(iDL)«SLASH 
IP(.t.Bl)5SLASm 
IP(IBR)»SLASH 
lPCiOk)»SlASH 
160  CoNfiNUt 

WRIIE  (6i*>50)  Cl-blflaltll'ltPo) 

Return 

c 

900  FOPMAl  UHlf5XfA3»8H  NOUtt  (»A3»18H  CUMPONtNT)»LANBOA*//f 

15(10H  lOu^X/A  tlOH  V^LA^^bPft  )♦/) 

910  FORMAT  {5(3X*I«f3x«F8,A»ax)) 

9i0  format  (IHl  ♦1oX»7HVmAx  s .FjO.b#//) 

930  FORMAT  (loX,io2AU 

9R0  FoMMAT  (5X.F5.2»102A1) 

950  format  UlXtlH*tl0(10M,,,,,,,,.»l|/»U/X»T3*5(17x»I3)»/t58Xf 

17hIOO*X/A) 

END 

•» 


oooooooooo 


subroutine  creek 

SUbROUTINE  TO  PLUT  GREEK  CHARACTERS 


DEFINTI0N8> 

A -X  location  of  character 

Y - V LOCATION  QP  CHAkaCTER 

H • Character  heicht 

T • plottinc  ancle 

.4  - character  NuMbER  (SEQUENCE  IS  GREEK 

ALPHAHETT 

DIMENSION  KUiO)»L(2Sj 

Data  K/77ai»a2?*>f i5o« *21  ♦‘457^1 ‘>26»S6aS, 5^77*1  «3«»R3«£« suit 

n277*770<if  l5?2»77lO,22R‘)t  ;7?a,lil?t25  32»3Jlb»l62K»i57  7f7703f3j77» 
2351B*R02,1  131,77U»21  Jl»1356»i^2Rtlto77»  770  l»I31Uf377»l3i5»aR20f 
3771  i»l«25i35«RfR3l3, 1221  ♦31fl2t«377W725t  1 1 21  .2277f77U*l*577,abli» 
a«l77f77uii2577f  I625t3l«l  f2tt77»R*i32»3l  77f  1221  »3?77,770bf  1511  »3«3b. 
S77l0«2U31ill02f l3URf l52B«3R77vl333*77l3»2R3R*R3tt2»3l?l*l2l3*7703t 
61tt54»77M«f3077»102R,1324*3RR3»32l2»775«*2Rl3*l22l»3lR2.R33tt»7703f 
71 a3a» 3377,241 o»7703»l«2a» 1120,3041 *5477 , 7720 • 3577 , 342a, 1 3l2» 2 l3l f 
b4243* J477,4477 *4 14 *3040 *771 0*2477,515 *312 »2233»4577,77 14 *302, 1131* 
V4243*34/7,2123/ 

data  L/1 *7* 14, 16*24, 29, 34,39*46*49 *53*57 *6? *65 *72,77 *a2»RS* 90, R4, 
199*106*109,115*121/ 

CALL  CALCKP  (xF»YF*lOUM*-4) 

Hx*H»xF/6.0 
HY»M*YF /6,0 

TT*U,0174533AT 

s*si'nTn 

cscus(n) 

ChX«C*hx 

ShX«S*HX 

ChY«C*HY 

ShY«8*HY 


IZ»99 
lA«LCi'i) 
lB*L(Mtl )•! 

CALL  CALCNiP  (X*Y*00,n 
DO  3 I«lA,Ib 
JaK(l)/loo 
DO  3 II«1,2 

IX«J/10 

IE  Ux-7)  2*1,2 

1 IZ*00 
GO  10  5 

2 IY«J-10*IX 

CALL  CALCMP  (X>CMXTlX«SMY*lY,Y4CHYAIYtSHX1‘IX*l7,n 

IZ»RR 

3 J«K(I)»100YJ 


CALL  CALCMP  {X4C<‘h»XF*YtS*H*XE,0n,U 

RETURf'j 

END 


if 
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or»oooor»r»o 


SUBHOUTINE  AXl8(XPA6t»YPA6EfIflC0»MCHAR»AXLfcN,ANGLE»FlR8TV»0ELTAV) 


• • • • • 

CONTAINS  CALLS  TO  CALCmP  INSTEAD  OF  PLOT 

CONTAINS  lAXlS  ENTRY  PqINT  • WRITES  VALUES  IN  INTEGER  FORMAT 
XPAGEtYPAGE  coordinates  OF  STARTING  POINT  OF  AXISi  IN  INCHES 

• • • • • 

IBLO 

AXIS  title. 

1 1 • • t 

NCHAfi 

number  of  Characters  in  title,  t for  c.c»w  sioe. 

• • ■ • • 

AXLEN 

floating  point  axis  length  in  inches. 

• 1 1 • • 

ANGLE 

angle  of  axis  from  the  x-oirection»  in  degrees. 

• • • • • 

FIHSTV 

scale  value  at  the  first  tic  mark. 

• t • • f 

UELTAV 

change  in  scale  BETv'EEN  tic  marks  One  inch  apart 

•%•••  UCUTAV  U 

oiMENsiuN  lecud) 


lb  NDECB2 
InT«2 
GO  rn  a 
c 

ENTRY  lAXlS 

C 

IP  (Aa8(n£LTAV),LT.l.}  GO  To  lb 

NOECBol 

InT«1 

B KNbNChAH 

V2CW»1, 

Aal.n 

lFfKN)lt2,2 
1 As*A 

i EXBO.O 

AOX«AbS(OfcLTAV) 

lF(ADX)3f7*3 

3 1F(AOX*99,0)6i«»<» 

4 AuXsAQX/lU.a 
Ex>EXtt  .0 

GO  TO  3 

b aoxbaux*iu«u 

EXBbX-1  ,0 

6 IMAnx-0,J)5»7«7 

7 XVALBFIKSTY410.U**(-EX) 

AOXbOELI AV*10,0YYC-l A) 

STHBAN6LE»0,0l7i»535 

CTHbCOS(STH) 

STHbsimisTH) 

CTHTIcbCTh*VI£w 

STHTIC»STH+vIEa 

PXB*-0,l 

POSNB.lb 

iFTvlE'^in  .f9JPU5Ns,2b 
OYBbPoSN4ia»(),05 
XnbXPaGE+OXB^CTm-DYbASTh 
YN«YPAfiE*DYej*CTM  + UXHA5Th 

nticbaxlfn+uo 

MTbnTiC/2 

DO  2o  IslfNTlC 

OXNbO, 

oynbo, 

GO  TO  (9»10)»1NT 
9 NbGbO 

IF  CxVaL.lt .0,)NfcG*l 
AXVALsAdS(XVAL) 

NDIGbB-'vEG 

IFCAXval,6E,9,5)ni;TG8  1-.nEo 


219 


If-  (AXVAL,GE.99.b)NDir,*0 

nxN»MDl6»0.0«5*CTh 

DYN*f'iUlG*U*095»Srh 

10  CAU.  l'lU^'Rt9  (XMf0XK|Y:«!-M;V".»n,U'bfXVAI.*ANGLtffU'ft3 

XvAL*XVALf ADX 

XMSXNtCTHT  IC 
YM*Yi'JtS7HllC 
lF(iYT)20.iU?0 

11  ZsKiM 
lp{LX)12«l3tl2 

U Z«Zt7,0 

13  DxB«*,07*Z4AXLEn^0.5‘'‘VU-a. 

PoSiMS.Sib 

IF  ( VlE^.l.T,,9)PUSi-vS,a<?S 
DVB*PuSKYA-0 ,07b 

X7sxPAGE*nXMYCTr1-OYb*STH 

YT*YPAGt^LYh*Cfr<-^l;X0YSTM 

Call  StMbOL  txT*  Y 1 fO.lu.IoCi'f  AI'M  F tr  ") 

IF(EX) ia,20»  1« 

14  ZsKN+2 
XTBxTfZYC7MAO,14 
YTSYTfZ^STH^O.ly 

Call  symbol (xt»vi  »i;.i4»3hAin,Aivr,LF»3) 

XiaxTf  (3,0»CTH-o,tt*i;TH)»o, 

YTaYT  + (3,0Fi»TH*O,a»L  ■ rt)*0, 1 /i 
CaLL  OUMHEH(xT»YTi0»O7»FX»A'^'lil.F*»n 
20  NT«nT*1 

CALL  CALCMP(xPAGt*AXLFi-*rTnTlC»  YF  AGPt  AXLt ''<♦«  I T U ♦ 0 » I ) 

nxBs»o,07^AARTH 

OYBB^OtOZAA^CIH 

AaNI lL*l 

XNSxPAGttAALTHriC 
Y\»YPAGE*A»STnnc 
DO  30 

Call  callmpcx^i  Y.vt99»i  3 

CALL  CALCMP(xNtpx«»Yr-*LYb»99. 1 ) 

call  CalCMPCxn*  Yt)f  Of  1) 

XM«XN.CThT ic 
YobYI.-STmT  IC 
3u  COMIi-jUt 
R£TuR/j 
EMO 

• • 
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SUBROUTINfe  CSIMEU(A»N«H*M»kS) 

C SCREEN, CSIWEO 

C CDC6700***C8IMEU 

c 

C PURPOSE 

C UBTAIN  SOLUTION  oF  A SET  HE  ST^»ULTANtOUS  LlNfcAK  EQUATIONS* 

C AXsB 

C 

C USAGE 

C call  CSIMF0( A»Nf 6»f''»KS) 

c 

C DpSCRIPTION  of  HAKAt^tTEKS 

C A • matrix  of  COfFETCUntS  stored  columnwise,  these  are 

C destroyed  in  the  CQMpUTATTr.N,  ThE  SIZE  OF  MfTkIX  a IS 

C N HY  N, 

C B • matrix  of  original  ConSTa^'TS  (LENGTH  n ry  n),  THESE  ARE 

C RFRLACtD  BY  FINAL  SOl»'TTGm  VALDES*  MATkTX  X. 

c N - number  of  equations 

C M - number  of  sets  of  SOL'JTIfNS 

C KS  - OUTPUT  DIGIT 

C 0 FOR  A normal  solution 

C 1 FUR  A SINGULAR  set  OF  EOUAlIONS 

c 

C remarks 

c matrix  a must  he  general. 

C IF  matrix  is  singular*  Si-LUTiom  VALULS  AkF  MtAMuGieSS, 

c 

C method 

C METHOD  OF  SOLUTION  IS  HY  ELlMTNAlluN  USING  LARGEST  PIVOTAL 

C DivISOR.  EACH  stage  OF  F.^ImInatTUN  CONSISTS  UF  T nTE  RLHaNG  I nG 

C ROWS  WHEN  necessary  TO  AvOID  MVISION  BY  ZERO  OR  SmaLL 

C ELEMENTS, 

C the  FORtNARD  SOLUTION  TU  OBTAIN  VARIABLE  N IS  OUNE  I'l 

C N stages,  THE  BAC<  SOLUTION  FOR  ThE  UThEK  VARIABLES  IS 

C calculated  by  SUCCESSIVE  SUBSTITUTIONS,  FINAL  SOLUTION 

C VALUES  ARE  OtVELOPtf)  IN  vFCTUR  B*  wITH  VARIABLE  1 IN  R(l)* 

C variable  2 IN  H(2)*, ,,,,,.,*  variable  N IN  H(N), 

C IF  NU  PIVOT  CAN  be  FOUND  FaCEFOI^'G  A TULtRANCF  OF  U.O* 

C THE  matrix  is  COnSIOERFO  SIN&UlAh  AND  •'S  IS  SFT  TO  ],  THIS 

C tolerance  can  be  modified  HY  REPLACING  ThE  FIRST  STATFmFnT, 

c 

c 

c 

c 

c Forward  solution 

c • 

complex  A,B*BIGA*SAVE 
DIMENSION  A(20*B0J *flC20*20) 

TOL*0,0 

K$SU 

DO  20u  Jsl*N 
Bl6As(0,*0,} 

C 

c SEAR.h  for  maximum  coefficient  in  column 

c 

DO  120  I«J*N 

If  (CABS(BIGA)»CABS(A(1*J)))i10*120*1?0 

no  bigaba(i*j) 

IMAX«1 

120  continue 
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I 


f 

I 

I 


I- 

i 


t' 

•? 

f 


i 


I 


c 

C TF8T  FOK  PivnT  LtSS  Than  fOtfcKANrE  lSi'M(*ULAK  MATRIX) 

c 

IF(ABS(KEAL(BIGA)  J4A8S(AlMAG(tJT0A)  )»TUU  I30t  130«  UO 

UO  KS«1 

RETUHn 

c 

C iN'TtPCHANGE  R0>N8  IE  NfcCtSSApY 

C 

IRO  CONllNUk 

Do  i^>0  I*J»N 
SAVESA(J«I) 

A(J«I)BA(lMAXf I) 

A(IMAX«l)8SAVb 

C 

C DIVIDE  EQUATION  BY  LEADING  cOtEFTclENT 

C 

A(J»I)«A(J»n/BlGA 
150  continue 

Do  Uo  l»l«M 
SAVE*B(J»I) 

B(jf  n«BaMAx»n 
B(lMAX»n«8AVE 
160  B(JtIlsb(J(I)/eiGA 
C 

C eliminate  next  vahtablf 

c 

IF(J-n)  I70fil0»l7l. 

DO  IsJlfN 
DO  IBO  (VaJlfN 

A(ItK)aA(l»K)*AlJtK)*A(l»J) 

DO  IPU  K»1»N 

B(I»Kj*BCI»A)*B(J*K)*A(If J) 

continue 

BACK  SOLUTION 


1 


I 


3 

i 


T 

I 


i 


j 


NV*N»1 

DO  22y  L2*lfNY 

JsN«L2 

JY«J^1 

DO  220  Kal»N 
DO  220  LajY,N 

B(J»KjaB(J»K)-A(J»L)YB(L*K) 


K.lNCTlt)l^i  S2(X) 

IF  (AHSlXJ,GT,l,fc-)b)  (lU  TP  lyo 
S2>U. 

RfTuWrs 

loo  S2«1*U»8InC(X)/X 
hfTURn 
EnP 

•• 


I. 


‘i 
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FU^tTlOiv  SilXj 

If  CAbS(XJ,GT,l  ,f«l5)  GO  r»f  too 

S3«0,o 

RtTuRN 

83*(Sli'iC(x)**2)/AbS{X) 

Rf  T URf\( 

END 


FUf^LTlOw  SXUX(KJ 

Sxuxsi  ,0 

U (4bS(O.LT,l,f-10J  HtTiiRM 
SX0X551NC  CX)/X 
Rt  T URw 
END 
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METRIC  SYSTEM 


BASE  UNITS: 


Quantity 

Unit_ 

SI  Symbol 

Formula 

length 

metre 

m 

mass 

kilogram 

kg 

time 

second 

s 

electric  current 

ampere 

A 

thermodynamic  temperature 

kelvin 

K 

amount  of  substance 

mole 

mol 

luminous  intensity 

candela 

cd 

SUPPLEMENTARY  UNITS: 

plane  angle 

radian 

rad 

solid  angle 

steradian 

sr 

DERIVED  UNITS: 

Acceleration 

metre  per  second  squared 

m/s 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration)/s 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

area 

square  metre 

m 

density 

kilogram  per  cubit,  metre 

kg/m 

electric  capacitance 

farad 

F 

A-slV 

electrical  conductance 

siemens 

S 

AIV 

electric  field  strength 

volt  per  metre 

V/m 

electric  inductance 

henry 

H 

V-s/A 

electric  potential  difference 

volt 

V 

W/A 

electric  resistance 

ohm 

V/A 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

I 

N-m 

entropy 

joule  per  kelvin 

)/K 

force 

newton 

N 

kg-m/s 

frequency 

hertz 

Hz 

(cycle)/a 

illuminance 

lux 

lx 

Im/m 

luminance 

candela  per  square  metre 

,, 

cd/m 

luminous  flux 

lumen 

tm 

cd-sr 

magnetic  field  strength 

ampere  per  metre 

Aim 

magnetic  flux 

weber 

Wb 

V-s 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  force 

ampere 

A 

power 

watt 

W 

J/s 

pressure 

pascal 

Pa 

N/m 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

)oule 

I 

N-m 

radiant  intensity 

watt  per  steradian 

W/sr 

specific  heat 

loule  per  kilogram-kelviii 

I/kg-K 

stress 

pascal 

Pa 

N/m 

thermal  conductivity 

watt  per  metre-kelviii 

W/m-K 

velocity 

metre  per  sec.onci 

m/s 

viscosity,  dynamic 

pascal-second 

Pa-s 

viscosity,  kinematic 

square  metre  per  sec;ond 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(wavo)/m 

work 

joule 

) 

N-m 

SI  PREFIXES: 


Multiplication  Factors 

Prefix 

SI  Symbol 

1 000  000  000  000  - 

10” 

tera 

T 

1 000  000  000  - 

10’ 

g(g* 

C 

1 000  000 

10* 

mega 

M 

1 000  = 

10’ 

kilo 

k 

too  = 

10’ 

hecto* 

h 

10  = 

10' 

deka* 

da 

0 1 =- 

10-> 

deci* 

d 

0.01 

10-’ 

centl* 

c 

(1001  -- 

10- ’ 

mllli 

m 

0 000  001  - 

10* 

micro 

M 

nooo(ioj(H)i  - 

10-* 

nano 

n 

0 0(M)  000  000  001  ^ 

10-  ” 

pi  CO 

P 

0.000  000  1 >0  000  001 

10 

fiimlo 

r 

0 000  000  000  000  000  001 

in-'" 

atto 

a 

* To  be  avoided  where  poaaibic 


